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Foreword 


After several friendly discussions of the pros and cons of tensors versus 
differential forms in the solution of engineering problems, I persuaded my 
colleague Dr. Flanders to prepare a number of lectures on differential forms. 
The result was an. outstanding series of lectures which was presented to a 
group of interested faculty members within the several schools of Engineering 
at Purdue University. 

It became obvious to those attending that the use of differential forms 
would give them another tool for the analysis and synthesis of engineering 
systems. There are certain problems, normally very difficult to solve by 
using tensors only, for which results are more quickly and directly obtained 
with differential forms. 

The author was encouraged to formalize his notes to the extent necessary 
for publication, to enable others to study this important subject. The text 
is recommended highly because differential forms and related concepts which 
have evolved from modern mathematics are new and powerful analytical 
tools for use by the engineer and scientist. 


GEORGE A. Hawkins, Dean 
Schools of Engineering and Mathematical Sciences 
Purdue University 


November 20, 1962 


Preface to the Dover Edition 


I have made the following changes to the 1963 edition. First, I have 
rewritten the proof of the Third Lie Theorem, starting on page 109, to better 
systematize the computations. Second, I have rewritten the derivation of 
some basic relations in phase space, pp. 164—165, to be less computational and 
more in the spirit of the book. Finally, there is an addendum to the bibliogra- 
phy on page 199 which should be useful. In particular, it mentions the 
forthcoming MAA Studies volume in which I give Kannai’s differential 
form proof of the Brouwer fixed point theorem. 


HARLEY FLANDERS 
dune 1989 


Preface to the First Edition 


Last spring the author gave a series of lectures on exterior differential 
forms to a group of faculty members and graduate students from the Purdue 
Engineering Schools. The material that was covered in these lectures is 
presented here in an expanded version. The book is aimed primarily at 
engineers and physical scientists in the hope of making available to them new 
tools of very great power in modern mathematics. Although none of our 
applications goes very deep, it is hoped, nevertheless, that enough ground 
is covered in each case to indicate the usefulness of this machinery. 

A word about the organization of the book is in order. The first chapter is 
introductory and sketches where we are going and why. Chapters II, III, 
and V include all of the theoretical material; a knowledge of this opens the 
door to the applications. Probably on first reading, one should aim more at 
developing some intuition for the subject and getting a firm idea of what the 
various different things which are defined look like, rather than at working 
out proofs in detail. Applications to questions in differential geometry (in- 
cluding many topics of considerable use in physical sciences) are mostly in 
Chapters IV, VI, VIII, and IX. Applications to various topics in ordinary 
and partial differential equations will be found in Chapter VII. Finally, 
applications to several topics in physics are in Sections 3.5, 4.6, 6.4, and 
Chapter X. 

What is presupposed of the reader is first of all a certain amount of scien- 
tific maturity, the precise direction not being too important. While the book 
is not really advanced mathematics, it is not exactly ground floor mathe- 
matics either, and a reasonable knowledge of the calculus of functions of 
several real variables is necessary, as is a working knowledge of linear algebra 
through the ideas of linear combination, basis, dimension, linear transforma- 
tion. Some exposure to a minimum amount of the ground rules of modern 
mathematics, sets, cartesian products, functions on sets, is helpful but not 
essential. This material is usually picked up by osmosis anyway, and the 
Glossary of Notation at the end of the book should be helpful. The reader 
should also know about the existence of solutions of ordinary differential 
equations. A passing familiarity with tensor methods is useful, but not 
essential. 

If our audience consisted of mathematicians alone, it would be in order to 
use somewhat more care in our formulations of definitions and proofs of 
theorems and to discuss in considerably more depth numerous technical 
points we here‘pass over lightly. Our goal, however, is to develop an intuition 
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and a working knowledge of the subject with as much dispatch as is possible. 
This perhaps could be done in less space except for our insistence on a degree 
of rigor matching that found in the better treatises on theoretical physics. 
This falls short of the extremely great precision which is customary in modern 
abstract mathematics and pretty much inherent in its nature. One who quite 
rightly is searching recent developments in mathematics for applicable 
materia] must find this precision a considerable barricade, overpedantic if not 
downright tedious—a very real factor in the great separation between modern 
mathematics and modern science. Making his craft available to science is not 
a light task for the mathematician and the extent to which this book makes a 
contribution therein must necessarily be its primary measure of success. 

In spite of all this, we do not hesitate to recommend this material to 
graduate students in mathematics as an introduction to modern differential 
geometry; indeed, a well-trained advanced undergraduate should find the 
book quite accessible. Considering the degree to which present day mathe- 
matical] training consists of one abstraction after another, some of the things 
in this book could be a bit of an eye-opener, even to a mathematics student 
who is well along. For example, one could envisage such a student meeting 
here a parabolic differential equation, or a matrix group, or a contact trans- 
formation for the very first time. 

It is my pleasant duty to acknowledge the substantial help and encourage- 
ment I have always had from my teachers, colleagues, and students. In this 
respect a special vote of thanks is due George A. Hawkins, Dean of the 
Schools of Engineering and Mathematical Sciences of Purdue University. 
Finally, I wish to express my gratitude to Elizabeth Young, whose beautiful 
typing of the manuscript was a substantial contribution. 


July 1963 HarLey FLANDERS 
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Introduction 


1.1. Exterior Differential Forms 

The objects which we shall study are called ezterior differential forms. 
These are the things which occur under integral signs. For example, a line 
integral 


| 4dr + Bay +04 


leads us to the one-form 
w= Adzx+ Bdy + Cdz; 
a surface integral 


[[ Pave: + Qdzdz + Rdzxdy 


leads us to the two-form 
a= Pdydz+Qdzdz+ Rdxdy; 


and a volume integral 
{Il H dx dy dz 


A= Hdxdydz. 


These are all examples of differential forms which live in the space E? of three 
variables. If we work in an 7-dimensional space, the quantity under the 
integral sign in an r-fold integral (integral over an r-dimensional variety) is an 
r-form in n variables. 

In the expression « above, we notice the absence of terms in dzdy, dzdz, 
dydx, which suggests symmetry or skew-symmetry. The further absence 
of terms dxdz,-:- strongly suggests the latter. 

We shall set up a calculus of differential forms which will have certain 
inner consistency properties, one of which is the rule for changing variables 
in a multiple integral. Our integrals are always oriented integrals, hence 
we never take absolute values of Jacobians. 


Consider 
|[4 (x, y)dudy 
1 


leads us to the three-form 
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with the change of variable 


" = 2(u, v) 


y = y(u, v). 
We have 





| | Ale, y)dedy = | | Afx(u, v), y(u, v)] = duds 


which leads us to write 





Ox Ox 

| Gu Ov 
dp ioe du dv. 

O(u, v) oy oy 

du dv 


If we set y= 2, the determinant has equal rows, hence vanishes. Also if 
we interchange x and y, the determinant changes sign. This motivates the 
rules 

dzdzx = 0 

dydzx = —dxdy 


for multiplication of differentials in our calculus. 
In general, an (exterior) r-form in n variables x', - ++ , x" will be an expression 


1 ; 
O= TVA. da‘! +++ dz", 


where the coefficients A are smooth functions of the variables and skew- 
symmetric in the indices. 

We shall associate with each r-form w an (r + 1)-form dw called the exterior 
derivative of w. Its definition will be given in such a way that validates the 


general Stokes’ formula 
| o= | dw. 
az x 


Here © is an (r + 1)-dimensional oriented variety and @X is its boundary. 
A basic relation is the Poincaré Lemma: 


d (dw) = 0. 


In all cases this reduces to the equality of mixed second partials. 


1.2. Comparison with Tensors 


At the outset we can assure our readers that we shall not do away with 
tensors by introducing differential forms. Tensors are here to stay; in a 
great many situations, particularly those dealing with symmetries, tensor 
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methods are very natural and effective. However, in many other situations 
the use of the exterior calculus, often combined with the method of moving 
frames of E. Cartan, leads to decisive results in a way which is very difficult 
with tensors alone. Sometimes a combination of techniques is in order. We 
list several points of contrast. 

(a) Tensor analysis per se seems to consist only of techniques for cal- 
culations with indexed quantities. It lacks a body of substantial or deep 
results established once and for all within the subject and then available for 
application. The exterior calculus does have such a body of results. 

If one takes a close look at Riemannian geometry as it is customarily 
developed by tensor methods one must seriously ask whether the geometric 
results cannot be obtained more cheaply by other machinery. 

(b) In classical tensor analysis, one never knows what is the range of 
applicability simply because one is never told what the space is. Everything 
seems to work in a coordinate patch, but we know this is inadequate for most 
applications. For example, if a particle is constrained to move on the 
sphere S*, a single coordinate system cannot describe its position space, let 
alone its phase or state spaces. 

This difficulty has been overcome in modern times by the theory of 
differentiable manifolds (varieties) which we discuss in Chapter V. 

(c) Tensor fields do not behave themselves under mappings. For 
example, given a contravariant vector field a’ on x-space and a mapping @ 
on z-space to y-space, there is no naturally induced field on the y-space. 
[Try the map t — (é?, t?) on E! into E?.] 

With exterior forms we have a really attractive situation in thisregard. If 


go: M— N 


and if w is a p-form on N, there is naturally induced a p-form ¢* q@ on M. 

Let us illustrate this for the simplest case in which q@ is a 0-form, or scalar, 
i.e., a real-valued function on N. Here ¢* w = w o d, the composition of the 
mapping @ followed by w. 


Mg 
4 


o*w=wod 


Reals 


(d) In tensor calculations the maze of indices often makes one lose sight 
of the very great differences between yarious types of quantities which can 
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be represented by tensors, for example, vectors tangent to a space, 
mappings between such vectors, geometric structures on the tangent spaces. 

(e) It is often quite difficult using tensor methods to discover the deeper 
invariants in geometric and physical situations, even the local ones. Using 
exterior forms, they seem to come naturally according to these principles: 

(i) All local geometric relations arise one way or another from the 
equality of mixed partials, i.e., Poincaré’s Lemma. 

(ii) Local invariants themselves usually appear as the result of applying 
exterior differentiation to everything in sight. 

(iii) Global relations arise from integration by parts, i.e., Stokes’ 
theorem. 

(iv) Existence problems which are not genuine partial differential 
equations (boundary value or Cauchy problems) generally are of the type of 
Frobenius—Cartan—Kahler system of exterior differential forms and can be 
reduced thereby to systems of ordinary equations. 

(f) In studying geometry by tensor methods, one is invariably restricted 
to the natural frames associated with a local coordinate system. Let us 
consider a Riemannian geometry as a case in point. This consists of a 
manifold in which a Euclidean geometry has been imposed in each of the 
tangent spaces. A natural frame leads to an oblique coordinate system in 
each tangent space. Now who in his right mind would study Euclidean 
geometry with oblique coordinates? Of course the orthonormal coordinate 
systems are the natural ones for Euclidean geometry, so they must be the 
correct ones for the much harder Riemannian geometry. We are led to 
introduce moving frames, a method which goes hand-in-glove with exterior 
forms. 

We conclude the case by stating our opinion, that exterior calculus is here 
to stay, that it will gradually replace tensor methods in numerous situations 
where it is the more natural tool, that it will find more and more applications 
because of its inner simplicity, body of substantial results begging for further 
use, and because it simply is there wherever integrals occur. There is 
generally a time lag of some fifty years between mathematical theories and 
their applications. The mathematicians H. Poincaré, E. Goursat, and FE. 
Cartan developed the exterior calculus in the early part of this century; in 
the last twenty years it has greatly contributed to the rebirth of differential 
geometry, now part of the mathematical main stream. Physicists are begin- 
ning to realize its usefulness; perhaps it will soon make its way into 
engineering. 


I 


Exterior Algebra 


2.1. The Space of p-Vectors 
Notation: 
R = field of real numbers a,b,c, ---. 
L. = an n-dimensional vector space over R with elements «, 8, ---. 


For each p = 0,1,2, --- , n we shall construct a new vector space 
A’t 
over R, called the space of p-vectors on L. We begin with 
K° L=R, N' L=L. 
Next we shall work out A’ L in some detail. This space consists of all sums 
>, aj(a; A B;) 

subject only to these constraints, or reduction rules, and no others: 

(41a, +4202) A B —a,(a, A B) — a2(a, A B) = 9, 

a A (6,8, + 6282) — 6,(a A By) — 62(0 A B2) =9, 


aAa=Q, 
anB+Bpaa=0O. 


Here a, 8, etc., are vectors in L and a, b, etc., are real numbers; « A f is called 
the exterior product of the vectors a and B. Ifa and f are dependent, say 
8B = ca, then 


an B=an (ca)=ctana)=cO0=0 


according to our reductions. Otherwise « A B # 0. 
Suppose a',---,o" is a basisof L. Then 


a=) a', B = > bo’, 
ar p=(> a') a (> 6,0’) = ¥ abo’ a o?). 
We rearrange this as follows. Each term a‘ A o'=0 and each o/ a ai = 
~o' voi fori<j. Hence 


ix} 
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The typical element of /A\? L is a linear combination of such exterior pro- 
ducts, hence the 2-vectors 


gag, lsi<jn, 


form a basis of N L. We conclude 


dim APL = MED. (*) 


In general, we form \? L (2s p< 7) by the same idea. It consists of 
all formal sums (p-vectors, or vectors of degree p) 


Yi ala, A +++ A Gp) 
subject only to these constraints: 
(i) (@a+ 5B) Aa, A ++ AG, 
= AKAM A *'* AG) + OBA A*** AO), 

and the same if any a; is replaced by a linear combination. 

(ii) a, A+++ Aa, =0 if for some pair of indices ¢ 4 j, a; = a;. 

(iii) a, A ++: A a, changes sign if any two a, are interchanged. 
It follows easily from (i) that «, A --: A «, is linear in each variable; we may 


replace any variable by a linear combination of any number (not just two) 
of other vectors and compute the value by distributing, for example 


a A (b,B, + 528, +6383) Ay AO 
=Bbi(an Bb, AYA) +5,(a AB, Ay AO) +4,(H A Bz Ay AO). 
It follows from (iii) that if x is any permutation of the {1, 2, ---, p}, then 
Oty Ao A Aap = (SEN Ma, A t's AG. 
Exactly as in the case p = 2, we can show that if 
at, a 


is a basis of L, then a basis of A? L is made up as follows: for each set of 
indices 

H = {h,,hz,-+-, hy}, 1Shy <hy<+++ <h, Sn, 
we set 


=a" A+++ Ao", 
Then the totality of o” is a basis of NV? L. We conclude that 


n 
dim A? L = ( ) 
A p 
the number of combinations of n things taken p at a time. In particular 


dim A\"L = |. 


2.2. DETERMINANTS i 
If Aisin A? L, then 


A = J ayo", 
H 


summed over all of these ordered sets H. One can also sum over all p-tuples 
of indices by introducing skew-symmetric coefficients: 


1 
A=— DP dae etp OE A TS AO 


! P 
Pp: Ai, eee, Ap 
where the 8,,...,, is a skew-symmetric tensor and 


b, -shp — 2H for H = {h,,-:-,h,}, hy <hay<-t+ <h 


sh cs ns 
This skew-symmetric representation is often quite useful. 

Let us note why we do not define A?L for p>. (Sometimes it is 
convenient to simply set KV? L=0 for p>n.) We express each « in a 
product a, A *** A a, as a linear combination of the basis vectors | a 


and completely distribute according to Rule (i). This leads to 
hp AS Ady = Da OPA TH AM, 


Each term a" a --- Ao"? is a product of p > n vectors taken from the set 
o',--+,o" so there must be a repetition; by Rule (ii) it vanishes. We are 
left with a, A :++ A a, =0 as the only possibility. 

We close with a very important property of the spaces /\? L. 

In order to define a linear mapping f on /\’L it suffices to present a 
function g of p variables on L such that (i) g is linear in each variable 
separately, (ii) g is alternating in the sense that g vanishes when two of its 
variables. are equal and g changes sign when two of its variables are inter- 
changed. Then 


Fay A +++ Ay) = G(%1,°7* » bp) 


defines f on the generators of A\’ L. 

It can be shown that this property provides an axiomatic characterization 
of h? L. In the next section we apply this property to define the deter- 
minant ofia linear transformation. 


2.2. Determinants 


As above L is a fixed linear space of dimension n. Let A be a linear 
transformation on L into itself. We define a function g = g, of n variables 
on Las follows: 


Ja(@1.°°*>%q) = Aa, A+++ A Aa, 


Ja: xK"L— A"e 
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where a L denotes the cartesian product. Since g is multilinear and 
alternating, there is a linear functional f = f, , 


‘vt NL—A’L 
satisfying 
faley N**" A G,) = G4(%,,°°*, &,) = Ag, A al A Ag, . 


But N’ L is one-dimensional so the only linear transformation on this space 
is multiplication by a scalar. We denote the particular one here by [4] 
and have 

Aa, A °°: A Aa, =|Al(a, A +++ A a,). 


This serves to define the determinant |A| of A. We must not fail to note 
that this definition is completely independent of a matrix representation 
of A. 
We observe next 
|ABl(a, A °°: A@,) = (ABa,) A +++ A (ABa,) 
= |A|(Ba, A --: A Ba,) 
= |A|-[|Bl(a, A +++ A O,), 


hence 
|AB| = |A|-| Bi. 
We can relate this to the determinant of a matrix as follows. Let 
o',---, 0" be a basis of L and |a,,|| an» x ~ matrix. Set 
x; = » a, ;07. 
Then 
a At AG, = |aylo' A --* Ao". 


In particular, if one obtains the matrix representation of A with respect to 
the basis (o') by 

Aa' = } a',o’, 
then Agin::: a Ao" =|aijol ass: Ao", [Al = at). 


2.3. Exterior Products 

We now observe that our spaces A’ L have a built-in multiplication process 
called exterior multiplication and denoted by A for obvious reasons. We 
multiply a p-vector p by a q-vector v to obtain a (p + q)-vector pz A v (which 
is 0 by definition if p+ q> 7): 


A: (APL) XK (AtL) — AP*tk. 


It suffices to define A on generators and use the basic principle at the end of 
Section 1 to extend it to all p- and q-vectors: 


(ap ATs AG) A(BL At AB) =A AGABLA** A B- 
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The basic properties of this exterior product are 

(1) AA pis distributive, 

(2) AN (HAV) =(AA Bp) A V, the associative law 

(3) pAA=(—I)PAA py. 
Property (3) simply says that any two vectors of odd degrees anticommute, 
otherwise vectors commute. The following will illustrate why this is the 


case. 
(4, AG, AGz,)AB= —(a, Aa, AB A4Q3) 


= (—1)?(a, AB Aa, A a3) 
= (—1)°B A (a, A & A m3); 
(%, A 2 A a3) A (By A B2) = (—1)°B, A (a, Aa Aas) A By 
= (—1)°(—1)7(By A Ba) A (ay A & A &3) 
= (—1)*7(B, A Bp) A (a A a2 A 3). 


Examples. We take for L the linear space based on the differentials 
dx, dy,::- and, as is customary, omit the exterior multiplication sign A 
between dx’s. Thus dxdy denotes dx a dy. 


|. (Adx+ Bdy + Cdz) a (Edz + Fdy + Gdz) 
= (BG — CF) dydz + (CE — AG) dzdz + (AF — BE)dzxdy, 


illustrating the vector-, or cross-product of two ordinary vectors. 


2. (Adx+ Bdy + Cdz) a (Pdydz + Qdzdzx + Rdzxdy) 
= (AP + BQ + CR) dxdydz, 


illustrating the dot-, or inner-product of vector algebra. 
3. Let « be any form of odd degree. Then 
ar7=arna=0. 
For if « and B are of odd degree p, then 


Baa=—aanB. 
We set B =a to have 


AAX=—-AAG, 2(a A «) = 0, avna=Q0. 


4. Here we take 
w = dp, dq' oe" + dp, dq’, 


@ form arising in mechanics. The two-forms dp,dq' all commute, hence 
a" = (n!) dp, dq’ dp,dq’ --- dp, dq" 
= (= 1)" Yat) dpy + dp, dg? ++ dg” 
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The product dp,---dg" is called phase-density. We shall discuss this 
further in Chapter X. 

We apply the exterior product to obtain the Laplace expansion of a 
determinant by complementary minors. 

Let ||a;;|| be an » x n matrix. For H = {h,,---,h,}, set 


Ain F1hp 
By hy eae a 


Setp+q=n. For K = {k,,--+,k,}, set 


Php 


By+t,ky  °  Uptiykg 
Cr — . 


a ere 


n,ky ny kg 


Thus if K = H’, the complementary set of indices to H (always arranged in 
natural order), then 6, and cx are complementary minors of |\a;,||. 
Now set 


_ i 
a,=) ao 
where (o/) is a basis of L. We easily see that 
a Att Aa, = > byo%, 
K 
Be, Att AG, = >) Cgar, 


hence 
Ap AAO = (hy Att A Oy) A (0544 A *t* A Oy) 

= V bycga" a aX. 

But 
OAc AG, = [ajl(o' A ++: AG") 

and 

0 if K # A’ 

tno =| 

eH (gl A +++ A a") if K = H', 

hence 


If H = {h,,---,h,}, H’ = {k,,-+>, k}, then 


] 2 e+ HL 
H,H’ __ 
sla EN | 


2.4. Linear Transformations 


In this section we deal with two linear spaces M and N with - 


dim M = m, dim N = n. 


2.4. LINEAR TRANSFORMATIONS ll 


Let us agree that when we need bases, a', --- , 0” will denote a basis of M and 
qi, ++, t" @ basis of N. 
Let A be a linear transformation, 
A: M—N. 
The mapping 
(&,,°°°,@,) —»> Aa, A+*+ A AG, 
sends 


x? M— I’ N. 
It is alternating multilinear, hence defines a linear transformation, denoted 
A? A on \’ M to NV? N. This exterior pth power of A is defined on 
generators by 


(A? A)(a, Pore kh a,) = Aa, AN “tt A Aa, . 
Suppose A is represented by the m x n matrix ||a' j|| according to 
Aa' = Yai ,w. 


The o” and c* form bases of (\? M and /\? N, respectively, where H and K 
are ordered sets of p indices. We have 


(A? A)o# = Ao a --+ A Ate 
ee yah, ee “ary, th A eee A he 
= yal tt. 
Hence A? A is represented by the matrix 
lor ell 
of all p x p minors of |ja‘,||. This is sometimes called the pth compound of 
oo one has three spaces L, M, N and this situation: 


.—_* »™ 


A 
AB 


We compute A? (AB): 
AP (AB)(a, A +++ AG) = (ABa,) A +++ A (A Ba,) 
= (A? A)[(Ba,) A «++ A (Ba,)] 
= (A? A)[(A? B)(a, A +++ A a,)] 
== [(A? A)(A? B)](a, A +++ A a), 


12 II. EXTERIOR ALGEBRA 


hence 
A? (AB) = (A? A)(A? B). 


It follows that the pth compound of the product of two matrices is the 
product of their pth compounds, a nontrivial result. 
We must consider one other matter. Again let d: M-—>N. Suppose 
@ isin A? M and 7 is in Ai M. Then 
(A?*4 A)(@ A 1) = (A? A)(@) A (At A)(y). 
For if we take monomials, @ =a, A ++: Aa,,n=B,A°*' A B,, then 
(A? *4 A) An) = (AP*4 AYO, A> Aa, A By A ++? A B,) 
=Aa,A+*: A AB, 
=(Aa, A ++: A Aa,) A (AB, A ++: A AB,) 
= (A? A)(w) A (At A)(n). 


2.5. Inner Product Spaces 


In the remainder of the chapter we shall study a space L which has an 
inner product («, 8). This is a real-valued function on L xX L. which is 
(i) Linear in each variable, 


(ii) Symmetric: (a, 8) = (f, «), 
(iii) Nondegenerate: if for fixed «, (a, 8) = 0 for all B, then « = 0. 


Example 1. The Huclidean inner product on E” is given by 

= (Ay, °°, Gq). = B = (Dy, +++ sy), 

(a, 8) = a,b, + «++ +4,5,. 

Example 2. The Lorentz inner product in four-space: 

= (4,,°°*, G4), B= (b,,-°*, 54), 

(a, B) =a,b, + azb, + 3b, — cab, 
where c is the speed of light. 
Condition (iii) is equivalent to the following. Ifa!,---, o" isa basis of L, 


then 
a’, o)| # 0. 


(The left-hand member is the Gram determinant, or Grammian.) For this 
determinant vanishes if and only if there is a nontrivial solution (a, ,--- , a,) 
of the homogeneous system 


¥ afa', a) = 0. 
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But this is the same as having the vector 
a=) a,o' 


satisfy the relation (a, 8) = 0 for all f. 
An orthonormal basis of L consists of a basis o!, --- , ¢” such that 


(o', o/) = + 6%. 
If there are r plus signs and s minus signs, then r+s=n, andt=r—<s is 
the signature of the inner product. It does not depend on the choice of basis. 
It is a basic fact that each inner product space L has an orthonormal basis. 
This is proved in several steps. 
1. If dimL > 0, there is a vector o in L such that (o, o) 4 0. 
For if (a, «) = 0 for all a, then 


O= (a+ B, a+ B) = (a, a) + 2(a, B) + (B, B) 
=2(a, 8), (a,B)=0 forall a, B, 
a contradiction to nondegeneracy. 
2. Pick a maximal sequence a’, ---, o” of vectors satisfying 
(a', of) = +64. 

Let M be the subspace of L these vectors span. Then dimM=r. [The 
a’ are independent since )' a,o' = 0 implies ) a,(a', o/) = 0, +a,;=0.] We 
suppose r <n. 

3. Let N be the orthogonal complement of M, i.e., N is the space of all 
vectors 8 such that (a, 8) = 0 for alla in M. Since N is determined by the 
r relations (a', 8) =0, dimN 2n—r. But obviously MAN =0 (ie., 
the only vector common to M and N is 0), hence dim N =  —7, M and N 
together span L,M + N = L. 

4. WN itself is an inner product space relative to the inner product of L. 
Only the property (iii) of nondegeneracy must be checked. Suppose f is in 
N and (y, 8) = 0 for all yin N. But (a, 8) =0 for all a in M, hence (a, £) = 0 
for all x in L since M and N together span L. Hence f = 0. 

5. By (1), there is a vector « in N such that (a, «) 40. We set 


of ** = a/|(a, a)|"/? 


and see that we have constructed a sequence a',-:-,o"*! longer than a 
maximal one. Since this is impossible we conclude that we must have had 
r =n in the first place, which completes the proof. 

There is another basic property of inner product spaces which we shall 
need below. 

Let f bé a linear functional on L. Then there is a unique vector B in L such 
that 


f(a) = (a, B). 
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This is easily established by taking an orthonormal basis o',---,0". We 
set 6, = f(c') and for 8 simply take 


B=> +b,0° =¥ (', o/b, 07. 


(0', B) = ¥ (6, ob, (0', 0) = b, =f (0"). 


For then 


2.6. Inner Products of p-Vectors 


Again we start with an n-dimensional vector space L with an inner product 
(a, B). We shali define an induced inner product on each of the spaces 


A’ kL. We set 
(A, nw) = (a; B;)I 


for A= a, A+** AG, H=B,A*+: AB,. This definition works because 
the determinant on the right is an alternating multilinear function of the 
a’s, ditto the B’s. This means the formula defines a scalar-valued function 
on (A? L) K (A? L) which is linear in each variable. Next (u, A) = (A, y) 
because interchanging the rows and columns of a matrix (transposing) does 
not change its determinant. 

The nondegeneracy of this inner product is most easily seen by comput- 
ing with respect to an orthonormal basis o!,---,o" of L. As usual the 
of, H = {hy <h, < ++: <h,}, form a basis of \? L. We have 


(o%, o*) = |(o", o%)]. 
If H # K, this is zero since the determinant has a row (also a column) of 
zeros. If H = K, all but the diagonal elements vanish and these are +1, 
hence 

(o", o*) = + bt K 
In other words; the oF form an orthonormal basis of Nd L, nondegeneracy 


follows free of charge. 
In particular ¢ = o' A --+ A o” is an orthonormal basis of A\” L and 


(c, 0) =(a', a) «++ (0", a") = (-1)"9? 


where ¢ is the signature of L. 
For another example, set 
1 i+ Atte an a", 


a=agian::-agd tag 


forming a basis of (\""' L. Clearly 
(a!, a!) = (a, 0)/(a', a') = (a, o)(0', 0’), 
(a,0', ¥ be) = (6, 0) ¥ (o', oa; 
= (6, 0)()) a, 0', )) bo"). 


hence 
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2.7. The Star Operator 


Again let L have inner product (a, 8). We shall take a definite orientation 
of L which will remain fixed. (This simply means we take one basis for L 
and only consider other bases which are expressed in terms of this one by a 
matrix with positive determinant. The space L has two orientations and 
we take one of them.) We only use bases coherent to the orientation. 

We shall define an operation *, called the (Hodge) star operator. This will 
be a linear transformation on /’ L. onto as L. This operator depends, 
of course, on the inner product and also depends on the orientation. Re- 
versing orientation will change its sign. 

We note that the orientation of L determines a definite orthonormal basis 
o of A’ L. 

Now fix / in \’ L. The mapping 


U— ANE 


is a linear transformation on /\"~? L into the one-dimensional space /A\” L. 
We may write 
Anp=f,(u)o 


where f, is a linear functional on yaad L. By our result at the end of 
Section 2.5, there is a unique (n — p)-vector, which we denote */ to indicate 
its dependence on A, such that 


AA p= (#A, po. 


This equation defines the *« map which is evidently linear on I’ L into 
APL. 

In order to compute */ for generators of \\? L, in view of the linearity, it 
is enough to compute *A where 4=a' A -*: A o? and where o!,---, o” is 
an orthonormal basis. Let K run over sets of g =” — p indices. Then 


An ok = (#1, o*)o. 

The left-hand side vanishes unless K = {p + 1, p + 2,---,n}, hence 

#A=coPT! A +++ Aa" 
and the constant c is determined taking K = {p+ 1,---,n}: 

o=ANG* = c(a*, ofa, 

c= (6%, o*) = +1, 

#1 = (a*, o*)o*. 

For definiteness, set H = {1,---,p},K = {p+1,---,n}. We have proved 


so! = (a*, o*)o*. 
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Since ok a og! = (—1)?"")a# an o*, we deduce, taking orientation into 


account, 
wok = (—1)PM "Pgh, gt)g!, 


hence 
* (ao) a (—1)P@— Pg, of\(ak, a*)o#, 


* (wo!) = (—1)?-?)(¢, a)o# 
= (—1)POrPt@-D/2gH 


where t 18 the signature. 
It follows that tf « 1s any p-vector, then 


xe = (— 1)P~ PIF O/2y 


Another consequence of these formulas is this result. 
If «, B are p-vectors, then 
aA «8B = BA ea = (—1)~%2 (a, Blo. 


For when Bf = c# as above, the only generator « = a” for which both sides 
do not vanish is « = o#, and then 


arn«B=a" a (6%, a*)ok = (6%, a*)a 
= (6%, o@)(-1)° 976 
= (—1)""?(a, Bo. 


Example 1. We take 4-space with coordinates so normalized that dz’, dz’, 
da*, dt is an orthonormal basis with (dz', dx’) =1, (dt, di) = —1. We have 
n=4, t=2, (—1)-9/? = —1. We shall study certain two-forms. For 
p=2,p(n-~-p)=4. Thus 


« (da! dt) = da dz" 
where (1, j, k) is cyclic order, 
* (dai da*) = —da' dt. 


Let E,; be the components of electric field strength, H, the components of 
magnetic field strength (all in free space) and consider the form 


@ = (E, dx! + E,dx* + E,dz*)dt + (A, da* dx’ + H,dz>dzx' + H,dzx' dx’). 
Then 
se = —(H,dx' + H, dx? + H,dx*)dt 
+ (E, dz* dx* + E,dz°dx'! + E;dz' dz’). 


We shall see the use of these forms in Maxwell’s equations later. 
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Example 2. E° with the ordinary metric. If f and g are functions, 


af = oF aa + ay + Fas oF a 
: of of 
= — dx d 
aaf Rp oye oe x dy 
and we have 
Of ag of og 22) 
d dx dy dz. 
Lae (72 dy oy ee 


2.8. Problems 

1. Let L be an n-dimensional space. For each p-vector «#0 we let M, be 
the subspace of L consisting of all vectors o satisfying a Ag=0. Prove 
that dim(M,) <y. Prove also that dim(M,)=p if and only if «= 


G, A ‘***: AG, where o,,°*:,6@, are vectors in L. 
2. (Continuation) Let « be any (n—1)-vector. Prove that a= 
oO; Att’ A On-1 


3. Let L be an n-dimensional space and « a 2-vector. Show that there 

is a basis o,, °°: , 06, of L such that 

a= Gd, AO2+ 63 A Gat bala + O4,-4 AN 62,- 
The number 27, which depends only on a, is called the rank of the 2-vector a. 
Show that a” # 0, a’*! =0. 

4. (Continuation) Let o,,-:-,o, bea basis of L and let A = |la,,|| be a 
skew-symmetric matrix. Show that the rank of the matrix A coincides 
with the rank of the 2-vector « = 4)’ a,,0, A ;. 

5. Weare given a linear transformation 


A: L—>L, 
Where dimL =n. Find the value of the determinant 
lA? Aj. 


6. Let A be an m x n matrix and B an n x m matrix, where m <n. 


Prove 


|ABl = > @uby 
where H runs over ordered sets, 
= {h,,ho,°°+, In}, lsh, <h,< +++ <h, Sn, 
and where 
Qin, 7 Bith,, Ont ean Ona 
ay = a , by = “i 
a b 


mM ,him fy, 
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Note that the special case 


a, @, a; 
A= ; B='A 
6b, 6, 


yields a well-known formula of vector algebra, 


ja x bj = [a|?|b|? — (a-b)?. 


7. Let A be ann x 7 matrix and denote by cof A the matrix of cofactors 
of A [so that A(cof A) = (cofA)A = |A|-Z]. Let by x be a typical element 


of A? (cof A). Show that 


bu K = +|(AlP- MOK 


where H’ is the set of indices complementary to H, ditto K’, K, and ay. x. 


is the corresponding element of 
A"? A. 


8. Express in terms of exterior algebra the formula from vector algebra 


a x (Bx y) = (a-y)B — (a-B)y. 


I] 


The Exterior Derivative 
3.1. Differential Forms 
Let P bea point in E”. The one-forms at P are the expressions 
}\a,dz', a, constants 
I 


These form an n-dimensional linear space L=L,. The p-forms at P are 
the elements of 


Ark = A? be. 


Y ay,dx"---da"?, ay constants. 


i.e., expressions 


Note that we are dropping the notation “A” so that differentials dz‘ juxta- 
posed will always be multiplied by exterior multiplication. 

Now let U denote an (open) domain in E”, A p-form on U is obtained by 
choosing at each point P of U a p-form at that point, and doing this smoothly. 
Thus a p-form w has the representation 


@ = ) ay(x',:-+, x")da™, 


where the functions a,(x) are smooth functions on U, differentiable as often 
as we please. 

The exterior algebra applies at each point of U and so may be interpreted 
on the differential forms on U itself. Thus if w is a p-form and 7 is a 
q-form on U, then w An is a (p + q)-form on U. (Of course w A 4 = 0 if 
P+q>n.) If 

@O= > a,x", n= > byda*, 
then 
wo An = > ayby dx" dx* 


So that the coefficients of w An are again smooth functions, being poly- 
nomials in the coefficients of w and yn. 
For example a one-form 


w= Pdr+Qdy+ Rdz 
19 
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may be identified with an ordinary vector field (P,Q, 2) in E°, a two-form 
a= Adydz+ Bdzdzx + Cdzdy 

may be identified with a polar vector field in E?. 


3.2. Exterior Derivatives 


We denote by 
F7(U) 


the totality of p-forms on U. In particular F°(U) is simply the set of all 
smooth functions on U. 
We shall now set up an operation d which takes each p-form w to a 
(p + 1)-formd@. In E? it will work this way. Fora 0-form f, 
of of of 


df = 4 ie 
fe A a 


For the one-form @ above, 


R 
dw = (= _- “Sy dy dz 


= | dede (= oP 
oy Oz 


Cee a aaa: 
bz Ox ax =) dedy 


while for the two-form « above, 


0A OB OC 
da = (= Tg: ao ay S:) ded 
Thus the operator d subsumes the ordinary gradient, curl or rotation, and 
divergence. 
It will turn out that d is completely independent of coordinate systems. 
This will be more or less clear when we axiomatize d. 
We shall establish the existence and uniqueness of an operator 


d: FU) —> FP*1(U) 
such that 
(i) d(a +n) =dw + dy 
Gi) dA a p)=di an pt (-1)E8Y) 2 A du 


(iii) For each w, d(dw) = 
(iv) For each function f, 


0 
af = Sas! 


Let us note the consistency of (iv) as it applies to the coordinate functions. 
For example z’ is a function on U and d(z') the effect of d on this function 
' is the symbol dz'. Thus from (iii), d(dz') = 0 once we have d. 
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First we prove there is only one such operation d. Suppose we are given 
such ad. We first show that 


d(da* ee - dz") a 


by induction on p. We have just noted this for p=1. If it is true for 
p — 1, then by (ii), 


d[a™ (da'2 - «+ dary] = dam --- date, 
d(da® +++ dat») = d{d(xdx™ --- dxt»)} = 0 
by (iii). Now if is a p-form, 
w= > a,(x)dx", 
dw = ) d(ay dz") 
= } (day) dx" 
=r oe t dal de®, 


which shows that the recipe (i-iv) completely determines dw. To prove that 
there exists such an operator d, we simply set 


do = SH db a 


for w = ) a,dz™ and check that the properties are satisfied. Properties 
(i) and (iv) are fairly clear; let us look at (ii) and (iii). Evidently if we can 
establish these for monomials, by summation they will follow generally. 
Suppose 
A=adz®, p=bdzx*. 
Then 
d(A A pt) = d(abda™ dz*) 


=) iad da' da™ da* 


da, i 0b ia eH gk 
= Lg, ode da™ da* +) a5 de dx™ dx 
-2(F di dx H) 1» (baek) 


a (—1)"e84) by (ada) A (= asda") 


= (dA) A pt (—1)9%”Y 1 A du. 
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The sign results from 
da’ dx = (—1)°®”) dat dz. 


This proves (il). 
Again, let @ =adz". Then 


da _. 
d oa ete i Bs | 
(dw) aly aa dz! dx ) 
Spi paca 
© Oa! dal - 


07a 07a 


1 | 
== ————— ‘dal 
3 (53 dat aa sai) de! de = 


I 
on 


} 


which verifies (iii). 

Property (iii) is nothing more than the equality of mixed second partial 
derivatives. It is the source of most “integrability conditions” in partial 
differential equations and differential geometry. It is usually referred to as 
the Poincaré Lemma. 


3.3. Mappings 
We study the following situation: U is a domain in E”, V is a domain in 
E” and ¢ is a smooth mapping on U into V. We write 


@: U-—-V. 
Also, we denote by z/,--+-, 2" the coordinates of E” and by y’, ---, y” the 


coordinates of E”. Then we can write 
yf = yllet,---, 2) 
to show that the point with coordinates x is transformed by ¢ to the point 


with coordinates y. The functions y'(x) are smooth. 
As before, R denotes the reals. Ifg is any real-valued function on V, 


g: V—-R 
then we may combine this with ¢ to obtain a function on U to R which we 
write 
o*9 = Goo. 
ot: F°(V) — FU). 


From the mapping ¢ on U to V we have constructed a new (induced) 
mapping ¢* on F°(V) to F°(U). 


Thus 
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o*g = Gog g 


R 


We are now going to define a map ¢* taking p-forms on V to p-forms on U: 
o*: F?(V) —> FP(U). 


(Strictly speaking we should index ¢* and write ¢,*, p = 0, 1,---+, but we 
shall skip this.) We have taken care of p = 0 already. The crucial case is 
p =1, after we do that, the algebraic considerations of Chapter II do the 


rest af the work. 
The basic idea is substitution of coordinate functions, replacing dy' by 


cy’ 
Se 


Thus if m = ) a,(y) dy’ is a one-form on V, we set 


7 dy? 
o* w => a,(y(x)) Agi da’. 
We now have 
g*: F'(V) — F'(U). 


By the method of Section 2.4, we extend this mapping to the exterior products 


to obtain 
g*: F'(V) — FU). 
As an example, 


p* (dy' dy”) = (¢* dy')(¢* dy’) 


(2%) eM 
i dy? 
da! Gx! 


1 2 
-55(45- wt) die da! 
dx! dai = Aa} Gar 


a! dae! 





_ ay’, y*) 4 i 
5D et ae dx. 
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We now list the basic properties of 6*. 
(i) p*(o+n) = GF wt o*n 
(ii) G*(A A p) = (G* A) A (G* H) 
(iii) If @ is a p-form on V, 
d(p* w) = p* (dw) 
(iv) If@: U—-+Vandy: V—W, then 


(Wo p)* = p* oy. 
The first property is evident and the second follows from the final formula 
of Section 2.4. 
Property (iii) is essentially the chain rule for partial derivatives. First 
we take a Q-form g on V. 


F F 
org =3 rn % oa 


5 << D) ax = ddrg. 








We proceed inductively, supposing we have verified (iii) for (p — 1)-forms. 
It suffices to verify (iii) for »-forms @ which are monomial since each p-form 
is a sum of such. Suppose then that 


w =gdy" = gdn 
where n = 7" dy" --- dy"? is a(p —1)-form. Then 
p* w = (p* g)(p* dn) = (G* 9) A (4 b* n), 
d(p* w) = d(h*g) A d(p*n), 


and 
dw = dg a dn, 


p* dw = (p* dg) A (p* dn) 
= d(p* 9) A d(p*n) =46* a; 
we have pushed through the next case. 


We now look at the final property (iv). 
For a 0-form (function) h on W we have 


[( o P)* AY(x) = AL(W o b)(x)] = hh Yld(x)]} 
= [P* A]G(x)] = {6* [p* A]} (x) 
= [(P* o W*)h](x), 
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hence 


(iso @) *h = (P* o W*)h. 


An induction similar to that above establishes the property in general. All 
it means is that one can substitute directly the expressions for the coordinates 


? ¢* 


U —_—____»- V F?(U ) —«¢—_—_———__ F*(V) 
wood Wy (Y o@)* y* 
= (p*o p*) 
Ww FW) 


z* on W in terms of the coordinates x' on U, or indirectly by first going 
through the coordinates y/ of V; the results are the same. 

What has really been seen in this section is that one can carry on fearlessly 
with the most obvious kind of calculations with differential forms. 


Examples. Consider the map ¢: t— > (x,y) on E'-—>+E? given by 
g=t?,y=t?. Ifw = xdy,aone-form on E’, 


)* w = (t”) 2 dt = 3t* dt. 


Take the map YW: (7%, y)—>t=2z—y. 
| y* (dt) = dx — dy. 
One finalremark. Suppose m < » and @ is a map on the domain U of E” 


into the domain V of E”. If w is a p-form on V and p > m, then necessarily 
o* w = 0. 


3.4. Change of Coordinates 


We apply the results of the last section to the special case in which U and 
V are both domains in E” and @ is a one-to-one mapping on U onto V with 
both ¢ and y = @~! smooth. (Note the map x — y= z° on E’ —> E' is 
one-to-one and smooth. But the inverse map y —> x = y'/? is not smooth— 
no derivative at y= 0.) In each figure, 1 is the identity map, i(x)= x. It 
follows that ¢* is a one-one map on F?(V) onto F?(U) and its inverse is y*. 
If we interpret the coordinates y of V as new coordinates on U, the result 


do* w = o*dw 


means that the exterior derivative of a differential form 1s independent of the 
coordinate system in which it is computed. 


26 Ill. THE EXTERIOR DERIVATIVE 


This inner consistency of the differential form calculus is most important. 
Later we shall base the global theory (forms on manifolds) on this. 


pee Vy ay 
; ) b 
U Vv 


We note in passing that with a proper formulation this independence of d 
on the coordinate system can be obtained as a consequence of the four basic 
defining properties (i-iv) of the exterior derivative in Section 3.2. 


3.5. An Example from Mechanics 

The following problem is taken from E. Goursat [15, p. 85]. We work in 
a region with coordinates (x, u) = (%,,°°',%,,%,;,°'*,4,). We are given 
a function 


which is supposed homogeneous of degree 2 in the variable u. (For example, 
a kinetic energy form )' a, ,(x)u,u;.) Define 


p= 0o/du,;. 


We assume that the mapping (x, u) —> (x, p) defines a regular change of 
variables. We then write 


$(x, U) = Y(X, Pp). 


The problem is to prove the relations 


ay ah Sy 


an, Ox; Op, 


= U,. 


The proof depends on two things, the Euler formula for homogeneous 
functions which in our case implies 


od 
y apy a D 
Ou, vk % 
i.e., 
> Pry 2¢, 


and the fact that exterior relations are independent of how they are derived. 


3.6. CONVERSE OF THE POINCARE LEMMA 27 


We have 


0d 0d od 
=) —dz, — =) —dz, d 
dp > dx, vj + La, du, > ox, ut; + > Pr Up 


and 
2db =) md + Duy aD,» 
hence by subtracting 


Now everything follows from ¢ = . and 


dy = ed Sd, +Ie ode. 


3.6. Converse of the Poincaré Lemma 
The Poincaré Lemma, d(dw) = 0, has these nee in 3-space: 


curl (grad f) = 
div (curl v) = 0 


according to the examples at the beginning of Section 3.2. In vector 
analysis one proves that a curl-free vector field is a gradient by line integrals 
and that a divergence-free vector field is a curl, usually by a brute-force 
method. We are now going to prove a general result. Jf w ts a p-form 
(p= 1) and dw = 0, then there 1s a (p — 1)-form a such that w=da. The 
result is hard if p > 1 because there are many solutions. Also the result is 
valid only in domains which are not too complicated topologically. 

The demonstration is based on a ‘cylinder construction.”’ We begin with 
a domain U in E”. We denote by I = (0, 1] the unit interval on the ¢-axis 
and consider the cylinder or product space. 


1X U. 

This consists of all pairs (¢, x) where 0 < ¢ < 1 and x runs over points of U. 

We single out the two maps which identify U with the top and bottom of 
the cylinder, namely, 


jr: UIXKY, file) = (1, x) 


jos U—IXK VY, — jo(x) = (0, x). 
Thus 

jit: FLX U)—> FU) (i = 0, 1). 
For example, to form j7,*w where w is a form on ix U. simplv replace 
t by 1 wherever it occurs in @ (and dt by 0 correspondingly). 
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We now form a new operation K, 
K: Ft x U)— FU); 
K is defined on monomials by the formulas 


K(a(t, x) dz") = 0 
1 
K(a(t, x) dtdx’) = (| a(t, x)at dz’, 
0 


and on general differential forms by summing the results on the monomial 
parts. Here is the basic property of K: If w ts any (p+ 1)-form on 
1 »4 U, then 


K(dw) + d(K@) =3,;* @ —jo* ow. 
It is enough to check this for monomials. 
Casel. w=alt, x)da*. 
We have Kw = 0, dKw =0, 


dw = = dt da" +[terms free of dé}, 


1 Oa 
Kdo = (| = a) dz = {a(1,x) — a(0, x)]dx™. 


0 
But j,* @ = a(1, x) dz", j,* @ = a(0, x) dx", so the formula is valid. 
Case 2. w= at, x) dtdz’. 


First j,*@ =jo*wa=0. Next, 
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Ca . 
x 


—_ (| Sorat) ata 


0 


1 
dKw =d i({ a(t, x) at) az" | 
o 


0 : 
=) re i a(t, x) at da;' da? 


1 a 


=X (| sre)artaer 


so the formula again works. 
Definition. A domain U is deformable to a point P if there is a mapping 


og: | xX U—U 
such that 
P(1, x) =x, 
o(0, x) = P. 
The boundary conditions may be interpreted in terms of the j, as follows: 
Poji=t Pojo=P. 


For a (p + 1)-form w on U we have as a consequence 
ji*lp*a)=o,  jo*o* wa] = 0. 
Now we can state and prove the main result. 
Let U be a domain in E" which can be deformed toa point P. Let w bea 
(p + 1)-form on U such thatdw=0. Then there is a p-form « on U such 


that 
aw = da. 


We merely substitute ¢* w in the formula above to have 
K[d($* w)] + d[K(o* w)] = o. 


But d(¢* w) = o* (dw) = 0, hence w = da with a = K(p*w). 

It is interesting to see how far the solution of the equation dx = w is 
determined. If B is another solution, then d8 = w = da, d(a— B)=0. If 
P 2 1, we conclude by the main result again that « — 8 = d/ where A is a 
(p —1)-form. In other words, given one solution a, the general solution is 
a — dA where / is absolutely arbitrary. (When p = 0, « and £ are functions 
and we conclude that « — B is constant.) 
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3.7. An Example 
We shall illustrate this whole method in the case n =3, p= 2. Thus we 


take a two-form 
wo = Adydz+ Bdzdzx + Cdxdy 


in E? for which dw = 0, i-e., 


Og OD oe 
Ox «Oy oz 


The space E* can be deformed to 0 by the map 
o(t, x, y, Zz) = (ta, ty, tz). 
The assertion is that wm = da where 


a= Kd* ow. 
First we compute ¢*q@: 


o* w = A(ta, ty, tz) d(ty) d(tz) +°-- 
= A(ta, ty, tz)(tdy + ydt)(tdz + zdt)+--- 
= A(ta, ty, tz)(ytdtdz — ztdtdy) + --+ + (terms free of df). 


Now we have 


1 
a= K(p*w) = (| A (tz, ty, te) td) yd — zay) 
0 


1 
re (| Bite, ty, t2) tat) (zdx — edz) 
0 


1 
+ (| Cita, ty, 2) (xdy — ydz). 
0 


One verifies after some calculation that indeed da = w. 


3.8. Further Remarks 


For 
w= Adydz + Bdzdx+Cdzxdy 
the problem of finding 
a= Pdx+Qdy+ Rdz 
so that 


da =a 


is that of finding three unknown functions P,Q, R of the three variables 
x,y, so that the system 
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oR @6Q 
eee aa 
oy = 
oP OR 
ae” 
La, 
Ox oy 


of three partial differential equations is satisfied, the given functions A, B,C 
being subject to the necessary condition 
0A OB OC 0 
re + Oy + Fleas 
It is remarkable that this system (and the more general ones covered in 
Section 3.6) can be solved by an explicit formula involving quadratures. In 
general, the theory of exterior differential forms exposes many types of 
systems of partial differential equations which are reducible to systems of 
ordinary differential equations and often solved by quadratures. 

Another point to be noted is this. If we are dealing with a (p +1)-form w 
such that dw = 0 and w happens to depend on several parameters smoothly, 
then we can find an « such that dx =w and « depends on the same parameters 
justassmoothly. This again follows from the explicit formulas of Section 3.6. 


3.9. Problems 


1. Let 
o=t) a,,dx' dz’, a,;+a,;; = 9. 
Prove that 
Ga,, Oa, O&:\,.,, 
dw = — + 4 | dx! dxi dx’. 
ss + (= : Gx +5 de 
2. Consider a linear transformation ¢: E”"—>E", ¢(z',-::,2")= 


(y',--+, y"), where 
y' =) a',zi+b', — a',,b' constant. 
What is o* (dz! --- dx")? 
3. Consider the mapping 
p: (z,y) — (xy, 1) 


on E* into E?. Compute $* (dz), o* (dy), and o* (ydz). 
4. Complete the unfinished calculation of Section 3.7. 


IV 


Applications 


4.1. Moving Frames in E° 


We first point out that in dealing with vectors in Euclidean space, no 
matter where we draw them for picturesque purposes, when we deal with 
them analytically, they always start at the origin. 


z 





x 


We attach to each point x of E* a right-handed orthonormal frame 
@,, €,, €, and suppose that the vector fields e; are smooth fields. 

What we shall do is express everything in sight in terms of the e,, apply d 
to these relations to derive further ones, and continue until we obtain no 
further results. 

First of all, dx is a vector with one-form coefficients, for example, dx = 
(dx, dy, dz) =dxi+dyj+dzk. We express dx in terms of the frame 
e,, @,, @; at the point x, which we certainly may do, say, by first expanding 
i, J, k in terms of the e; and then collecting terms: 


dx = 0,e, + 0,€, + 03e;, 


32 
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where the g; are one-forms. We do the same with each e;: 
de; = WiC, + Wir. + W303 (2 = 1, 2,3) 


where the w,,; are one-forms. 
Since e,°€, = 0;,, we have 
de,:e, + e, de, = (), 
that is, 
On, + O,; = 9. 
In particular, @;; =0. 
It will be convenient to introduce some matrix notation. We set 


e; 
e = [e, f, g¢= (oO; » 02, a3), Q = Ilo, ll 
e; 


and have these structure equations: 


dx = oe, 
de = Oe, 
Q+/'Q=0. 


Here applying d to a matrix means simply applying it to each element. In 
the last equation, the left-hand superscript ¢ denotes transpose of the matrix, 
ie., interchange of rows and columns, so this equation expresses the skew- 
symmetry of 0. 

From d(dx) = 0 we have 


doe —ode=(Q, 
doe — oQe = 0, 
(do — oQ)e = 0. 


Because the e; are linearly independent, this means 
do = 02. 
Similarly, from d(de) =0, we have 
0 = dQe — Qde = (dQ — N*)Je, 


dQ =’. 
In summary, then we have 
Structure equations Integrability conditions 
dx = ox de =aQ 
de = Qe ms = | 


QO +'Q =0 
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Further differentiation does not lead to new results. We shall see in our 
study of Riemannian geometry that the equation dQ — Q? = 0 expresses the 
lack of curvature of Euclidean space. 

A point to be noticed is that the three-form o, A 0, A @; is precisely the 
element of volume in E?: 


6, N06, A063 = dxrdydz. 
We shall verify this in the next section. 


It will be observed that the calculations of this section work equally well 
in E”. 





x 


Example. Spherical coordinates. The orthonormal unit vectors e,, e,, e3 
are taken in the directions of increasing r, ¢, 6, respectively. From 


x = (rsin ¢cos@, rsing sin 8, rcos @) 
we have 
dx = (sin ¢ cos 8, sin ¢ sin 8, cos d) dr 


+ (rcos cos 0, rcos@sin#, —rsin gd) dd 
+ (—rsin dsin 0, rsin ¢ cos 8, 0) dé 
= (dr)e, + (rdp)e, + (rsing dé)e, 
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with 
e, = (sind cos @, sin d sin 8, cos ) 
e, = (cos¢cos§, cosdsind, —sin d) 
e, = (—sin 8, cos 8, 0) 
and so 
o, = dr, o, =rdd, g,=rsingdé. 
Differentiating, 


de, = (dd)e, + (sin¢@dé)e, 
de, = (—dd)e, + (cosddé)e, 
hence since Q is skew-symmetric, 
0 dd sin @ d@ 
Q = —dd 0 cos @ dé 
—singd# —cosddé 0 
The volume element is 


a, A060, A0;,=9r'singdrdd dé. 


4.2. Relation between Orthogonal and Skew-symmetric Matrices 


It is no accident that Q turns out to be skew-symmetric. This is a con- 
sequence of the principle that the first-order approximation to an orthogonal 
transformation is a skew-symmetric one. We shall look at this from several 
viewpoints. 

A matrix B is orthogonal if its transpose equals its inverse, ‘B= B™', or 
B'B='BB=1I1. Suppose A is skew-symmetric, A+‘4 =0. Then for 
small ¢ we set B= I + €A and have 


B‘B = (14+ eA)(I — eA) =I + O(e) 


so that B is orthogonal up to first-order terms. 

Here is another approach. Let A be skew-symmetric. Since the 
characteristic roots of A are pure imaginary, J+ A and J — A are non- 
Singular. Set 


Then 


So that B is orthogonal. 
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Next we re-examine the calculations of the last section. Let 
i, 
i=] i, 
i, 
where the i j are the fixed unit vectors in the z, y, z directions, respectively 
(i, j, k in usual vector notation). Then 
e; = >: 5, ji;, e= Bi 
leading to a matrix B = ||b,,|| which is clearly orthogonal: 
I=e'e= Bi'i'B = BI'B= B'B. 
(Now we can prove the fact drdydz =a, A a, A @3 mentioned at the 
end of the last section. We have 
dx = (dz, dy, dz)i = ce = o Bi, 
(dz, dy, dz) = o B, 
hence 
dzdydz=|Bla, Ag, A@3. 
But from ‘B B = I we have |B|? = 1, |B| = +1. Since we are supposing e 
is a right-handed system, |B| = +1, 
dudydz=6,; N0,A63;) 


Then we have 
de=dBi= (dB) B™'e 


so that 
Q = (dB) Bo. 
We note this general result: If A is an orthogonal matrix whose elements are 
functions of any number of variables, then 
(dA) A7! 


is a skew-symmetric matrix of one-forms. 


For we have 
"'AA=T, 


'dAA+‘AdA =0, 
'‘A~1"d44+dAA™'=0, 
'dAA~')+dAA™'=0. 

There is also a converse which is important. Suppose A ts a matrix of 


functions defined on a domain U. Suppose A 1s orthogonal at a single point 


of U and that 
dA=AA 
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where A is a skew-symmetric matrix of one-forms. Then A is orthogonal on 


all of U. 
We set C = ‘A A and have 


dC = (dA) A +'A(dA) = (—'A A) A +'A(AA) = 0, 


hence C is a constant matrix onU. But we are assuming C = J at one point 
of U, hence C = I on U,'A A =TJ on U, A is orthogonal. 

Another point is this. If A is a variable orthogonal matrix (transforma- 
tion), each point v, of space is sent by the general A to 


v= AV. 
We then have 
dv =dAv, = (dA) A'v 


so that one passes from v to the “infinitely near’ vector v + dv under the 
action of the general A of our family by means of 


v—v+dv={I+(dA)A™'*]v 


with the skew-symmetric (dA) A™~! representing this ‘infinitesimal trans- 
formation.” 
All of these considerations work equally well in E”. 


4.3. The 6-dimensional Frame Space 


We consider the space of all right-handed orthonormal frames E,, E,, E; 
at all points x of E°. This space is 6-dimensional because we have three 
degrees of freedom in choosing x, two degrees of freedom in choosing the 
unit vector E,, one degree of freedom in choosing the unit vector E, per- 
pendicular to E, and then E, is determined. 

We write 


E, 
E={E, 
E, 
and have 
E=-Ae 


where A is a variable (three parameter) orthogonal matrix and e = e(x) is a 
definite moving frame. 


Then 
dx =ce=oA~'E, 
dE = (dA)e + Ade =[dA + AQ]e 
= {dA + AQ]A™~'E. 
We set 


a¢=0A', = (dA)A~1+ ANA}, 


38 Iv. APPLICATIONS 


These are matrices of one-forms on the 6-dimensional frame space and we 
have 


Structure equations Integrability conditions 
dx = cE dé = 6 
de = GE = 2 “4 
O + 'Q=0 


To check the integrability conditions we note 
0 = d(dx) = deE — sdE = (da — 6Q)E, de =G0, ete. 


In making a penetrating study of the differential geometry of E° one is 
necessarily led to this 6-dimensional frame space and its differential forms 
&', &,; which, it will be noted, are entirely independent of the choice of the 
moving frame e on E?. 


4.4. The Laplacian, Orthogonal Coordinates 

We continue the considerations of Sections 4.1 and 4.2. The forms 
dx, dy, dz make up an orthonormal basis for the Euclidean geometry of the 
space of one-forms at each point; these are related to the fixed (absolute) 
frame i. From 


e= Bi, dx= oe = (dz, dy, dz)i 
we have 
= (dz, dy, dz) 


as already noted. As Bis orthogonal, we see that ¢,, ¢,, 03 18 an orthonormal 
basis for one-forms at each point. 
Let f be a function on E*. Then we have 


of of of 


df=—d dy dz, 

a ae rac eer 
sdf = oF ay de + ot dds + F aedy, 
af af af 


dadf = (= zm ae 2+ 53 7) de adydz = (Af)dxdydz. 


The Laplacian Af of f is known as soon as the three-form d «df is known, for 
this has turned out to be the Laplacian multiplied by the volume element 
dix dy dz. 
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Now we know that the * operator can be computed equally well in any 
orthonormal coordinate system. Also o, Ao, Ao; =dzxdydz, so our 
procedure is this. We express df in terms of the a,, 


df =a,0, +a,0, +4303. 
Then 
#df = 0, 0,03,+4,030, + 436,02, 
dadf = (Af)o,0,03. 


A coordinate system u,v,w in a domain in E® is called an orthogonal 
coordinate system if the vectors 


ox Ox Ox 


—— 


au’ dv’ dw 
are mutually perpendicular. This means that for suitable functions A, LL, Y, 
the vectors 
7 l ox i 1 ox “ 1 ox 
tT) Ou’ 2p bv’ a) ow 


form an orthonormal, or moving frame. We shall presuppose that this is a 
right-handed one. (Otherwise we merely permute w and v.) We have 


ox ex ox 
cary, pia == 
dx gC Lie - + dw a 


= (Adu)e, + (udv)e, + (vdw)e, 
so that 


a, =Adu, o, = pdv, o, = vdw 
build an orthonormal frame for one-forms. Now we compute the Laplacian: 
af = f,du + f,dv +f, dw 
= (fula)or + (folk) O2 + (ful ¥) os - 
df = (f,/A)o203 + (film) oxo; + (flv) 0; 02 
= (pv fA) dudw + (Av f,/p) dwdu + (Au f,/v) dud», 
We compare this to 


dadf = (Af)o,0,0, = ApwAf)dudvdw: 


nue = E (7 iu) oe a3 =e * Be oe i) } 
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Let us apply this to spherical coordinates r, ¢, @: 
x = rsin @ cos@ 


y= rsin gd sind 


z=rcos@. 





x 


The orthogonality is easily checked (it is obvious geometrically) and we have, 


oa, = dr, o, =rdd, o,=rsingdé, 


0 — sang Lar (8# a) + 39 (6 ap) * a5 (ga) 


4.5. Surfaces 


We study a smooth surface £ in E>. We choose a moving frame e at each 
point x of Zin such a way that e, is the normal to the surface. Then e, 
and e, span the tangent plane at each point. We shall see how the equations 
of Section 4.1 specialize. 
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Since x is constrained to move in the surface, dx must lie in the tangent 


plane, 63 = 0: 
ax = 0,@, + 02€,. 


It is clear that the two-form 0,0, represents the element of area of E. 


e; 






unit normal 


tangent plane 


e, 


We exploit the skew-symmetry of Q by writing 


0 @ —@, 
Q=f-w 09 —-a, 


O, @, 90 


The structure and integrability conditions now reduce to 


Structure equations Integrability conditions 
aK = 0,€; + O,€, do, = 00, 
de, = we, — ae, do, = —wo0, 
de, = —@we, — @2e; 6,0, + 0,0, =9 
de, = wW,e, + @2e, dw + W,0, =9 
dw, = UM, 


dw, —— Ww, 
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In a certain sense, all of local surface theory is contained in these equations, 
It remains to interpret them in terms of curvatures, curves on the surface, 
etc. We illustrate a little of this. 

As already remarked, o,¢, is the element of area on XZ. As X moves over 
x, @,; moves over a region on the unit sphere $7, called the normal, or 
spherical, vmage of &. Since e, and e, are orthogonal to e,, they lie in the 
tangent plane to the spherical image and form a frame there. We see that 
the equation de, = m,e, + we, plays the same réle for the spherical image 
as dX = ¢,e, + o,e, does for XZ, hence w,w, represents the element of area 
of the spherical image. 

Since there is only one linearly independent 2-form on the 2-dimensional 
space L, we have 

10, = Ko,o, 


where K is a scalar called the Gaussian curvature. We shall see shortly that 
it is entirely independent of the choice of e, and e,. 
Similarly ¢,@, — ¢,@, is a 2-form on %, and so 


01@> — O2, = 2Ho,¢, 


defines a scalar H called the mean curvature of Z. 
The one-forms @,, @, are linear combinations of ¢, anda,. Because of 
the relation 
010, + 0,0, =0 


we have a symmetry in the coefficients: 
O, = po, + qo 
M2 = 90; + 102. 
We easily have from this 
2H=p+r, K=pr—q’. 
The characteristic roots of the symmetric matrix 
Pq 
a 
are called the principal curvatures K,,K, of Z. We consequently have 
2H =k, + Kp, K =k,k,. 
From the relation dw + @,w, = 0 we have 
do + Ko,0,=09. 
This relation gives us K once we know o,,¢, anda. But the relations 


do, = Wd, ; do, = —TWO; 
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suffice to determine m once o, and g, aregiven. (For then do, = ao,0, and 
da, = bo, are determined and we must have w=ao,+06c,.) In total 
then, K is completely determined analytically from o, and o,. This contains 
the theorem of Gauss that the curvature XK is an intrinsic invariant of &, 
independent of how = is imbedded in E°, so long as the distance between 
points of & measured along Z (on geodesics, or shortest paths) is preserved 


locally. 
When we apply vector operations to vectors with differential form 


coefficients, we must always combine the coefficients according to the rules 
of exterior algebra and pay strict attention to the ordering of the factors. 
With this we form vector (cross) products: 
= 6,7(e, x e;) + 0,7(e, x ez) + 0,0,(€; x €2) + 020,(€2 x e;). 
Now o,7 = 0 (and e, x e, = 0), ete. Also 
O20,(€, X €,) = (—o,02)(—e, X e,) 
= (0102) e3, 
so finally 
ax x dx = 2(0 102) e3 


and we have obtained the vectorial area element. 
Precisely, the vectorial area element is 
(5,02) @3, 
a vector directed along the normal with magnitude ¢,¢,, the element of 
area of £. Since 
dx x dx = (dx, dy, dz) x (dx, dy, dz) 
= 2(dydz, dz dz, dx dy) 


we have 
(dy dz, dzdx, dxdy) = (0,0,)e;. 


If v = (P,Q, R) is a vector field, then 


| (Pdydz + Qdzdz+ Rdxdy) = | ¥-(o.eses 
z 


= | (v-e3)(0,0) 
r 


is the flux of v through =. 
Similarly we have 
dx x dx = 2(0,0,)e; 
ax x de, = 2H(¢,02) eC; 
de, x de, — 2K (0,0) e; 
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which shows the independence of H and K on the tangent vectors e, , e, . 
If f is a function on Z with 
df = a,0; +2504 ; 
then on 2, 
sdf = —4,0;, ++- a,F4 ; 
dxdf = d(—a,0, + 4,02) = (Af)o,o, 
defines the Laplacian of f on the surface or the second Beltrami operator A. 
The same works for vectors and we have 
dx = a,e; + Ge, 3 
xX — 0,e, — G,e@, : 
We notice that 
= ge, — a; e€, ; 
hence 
dx = dx Xx @;, 


d«dx = —dx x de, = —2H(o,0,)e, 
and so 
Ax = (Az, Ay, Az) = —2He,. 


A minimal surface (surface of stationary area) is one for which the mean 
curvature vanishes, H=0. We have proved: The coordinate functions 
z,y,z are harmonic on each minimal surface. (That is, they satisfy Az = 
Ay = Az =0.) 

In this section we have given a sample of how the exterior calculus fits 
into the classical differential geometry of surfaces. Further materia] will be 
found in Sections 8.1 and 8.2, but there is much of the subject that we cannot 
cover in this text. A treatment from this point of view of exterior calculus 
which is not quite completely satisfactory and which unfortunately is em- 
bellished with historical comments often in bad taste is found in Blaschke [3]. 


4.6. Maxwell’s Field Equations 


In classical electromagnetic field theory one deals with the following 
quantities: 


E = electric field H = magnetic field 
B = magnetic induction J = electric current density 


D = dielectric displacement p = charge density. 
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These are all functions of the space variables x', x*, x? and the time t. The 
basic Maxwell equations in ordinary vector language are 


1 
(i) curlE = —- cha (Faraday’s law of induction) 
c at 
1 oD 
(ii) curlH = — = J + a ar (Ampere’s law) 
(iii) divD = 4xp (continuity) 
(iv) divB =0 (nonexistence of true magnetism) 


Here c is the speed of light. We shall put these equations into the language 
of exterior forms. To this end, we set 


a = (EH, dz! + H,dx* + E,dz*)(cdt) 
+ (B, dz? dz? + B, dz? dz! + B,dx' dz’), 
B = —(H, de! + H, dx? + H,dx3)(cdt) 
+ (D, dx? dx® + D,dzx3 dz! + D,dz! dx’), 
y = (J, dz? dx3 + J, dz? dz! + J,dx' dz?) dt — pda dx? dx’. 
Equations (i) and (iv) become 


da = 0. 
Equations (ii) and (iii) become 


dp ++ 4ny = 0. 
Applying d to this last equation yields 


dy = 0, 
in vector notation 


0 
div J + = = 0. 
From the equation dx = 0 one concludes, at least in any region of space-time 
which can be shrunken to a point, that there is a one-form A such that 
dA =a. 
We introduce the vector potential A and a scalar A, by writing 
A= A,dzx! + A,dz? + A,dx? + Agedt. 

The equation dj = a in vector form is 

curlA = B 


1dA 
grad 20 say 


46 IV. APPLICATIONS 
In free space, everything simplifies according to 
E=D, H = B, 
J = 0, p=0 
so that the Maxwell equations become 
curl E = — = divE = 0 


1 
Pt hla div H = 0. 
c ot 


We introduce the Lorentz metric into 4-space whereby 


dx', dx*, dx, cdt 
is an orthonormal! basis: 


(dzx',dxi)=65, (dz', cdt) =0, 
(cat, cdt) = —1. 


The signature is 3 — 1 = 2. 
According to the formulas of Section 2.7, 


«(dx! dz”) = —dz3(cdt), etc., 


*(dx' cdt) = dx? dx’, etc. 
We see that 


a = (EB, da! + -+-)(cdt) + (H, dx? dx? + ---), 
B = —(H, dx! + -+-)(cdt) + (E, dx? dx? +--+) 
= *O. 
Consequently Maxwell’s equations in free space are simply 
da = 0 
dxg = 0. 


We return to the general situation and refine our analysis by introducing 
one-forms: 


w, = E,dz' + E,dzx* + E,dx° 

, = B, dz? dz? + B, dx dz' + B,dx' dx* 
wo, = H,dx' + H,dx* + H,dx° 

w, = D, dx* dz? + D,dzxi dz' + D,dzx' dz? 
wo, =J,dx*dx° +J,dx° dx! + J,dz' dz’. 
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These involve space variable differentials only. Now we interpret d’ to 
denote the exterior derivative with respect to space variables only. We intro- 
duce 0/ét in this form 


0 : 
5 (@1) = Oy = H, dx + ---, ete. 


Now the Maxwell equations are 


d‘'o, = —-@ 

1 C 2 
Bi L 
Ve ee 
d'w, = 0 


d'w, = 4npdz' dx? dx’. 


The Poynting energy-flux vector S is introduced by 
S = (=) Ex H 
4n 
c 


(=) @, Aw; =S, dx? dx? + 8, dx> dx! + 8, dx! dx’. 
tt 


Poynting’s theorem, 


that is 


(=) B-H+E-J+ (=) E-D + divS =0, 


follows from 


a'(w, A@;3) = @'a,A@;-0,Ad@, 
Ll. (" ear 
—-@,)/A@,.—0@,A[—@ -@ 
Ae 3 1 Pi 


1, 4n ] ; 
a a Es Ch ORE ak 
Cc Cc 


For bodies at rest, one assumes D = cE, B = »H where the dielectric constant 
kK and the permeability yp are constant in time. Then Poynting’s theorem 
becomes 


Ou 
——=div8$+E:- 
ee 
where 


1 
u = — (KE? + pH?) 
8x 
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is the energy density of the field. The quantity E-J is called the therma, 
chemical activity. 


4.7. Problems 

1. Develop the formula for the Laplacian in cylindrical coordinates. 

2. A complex matrix A is unitary if A A*=J, where A* ='A, the 
transpose conjugate of A. We call A skew-hermitian if A4*+ 4 =0, 
Discuss the connection between unitary and skew-hermitian matrices. 

3. Show that e4 is orthogonal if A is skew-symmetric. Here 


, ir a) nt 
eto=] + py yl 
for the real matrix A. 

4, Set up a frame and the structure equations for a sphere of radius R 
Compute the curvatures. 
5. Find Gaussian curvature of the surface of revolution obtained by 


revolving the curve 
x = cos@ + In tan (6/2) 


y = sing 
5 <O<n 


about the x-axis. 

6. Given a surface in the form z = f(z, y), develop formulas for H and K 
in terms of f and its partial derivatives. 

7. Let & bea surface. Let e,,e, and e}, e, be two moving frames o: 
tangent vectors to X£. Determine the relation between the corresponding 
@ and a’ and verify that dw = da’. 

8. Lete,,e,,e, and e/, e5, e, be two moving frames in E°. Set ut 
the orthogonal matrix relating these frames and determine how the corre. 
sponding Q and 1’ are related. 


Manifolds and Integration 


5.1. Introduction 


An n-dimensional manifold is a space which is not necessarily a Euclidean 
space nor is it a domain in a Euclidean space, but which, from the viewpoint 
of a short-sighted observer living in the space, looks just like such a domain 
of Euclidean space. A case in point is the two-sphere $*. This cannot be 
considered a part of the Euclidean plane E*. However our observer on $2 
sees that he can describe his immediate vicinity by two coordinates and so 
he fails to distinguish between this and a small domain on E?. 

We have the technical problem of describing an »-manifold with sufficient 
precision so that we can define functions, tensors, and differential forms on 
such a space. The definition which follows is motivated in this way. Each 
observer on the manifold has an immediate neighborhood (local coordinate 
neighborhood) described by » coordinates. Each point of the space must 
lie in at least one of these observed neighborhoods. Now if we consider 
simultaneously two observers, their immediate neighborhoods may overlap, 
and we must specify what happens in each such overlap. In the next three 
sections we go over these matters with some care. 

After this is accomplished we tackle the problem of defining the integral 
of a differential form. In Sections 5 and 6 we lay the groundwork by 
defining chains, the geometrical sets over which forms are integrated, and 
in Section 7 we define the integral. 


5.2. Manifolds 


An n-dimensional manifold consists of a space M together with a collection 
of local coordinate neighborhoods U,, U,,°-- such that each point of M lies 
In at least one of these U. On each U is given a coordinate system 


lees ot 
So that the values of the coordinates 
(a'(P), oes © , x"(P)), 


Where P ranges over U, make up an open domain in Euclidean n-space E”. 
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Suppose that U with coordinate system 
zx}, vee gn 
and V with coordinate system 
1 n 
y , eee , y 
overlap (intersect). We may express the V coordinates y of a point P in 
terms of the U coordinates x of this point: 
yi = yal, ---, 2") (i= 1,--+,7n). 


As part of the definition, we assume that these functions are smooth (differen. 
tiable as often as we please). 


Having this formal definition out of the way, we explore some consequences 
First of all, on the overlap of U and V above we may interchange the rdéle 
of U and V to write smooth functions 


aay, --,y") (j= 1,+++,n). 
Substituting yields 
y= yi (z'(y), *-, ay) 
and we may differentiate by the chain rule: 
dy! Oxi 
k= Lae aye 
which has the matrix interpretation 
oy |). |" I 
Oa!) dy*|) 
We take determinants by the product rule: 
ay’, ree y") O(z', vee ah) 
A(z), +++, 2") ; ay’, oe y") 
It follows that the Jacobian 




















ay’, ++, y") 
O(x', +++, x") 


it is different from 0 at each point. 


#0; 
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A manifold is called orventable (two-sided) if it is possible to choose the 
jocal coordinates in the first place so that each such Jacobian (on an overlap 
of local coordinate neighborhoods) is positive. 

Example. We make the two-sphere S* into a manifold by using six 
coordinate neighborhoods. We set 


S? = {(z, y, 2) where a? + y? +27 = 1}. 
The neighborhoods are 

i ={x>0}, coordinate system y, z. 
Ur = {x < 0}, coordinate system z, y. 
Us ={y>0}, coordinate system z, x. 
Uz ={y <0}, coordinate system z, z. 
Us ={z>0}, coordinate system z, y. 
U; ={z<0}, coordinate system y, z. 


In comparing the overlap of two of these, we shall not be pedantic and intro- 
duce different letters, hoping the reader will forgive this sloppy notation. 
On the intersection of Uf and Uj we have the coordinate transformation 


ae 











Z= 2, zt>0, y>O 
and so 
Oy, 2) >= | 2 
Bz, 2) a, ao ils 
] 0 
On the intersection of Us and U;, 
aoe) 
_ —y? —z?, z>0, z<Q, 
aa ] 0 
»2 oa 
aya) 5. 5 rape 
On the intersection of UZ and U;, 
2=2 
—— y<0, z2<0, 
; 0 1 
aaa @ 
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Thus the two-sphere is a two-manifold, and our choice of local coordinate, 
proves it to be orientable. 

One could also cover the sphere S? with a system of only two local oo. 
ordinate neighborhoods by taking two opposite hemispheres, each extendeg 
slightly to make open overlapping neighborhoods. 





The sphere S? has two opposite orientations (outward or inward normal, 
corresponding to counterclockwise or clockwise sense of rotation). Similarly 
an orientable n-manifold has two opposite orientations. A definite one of 
these is determined by the order in which local coordinates 2, --- , 2" are 
given, up to an even permutation of this order. Making an odd permutation of 
local coordinates gives the opposite orientation. 

Let M be an n-manifold. To say that a real-valued function f on M is 
smooth at a point P of M means the following. Let U be a local coordinate 
neighborhood containing P with coordinates z',-+-,2z*. We require that 
f(x!,--+,2") be smooth near P. This restriction on f is independent of the 
particular U one chooses, since two coordinate systems whose neighborhoods 
overlap on a region including P are themselves related by smooth functions 
(from the definition of manifold). A real-valued function f is smooth on M 
if it is smooth at each point of M. 

Similarly, if M and N are manifolds of dimensions m and n, respectively, 
one defines a smooth mapping 


@: M—N 
by the requirements that in local coordinates z',---,2" on U in M and 
y',--+,y" on Vin N, we have ¢@ represented by smooth functions 


yi = y'(x', +>, 2") (v=1,-++, 2) 
on that part of U which ¢ maps into V. 
A manifold M is called a submanifold of a manifold N provided there is & 
one-to-one smooth mapping 
7; M—N 
which has this regularity property: in local coordinates (as written above); 
the matrix 
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has (maximal) rank m at each point. We refer to j itself as an injection or 
imbedding of M in N. 

This applies in particular when N = E” so that we may refer to submani- 
folds of Euclidean spaces. It is an established result of manifold theory 
that each m-dimensional manifold which is not too large may be imbedded 
in E" with 2 = 2m + 1. 


5.3. Tangent Vectors 


We study a manifold M and a point P on M. Our job is to define the 
tangent space at P, an n-dimensional vector space whose elements are the 
tangent vectors at P. Because we are not within the simple terrain of 
Euclidean space we cannot merely draw arrows emanating at P. We need 
a way of considering ordinary Euclidean vectors which depends in no way 
on arrows, or directed line segments. The answer is simple. We may 
identify Euclidean vectors with directional differentiations. Thus in case P 
is a point of E* and v = (a, b, c) is a vector at P, we may identify v with the 


operator 
(« d +6 : +c =) 
Ox oy Oz 


This does the usual things to sums and products, which motivates the 
following definition. 
First some notation. If M is a manifold, we denote by 


F°(M) 
the space of all smooth real-valued functions on M. 


Let P be a point on a manifold M. A tangent vector v at P is an operator 
v: F°(M)—> R, the reals 





P 


satisfying 

(i) vi(af + bg) = av(f) + bv(g), a,b constant. 

(i) v(f-9) = 9(P)-v(f) + f(P): vig). 

Thus v assigns to each smooth function f on M a real number v(f). 

We shall first observe that if we take a constant function c, then v(c) = 0. 
For setting f = g =0 in (i) yields v(0) = 0, settang tf = g = 1 in (ii) yields 
Vv(1) = 0, and setting f= 1, a =0, and g = 0 in (i) yields v(c) = 0. Observe 
that 

vicf) =cv(f) 


for any f and constant. 


Next, suppose z!,--:,2" is a local coordinate system, valid in some 
neighborhood of P. Then each of the operators 
7) 
Vi 


== 
Ox P 
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(the vertical bar means “evaluated at P’’) is a tangent vector, as one easily 
verifies. 

The totality of tangent vectors at P makes up a linear space, T,, calleq 
the tangent space to M at P. We shall show that these vectors v; form q 
basis of this tangent space. We set 

(xi, +++, a")| > = (cl, -+-,c”). 
If v is any tangent vector at P, set 
v(2') = v(x? — cc!) = a’, 
Now if f is any smooth function on M, we expand f in a Taylor series up to 
first-order terms with the integral form of remainder: 


f(x) = fle) + > (2 — ¢') 9,(x), 











of 
g(¢) = Bat 
Then _ 
vif) = vif(e)] + )) 9:(e) v(x’ — ¢') + >, (ce! — ¢') v(g;) 
. Of . of 
— 0 : —_— 0 = ‘ —_— . 
+e Sy] +0- 50S) 
hence 
7 
v= ; —_ 
ya 5 
which establishes the result. We refer to a’, --- , a” as the components of v 


with respect to the coordinate system x. If y is another coordinate system 
valid at P, and 


0 
oS Bh 
7 d dy’ P 
we find, by the chain rule, 





dy' 

b=) al — 
2 5a, 
the usual transformation law for contravariant components of a vector. 
Note here that we are working at a single point so that a and b are constant. 


A vector field on M consists of a smooth assignment of a tangent vector to 

each point of M. In local coordinates, 

0 

v=) a‘(x) Sgt? a‘(x) smooth. 

z 

On an overlap, 
ee 
v= b' =; 
» By) 3y) 


; oy’ 
b(()) = Dax) SS. 
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5.4. Differential Forms 

The smooth functions on M will also be called 0-forms. They form a 
space F°(M), the space of forms of degree 0 on M. 

We now define a one-form at a point P of M. We must have an 


expression 
a, dz", a, constant 
i i 


for each local coordinate system (z') valid in a neighborhood U which 
includes P and such that any two such expressions 


Yia,du', Yb, dy’ 
at P are related by 


dy’ 
ee ae 


the usual transformation law for covariant vectors. Evidently this is 
completely consistent with our local study in Chapter ITT. 

Having this, we may form sums of exterior products of one-forms at P to 
construct p-forms at P. Now we can define a p-form on M. This is a 
smooth assignment of a p-form to each point P of M. If U is given with 
local coordinates (z'), then on the neighborhood U we have the representation 


wo = Y ay(x) da? 


with smooth functions a,(x) on U, H = {h,,---,h,}. 
If we have the representation : 


o= > bx(y) dy* 


with respect to a second coordinate system which overlaps the first, then the 
relation between the b’s and a’s is given by substitution of y' = y'(x) for y! 
and 


dy’ 
3 aa dai 


for dy’. Asa consequence of our study of coordinate changes in Section 3.4, 
we see that the space 


F?(M) 
of p-forms on M is completely defined, that exterior multiplication 
OMAN 


of two-forms on M is accomplished by operating with one point at a time, and 
that the exterior derivative of a form on M is defined by working it out in 
each local coordinate system. 
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All the rules of Chapter III are readily verified, 

do An) =dw an + (—1)8% wa dy 

for example. 

If M and N are two manifolds and 

od: M—N 
is a smooth mapping, then there is a natural induced mapping ¢*, 
o*: FP(N)— FP(M) 


which again is defined by applying the local construction in one local co. 
ordinate system at a time and piecing together the results. As in the local] 
theory, we have the results 


(i) O*(@ +n) = ot w+ $*n. 
(ii) O*(A Ap) = (6*A) A (OB). 
(iii) d(¢* w) = $* (da). 


The last of these can be expressed by means of a commutative diagram 


o* 
F?(M) ~=&——— F?(\N) 


d d 


* 
FP+}(M)—e—— F?+1(N) 


Each of the two possible paths from F’(N) to F’*'(M) leads to the same 
result. 

In practice, one often constructs differential forms on a manifold this 
way. One knows in advance several smooth functions f,g,--- on M. 
From these one constructs one-forms df, dg, --- and from these in turn forms 
of higher degrees by taking exterior products. 


Example 1. On the two sphere S$? considered in Section 5.2, the functions 
x,y,z are smooth 0-forms. Thus dz, dy, dz are one-forms and dzxdy, dydz, 
etc., are two-forms. On the neighborhood U/ we have 


g*= 1 —y* —2?, 


a ee 
x 
ee re 
dx dy = bat is fa Bi dy = = dy dz, etc. 
x 
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Example 2. The circle S' = {a7 +y? =1}. Here x and y are functions on 


s! and 
adz+ydy=0. 
This means we can define a one-form a. At a point where x # 0, 
_% 
a 


At a point where x = 0, we have y # 0 and 


C=. 
y 
On any are of S' which is not the complete circle, we can find a function 6 
such that 
x = cos 8, y = sing, 
hence 


a = dé. 


It must be emphasized that no such function @ exists on all of S'—it would 
have to jump by 22 somewhere. 


5.5. Euclidean Simplices 


In this section we shall describe the standard building blocks which we 
later piece together to form fields of integration, p-dimensional spreads in a 
manifold over which we can integrate p-forms. These building blocks will 
be called Euclidéan simplices of various dimensions—we shall omit repetition 
of the adjective Euclidean in this section, but we understand that everything 
takes place in Euclidean space. 

A 0-simplex is a single point (Pp). 

A l-simplez is a directed closed segment on a straight line. It is com- 
pletely determined by its ordered pair of vertices (P,P,). 

A 2-simplex is a closed triangle with vertices taken in some definite order. 
It is completely determined by its ordered triple of vertices in the proper 
order, 


(Po, Py, P). 
Similarly one has a 3-simplex based on an ordered quadruple 
(Poy Pica) 


of four points, no three collinear. Geometrically it represents a tetrahedron. 
Finally, an n-simplex is the closed convex hull 


(Pog 084: Be) 
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of (n + 1) independent} points taken in a definite order. The geometrical 
set so spanned consists of all points 


P=to)Pot--'+t,P,, ¢:20, t,=1, 


i.e., all possible centroids of systems of nonnegative masses f,,-°°-,t 
located at P),°::, P,,, respectively. 

The boundary ds of a simplex s is a formal sum of simplices of one lower 
dimension with integer coefficients: 


OPo, Py.°°+s Pr) = Ye (-WPos Paves Prone Pron *** Pad 
An examination of the lower dimensional cases convinces one that this is 
consistent with the customary ideas on boundaries of oriented regions. 
OP, Py) = (P,) — (Po); 
OPo, Py, P2) =(P1, P2) — (Po, P2) + (Po, Pi), 
OPo, Py, P,, P3)=(P,, P2, P3) — (Po, P2, P3)+ (Po, Py, Ps) 


—(Po, Pi, Pa). 

Po P, 
o_O 
P, 

Po P 


In the triangle, the ordering of the vertices gives a sense of rotation of the 
triangle. In the tetrahedron, the ordering of the vertices gives a right- 
handed screw sense in space and induces a positive sense of rotation in each 
triangular face (outward drawn normal). One thinks of each minus sign 
in 0S as representing a reversal in this rotation sense. The result is that 
O(P,,-::, P) represents the oriented geometric boundary of the tetrahedron 
according to the outward drawn normal. 


J This means that the n vectors (P1— Po), (P2—Po),...,(Pna— Po) are linearly 
independent. 
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An n-chain is a formal sum 
where the a’ are constants and the s, are n-simplices. Its boundary is 
defined by 


dc = } a'‘(ds)). 





A basic result is that the boundary of each chain itself has zero boundary: 
o[dc] = 0. 
It stffices to check this for simplices. Let us try low-dimensional cases: 
OO(Po, Py, P2)) = O(P,, P2) — (Po, P2) + Po, Py) 
= {(P2) — (P;)] ~ {(P2) — (Po)] + (Pi) — (Po)] = 9, 
O[HPo, ++, Ps] =((P2, P3) — (Pi, P3) + (Pi, P)] 
~-[(P2, P3) — (Po, P3) + (Po, Pad] 
+[(P,,P3) — (Po, P3) + (Po, P1)] 
—[((Py, P2) — (Po, P2) + (Pos Pi) =9, 


which illustrates the general idea; each face occurs twice with opposite signs. 
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More generally, in computing 


OLA(Po.°**, Pad], 
one obtains 
(Poe Lists Faget a Fes pene 


twice, with opposite signs, once each from 


OP gs ok eis ass ee) 
and 
OP 8 o2 faie f jae sews 


so that everything cancels. 

Given two n-simplices (P,,---, P,),(Qo,°°*»@,), there is a unique linear 
correspondence between them which preserves the ordering of the vertices. 
It is given by 


t, = 1). 


i 


y t,P;<——> )tQ; (¢; 2 0, 
ry 0 


oM@.= 


It is convenient for defining integrals to have standard models of the 
simplices of each dimension. We define the standard n-simplex 


s"=(Ro,-+:, R,) 


as the simplex in E” based on 


R, = 
R, = (10---0) 


*®eeseeeereeere 





x 


We must now agree on a certain convention for integration. Let m be an 
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n-form defined on a domain U of E” which includes §". We wish to define 


We do this by writing w in the unique way 
w= A(a',---, 2")da! dx? --- dz" 


with the variables in their natural order, and then setting 


iz = {. A(x) dz! da? --+ dz", 


where the right-hand side is now the standard ordinary 2»-fold integration, 
which may be evaluated by any scheme of iteration, regardless of what order 
in which the variables are taken. 

For example, if wm = dzdydzx, then 


1 i=y 1-x-y 
| o=-| day a= — | ay | ax | dz = —1/6. 
s° > 0 0 0 


5.6. Chains and Boundaries 


Now we consider a manifold M and we shall define an n-simplex in M. 
As a preliminary definition, this consists of three things: a Euclidean n- 
simplex s", an x-dimensional neighborhood U of s" in Euclidean space,} and 
a smooth mapping @, 


od: U—M. 
We denote this preliminary simplex by 
(s", U, ¢). 
If we are given a second one, 
(t", VW), 


it will be considered the same as the first provided 


(5 uP = (> 101) (12 0, Y t= 1), 
1) 0 i) 
where 
s8=(Po,P1,°°*, Pa) t" = (Qo ,Q15°°° 5 Qn): 
In other words, if we set up the natural order-preserving linear equivalence 
+ That is, a neighborhood in the smallest flat submanifold of Euclidean space contain- 


ing s*. This smallest flat submanifold is the totality of points >3 &%P:, t real, > %=1, 
where s* =(Po,..., Px). 
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between s” and t”: 
s"<—> t", 


then $(P) = W(Q) whenever P and @ are corresponding points. This is also 
expressed by the commutative diagram 


Se 


® y 


The totality of these preliminary simplices (s", U, @) which in this way 
are identified with a single one make up an object which we call an n-semplex 
in M, denoted by a symbol o”. 

The open neighborhoods U we have introduced merely serve to eliminate 
difficulties with differentiability on the boundary. 


If o”" is a simplex represented by (s", U, @), then s” has faces ty, ---, t,, 
each a Euclidean (x — 1)-simplex, where 
ds" =) +t;. 


By restricting ¢ to the various t;, each extended a little in U to make open 
neighborhoods V,, we define the faces of o”, each represented by 


t= (t;, V', p) 
and the corresponding boundary 
do" = y +t; . 
U 
to 


This is an (n — 1)-chainin M. By an n-chain c of M we mean a formal sum 
c= > a; G;" 


with constant coefficients a, and »-simplices o7. Chains may be added and 
multiplied by constants. We denote by 


C,(M) 
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the set of all n-chainst on M. We set 


dc =) a,007 for c=) aj,oj. 
Thus 
0: C,(M) —> C,,_ ,(M) (n= 1,2,---). 


The basic property of the boundary operator 0 follows readily from the 
corresponding Euclidean situation: for each n-chain c, 


o(éc) = 0. 


A cycle is a chain z whose boundary vanishes, 0z = 0. 

A bounding cycle (or simply boundary) b is a chain which is the boundary 
of a chain of one higher dimension, b = 0c. 

Each boundary is a cycle, for if b = dc, then 


d(b) = 8(de) = 0. 


One further thing to be noted is this. In our preliminary definition of a 
simplex (s", U, @) we do not require that the smooth mapping ¢ on U 
into M be a one-to-one mapping. Indeed, it may happen that it takes all of 
s” into a lower dimensional space, even into a single point! A close analysis 
shows that not only is there no harm in allowing such “‘bad’’ mappings but 
that there are very great technical difficulties involved in attempting to avoid 
them. 


5.7. Integration of Forms 


Our data is a manifold M of any dimension, a p-form w on M and a p-chain 
con M. We must define 
| c. 
¢c 


c= ya, 0; 


where the a; are constants and the o; are p-simplices and write 


Jo-Zal 


so it remains to define the integral of w over a p-simplex 6. Now we can 
represent o in the form 


First we set 


(s’, U, ¢) 


where §? is the standard p-simplex in E? and ¢ is a smooth mapping of the 
neighborhood VU of Ss? into M. Our definition is 


t Precise topological terminology: ordered singular differentiable n-chains. 
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[o= [oro 


Since o* w is a p-form on U, this is an ordinary p-fold integral, as discussed 
in the next to last section. 

In application, one often does not bother to spell out in detail how a given 
geometrical region may be considered as a chain, but rather relies on the 
usual combination of experience and intuition, the latter an excellent guide 
in geometry. Forexample, suppose wisa 2-form on S* = {z? + y? + 27 = ]} 


and one seeks | w taken over $*. There will usually be a more effective 


procedure than using the coordinate planes to decompose the surface $? into 
eight spherical triangles, setting up mappings of the standard triangle onto 
each of these, etc. 

What then is the value of this rather long story on chains, boundaries, and 
integrals? In this age, it hardly seems necessary to defend the placing on a 
logical and rigorous basis things which are only understood in an intuitive 
sense. In addition, we have here a powerful theoretical tool as we shall see 
immediately in the following section on the general Stokes’ theorem. 

As an exercise, one could check that each of the standard tricks used to 
evaluate surface integrals, etc., fits into the above scheme of things. It 
hardly seems worth our time here. 


5.8. Stokes’ Theorem 


The general result we establish now includes all known formulas which 
transform an integral into one over a one-higher dimension spread. 
Let w be a p-form on a manifold M and ca(p+1)-chain. Then 


| w = | ae 
ec c 


Since ¢c is a sum of (p + 1)-simplices with constant coefficients, it suffices 


| | 
eo ¢ 


where o is a (p + 1)-simplex. According to a representation 
(s?*", U, 9) 
of 6 we have from the definition 


| a0 . o* (to) = | d(pb* w). 


This reduces the problem to a Euclidean one. Let 4 be a p-form on a 
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neighborhood U of s?*! in E?*+!. To prove 


font fot 
ogp+i eet. 


) = >, A(x) da? wer dai dyitt ... dyPt! 


so that it suffices to check the formula in case n is a monomial only. Since 


we may permute coordinates provided we are careful about signs, it suffices 
to take the case 


Now 


y= Adz’ ---dz?. 








Then 
0A 
dy = (—1)? BaP tt dx} es ~dgPt}, 
We remember that $’*' consists of all points (x!, --- , 2? +1) satisfying 
: pri 
“v“20 Yash. 
i 
We have 
OA 
| dn = (—1)? | Ageti dx} +++ daP*} 
gett gP +1 
(1- Lf x!) 
0A 
=(-1) | dz! + da? ( | sori de") 
{x'2 0, Lexis i} 0 
p 
=(-1)? JAets == 2%] — ¥'z') 
1 


{x'20, Laxist) 
— A(t, SRR. xP, 0) da eee dx? . 


We must next investigate ds’*'. We write 
gPtt = (Ro, R,, see Ry+1)) 


R, = (10---0) 
points in E?*?. 


es ¢ @ @ ee & & & 


Rya1 = (0++- 01) 
We have * 
ase*! = (R, 9°" y Ro+1) + (—1)?**(Ry, Rk, 97" y E,) 


+ other faces, 
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where 7 = 0 on each of the other faces since some one of 2’, - -- , z? is constant 
there. Thus 


| n= n+(—1)?*! N. 


égrt} (Ry,* °°, Rp+1) (Ro, Ry,***,Rp) 





The face (Ry, R,,--:, R,) is the standard s’. On it z?*! = 0 and so 


(—1)?*! | n= (—uptt | ae +++, a?, 0) dal ++ da? 


(Ro, ***; Rp) s? 


which is precisely the second term in the expression for | dy above. The 


first term is obtained by projecting downward in the x?*! direction: 


1 


P 
| n= A(t, -+-,2%,1— Ya!) dat «der 


(Ri, +++, Rp+1) (Ri1,°*+*) Rp, Ro) 
= (—1)? A(at,---,2%,1 = Drat)dat der 
(Ro, Ri, ++, Rp) 
= (pr [ale ee =‘) dz' +--+ dz?, 
i 
dp 


and this is the first term in the expression for | dn. The proof is completed. 


5.9. Periods and De Rham’s Theorems 


We consider an example. The manifold M consists of E> with the origin 
removed, 
M = E* — {0}. 
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Suppose w is a one-form on M such that dw =0. Then is w exact? That 
is, is it the differential of a function on M? The proof in Section 3.6 will 
not avail here because M cannot be shrunk toa point. Nonetheless, w = df, 


where 
x 
f(x) -( w, 
(1,0,0) 


the integral taken along any path ¢ which avoids0. That this is independent 
of the path follows from Stokes’ theorem. For if ¢’ is another path in M 
from (1, 0, 0) to x, then the chain ¢ — ¢c’ is the boundary of a piece of surface 
x (2-chain) in M and 


[o-| o=| o=| wm | domo. 
c e’ c-c’ Fb 3 £ 


Next suppose « is a two-form on M such that dx=0. We seek a one- 
form 1 on M such that a=d/. By the converse to Poincaré’s lemma in 
Section 3.6, such a form J exists locally. But we are asking the global 
question: Is there such a form A on all of M?. The answer to this one is no 
in general, we shall have explicit examples later. For if there were such a 
one-form A with dd = « we would have 


fake Le 


since the unit sphere S* has no boundary. But there is no reason a priori 
for assuming that 
| a= 0. 
$2 


The correct result is this. If «a is a two-form on M = E? — {0} with da = 0 


and 
[ens 
§2 


then « = dA for some one-form 1 on M. 

This result is contained in De Rham’s theorems which we shall formulate 
now without proofs. 

We deal with a fixed manifold M about which we assume only some mild 
limitation on its size, for example we may suppose it can be imbedded in a 
sufficiently high dimensional Euclidean space. 

A closed form is a differential form w on M satisfying dw = 0. 

An exact form is a differential form w on M satisfying w = dn for some form 
yon M. 

Each exact form is closed: 

dw = d(dn) = 0. 
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Let w be a closed p-form. To each p-cycle z on M corresponds a periodt 


is 


If z happens to be a boundary b = dc, the period vanishes, 


[o=| o=[ao-[o=o 
b éc c c 


Because of this there is a relation between periods: 


of Qo, 


Whenever cycles z, ,--- are related by 
>. a;z; = boundary, 
(t){ then 


Ea| ono 


Dr Ruxam’s First THeorReM. A closed form is exact if and only if all of 
its periods vanish. 


Dr Ruam’s SECOND THEOREM. Suppose to each p-cycle z is assigned a 
number, per(z), subject to the consistency relations 


whenever 
> a,z; = boundary, 
(f) 
then 
>. a; per(z;) = 0. 


Then there is a closed form won M which has the assigned periods, 
| @ = per(z) for each p-cycle z. 
Zz 


On many spaces one is able to apply these results because there is a finite 
set of independent p-cycles which spans all p-cycles, up to boundaries. For 
example, on the n-sphere $" it is known that each p-cycle is a boundary for 
p> 0, p # n, and that in dimension x there is a single n-cycle (S" itself with 
outward normal for orientation) such that each n-cycle is a multiple of this 
one plus a boundary. These things are established by algebraic topology. 

A complete analysis of De Rham’s theorems reveals the following result, 
which has considerable attraction in itself. 

Suppose we consider only chains c = ) a,6, which are sums of simplices 


. + The nomenclature derives from the periods of elliptic integrals and the corresponding 
differentials for algebraic functions. 
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with integer coefficients. Then we may talk of these as integer-chains and 
have integer-cycles and integer-boundaries. The integer-periods 


Ie 


of a closed form w are the periods taken over integer-cycles only. 
Let m and n be closed forms of degrees p and q respectively. Suppose that 
the unteger-periods of w and n are all integers. Then the same is true of w A n. 


5.10. Surfaces; Some Examples 


It is shown in topology that each closed surface in E*? may be smoothly 
deformed into a sphere with h handles, or alternatively, a button with h 
holes. Let us consider the case A = 2 and orient this surface with the out- 
ward drawn normal. The only significant two-cycle is L itself. By De Rham’s 
First Theorem, a two-form « on this surface is an exact differential if and 
only if 





normal 


There are four significant one-cycles, ¢,, ¢}, €,,¢,. Here ¢, and ¢, intersect 
once and cross, the same for ¢, and c,. But c, and ¢, intersect once 
without crossing. To see the geometric plausibility of the statement that 
each one-cycle ¢ on £ is a sum of multiples of the ¢, and ¢; plus a boundary, 
one cuts the surface LZ along these basic cycles. Having done this, & may be 
Smoothly deformed into a plane domain without holes. 

De Rham’s First Theorem now asserts that if q@ is a closed one-form on Z, 
then @ is an exact differential if and only if 
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| o=| o=| o=| o=0. 
cy Cc’; €2 C’2 


Applied to dimension one, De Rham’s Second Theorem asserts that if real 
numbers @, , 4; , 2; , @, are given, there exists a closed one-form w satisfying 


| @=4a,, | O=4a,, [oma | @ =a). 
Ci C1 C2 C’> 


Cc; C2 


Cc; & C, 


It is also interesting to consider non-orientable closed surfaces. These of 
course cannot be realized in E?. Perhaps the simplest is the projective 
plane P?. This is defined by pasting the edges of a rectangle together in the 
order indicated. The boundary relations are 

6(P?) = 2¢’ — 2e, 
de = (P) — (Q) 
dc’ = (P) — (Q). 
This means first of all that there is no effective two-cycle, each two-form is 
exact. The only effective one-cycle is c’ —c, and this actually bounds, 
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c’— c= 3428P*. Thus each closed one-form is exact. 





Another interesting example is the Klein bottle K*, again defined by 
pasting edges together. The boundary relations are 


OK? = —2e¢ 





The one independent one-cycle is ¢’. 


5.11. Mappings of Chains 
Suppose M and N are manifolds and f is a smooth mapping: 
f: M—N. 


Then to each p-chain ¢ on M there corresponds in a natural way a p-chain 
f,con N. 

It suffices to explain this when c is a simplex o”. Such a simplex is 
represented by (s?, U,@) where U is a neighborhood of the Euclidean 
simplex s? and @¢: U-—>M. We merely compose f and @¢ so that f, ¢ is 
represented by 


(s?, Ufo p). 


Uy 


? 
M 
f 
fog | 
N 
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We illustrate the process for the case of a two-simplex 





This induced map f, takes the space of chains onto the space of chains: 
M—N 
C,(M) = C,(N). 
We observe that if ¢ is a p-chain in M, then 
fx (00) = O(f, ¢), 


which leads to the commutative diagram 


7) é 


C,-.(M) > G -(N) 
which is certainly analogous to the corresponding diagram in Section 5.4 
for f* and d. The validity of the result is established by looking at in- 
dividual simplices. 

We now see what happens with two mappings. Let 


f 
al. | 
g 
gof 
P 


Then the assertion is 
(Go f)x =Ge otn 
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which again follows for a simplex almost directly from the definition of f, . 
Finally we consider this situation. Let 


f: M—N. 


Suppose that w is a p-form on N and Cc is a p-chainon M. Then f*q@ isa 
p-form on M and f, ¢ is a p-chainon N. We have 


[fr -| Q. 
c fee 


This important result also follows directly from the definition for a simplex 
and is obtained for a general chain by summation. 


5.12. Problems 


1. Show that the totality of unit tangent vectors to the sphere S$? is a 
three-manifold. Construct local coordinates. 

2. Show that the set of all directed lines in E? is a 2-manifold. Discuss 
orientation. 

3. More generally, consider the set of all oriented r-dimensional planes in 
E”. Show that this is a manifold, compute its dimension, and discuss 
orientation. 

4. Projective n-space P" consists of all (m + 1)-tuples (ag,--+,a,) of 
real numbers not all zero, where proportional (n + 1)-tuples are considered as 
representing the same point. Show that P” is an »-manifold. 

5. Complex projective n-space CP” consists of all (n + 1)-tuples 
(@9,°**,@,) of complex numbers not all zero, where two such n-tuples are 
considered the same if they differ by a (complex) proportionality factor. 
Show that CP” is a manifold and determine its dimension. 

6. Show that the manifolds of Examples 4 and 5 are closed (compact). 

7. Let M be the manifold of Example 3. Show that the set N of all 
oriented r-planes in E” which pass through the origin is a closed submanifold 
of M. 

8. Denote by U the open region 

U ={xit+---+22>1} 
in E”. Suppose w is an r-form in E” which vanishes identically on U. 
Under what conditions does there exist an (r — 1)-form « on E” which also 
vanishes identically on U and which satisfies dx = @? 

9. Show By direct calculation (i.e., without De Rham’s theorems) that 


if @ is a two-form on S? whose integral over S’ vanishes, then w = da for a 
suitable one-fornt'x on $7. 


Vi 


Applications in Euclidean Space 


6.1. Volumes in E" 


We denote by 
@o=dz,:--dz 


the element of volume in E”, an n-form, and set 


2 2 
v= | Q, r=) fi 
rsil 


so that V,, is the volume of the unit ball. Next we denote by o’ the element 
of (n — 1)-dimensional volume on the unit sphere $"~' = {x |r = 1}, and set 


A,-4 -| Oo’. 


Thus A, = 2x, A, =4n, V,; =2, V,=n, V3;=4n. It is clear that the 
volume of the sphere of radius r is r"~'A,_, , hence 


1 
1 
V, -| Pl dr=—A,_1. 
0 nr 


One may evaluate V, by integrating over slabs: 


1 
V, | (l-22)@-D2y dx 
=—{ 


= Vp Yn 
where 


1 
J, -| (1 — x2)" D2 dy. 
ase | 


Integration by parts once leads to 





1 
J, -| x(22)(” 3 -)a — 27)" 9)? day = (m ~ 1)(—F, + Iya); 
al | ant 
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AF\ 
ares 


These recursion formulae lead to the standard result 


n{2 
Te 
V.= 


" (44) 

(5+) 
Next we obtain an explicit formula for o’ in terms of the Euclidean co- 
ordinates z,,°°:,2z,. We begin with the form 


rdr =) x,dz,, 


a one-form in E” which is invariant under rotations (orthogonal transforma- 
tions) of E”. Consequently 


YN 
ardr = ) (—1)'"'2,dx, ---dx,-++dax 


n 


(the “hat”? denotes a missing factor) is an (nm — 1)-form in E" which is in- 
variant under rotations. It follows that on $"~/, 


o’ =cx«rdr, 
where.c is a constant. 
Next we note that 


: YN 
d(*rdr) = ) (—1)'"'dz,dx, ---dx,-+- dx, = nw, 


4,1= | o=e| vrdr=e| d(« rdr) 
n-1 gr~i rsl 


4 =| na =cnV, =cA,_1, 
rs1 


hence 


c=l,o’ =«rdronSs"'. 
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Summarizing, we set 
n “N\N 
o=*rdr=) (—1)' 'x,dx,---da,---dz,, 
i= 
defining an (x — 1)-form o in E”. Then do = nw, and if o is restricted to 


S"—! the result is the (n — 1)-dimensional volume form o’ of $"~'. 
Next we consider the natura] projection 


nm: E"—{0}—+$""! 


defined by (x) = x/|x|. 
We seek x* o’, an (n — 1)-form on E" — {0} satisfying 


d(nx* o’) = 0 since d(x* o’) = n* (do’) = n* (0) = 0. 
(do’ is an n-form on $"~', hence 0.) We shall prove 
ee eed 
wg =. 
We could prove this by directly substituting 


_ A 
y, = 2,/7 in c= ¥ (-1)! ly,dy,:+:dy;°+- dy,, 


but we prefer to proceed indirectly by exploiting the symmetries present. 
We set 


Then 
l 
dt =< do — 5 (rdro = - 5 = 0. 


Now we observe that *(z*o’) is a one-form in E" — {0} which is invariant 
under rotations, hence dependent on r alone. We may write 


«(7* a’) - (rdr). 
From this we have 
n¥o! -— o = fir)t, 
edinto\an:. We : 
Veen eS t a J =, 


a constant, z*o’=ct. To evaluate c, we simply note that on $"~!, both 
n* co’ and t collapse to o, hence c = I, 


n* o' =T. 
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6.2. Winding Numbers, Degree of a Mapping 
A basic result of topology (Seifert und Threlfall [20], p. 283) asserts that 


if M and N are closed oriented n-manifolds and f: M— N, then the 
chain f, M is an integral multiple of N plus a boundary. This integer 
multiplier is called the degree of f and written deg f. 

Now suppose that Z is a closed oriented (m — 1)-manifold in E" — {0}. 
Then by the Jordan—Brouwer theorem of topology, £ decomposes E” into 
exactly two regions. We assume & is oriented by the outward normal. The 
projection mapping z of Section 6.1 sends £ into $"~!. It is true that 
deg x = 0 or I; our point is that this can be determined by an integral. Let 


d = deg z. 
[== [=o =| = | o = 6A,_,, 
x x n(E) s"-1 


Then 


hence 





More generally, let M"~’ be a closed oriented manifold, 
f: M"-'—+E" — {0}. 


Essentially we are thinking of f(M"~*) as a hypersurface in E” — {0} which 
may intersect itself. We look on this hypersurface as winding around the 
origin and we want to count how many times it encircles. This winding 
number is given by the Kronecker integral 


/ * 
|, - 


We may justify this as follows. Setg=2.f: M"~'—+S""'. What we 
are after is degg. Now 





gx (M) = (deg g)S"~* + (boundary), 


mM ge M 


= eg) | o = A,_, degg, 
sn 


1 
g* a’. 
A,- 1 2 


hence 





degg = 
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Me —_2—> * ~ (0) 
TT 
= 7 of 
g"-1 


But g* o’ = (f* o n*)o' = f*t, so we have 


l 
degg = | Die 
Ay 2s M 


The most general situation is this: 





f: M’— N’. 


Let B be the volume form on N taken so that | B=1. Then 
N 


deg f= | f* B. 
M 


For f, M = (deg f)N + (boundary), hence 


| rp=[ p = (aes f)| B = deg f. 
M SeM N 


One interesting example: let T” be the n-torus, f: S"—> T" where n = 2. 
Then deg f = 0. 

Because the integrals involved are integer-valued, they remain constant 
when the mapping in question is subject to a deformation. Precisely, let 
f,: M—+N bea one-parameter family of maps. Then 


dea f= | f,*B 
M 


is a smooth function of t, always an integer, hence constant. It follows that 


deg fy = deg fy. 
One other remark. Suppose we have 


f: M—N, g: N—P so that h=g.of: M— P. 


Then 
deg h = (deg f) - (deg g). 
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6.3. The Hopf Invariant 


For each sphere $", let a, denote the element of area, normalized so that 


Consider first a map f: S*—>S$?. Then f*a, is a 2-form on $*. Also 
d(f*o,) =f*(do,) =0. Since S* has no nontrivial 2-dimensional cycles, 
we deduce that 


f*o, = da, 


where the one-form a, on §$° is unique up to the differential of a function. 
The 3-form «, A f*a, has an integral 


(ct A f*o, ’ 


which has the remarkable property of being an integer, called the Hopf 
invariant of f. It is invariant under deformation of f. More generally, let 


f: g2"-1 —>» §". 


Then f*o, = da,_,, and the Hopf invariant of f is 


\, n= Af*o,.- 


We may represent S° by pairs of complex numbers 
(z,w), fz]? + foo]? = 1. 


The mapping (z, w) —— z/w provides a mapping of S$? into the closed complex 
plane, i.e., the Riemann sphere $?._ This map has Hopf invariant +1, hence 
it is essential in the sense that it cannot be deformed to a trivial map, every- 
thing going to a single point. 


6.4. Linking Numbers, The Gauss Integral, Ampére’s Law 


Let M’, N* be oriented closed manifolds in E”, where r + s = n — 1, and 
suppose these have no common point. (Best example: two disjoint closed 
curves in E?.) We want to count how many times they link. To do this, 
we form the product space M x N which is an oriented manifold of dimen- 
sionr+s=n—1. Weconsider the mapf: M x N — E" — {0} defined 
by 

f (x,y) =y— x. 
Now we set 
link (M, N) = deg f. 
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Thus if 7 is the n-form in E” — {0} we considered above, 








link(M, N) = 
MXN 


We shall work this out in E? for a pair of closed curves M, N: 
l dz):d 
iz = EF 2.dz jdt = iz eee. 


We let x, y be the moving points on M, N, respectively. Then 


t= 


z=f(x,y)=y— 
so that 
pe ge ee 
1 
= oy — x {[—(y — x) x dy]-dx —[(y — x) x dx]-dy} 
1 
ss ly os LY — x) x dy]-dx, 
link (M, n=s | «(2 Aaa SL 
ly —x| 


(In this computation dy x dy = 0, etc., since dy involves only one variable.) 

Imagine a steady unit electric current flowing around the closed loop N. 

By Ampere’s law, the magnetic field at a point x due to the current in a 
segment dy is 

1 (y — x) x dy 

4n jy—x|> ” 


hence the total magnetic field at x is 


_ ol (y — x) x dy 
r= aa) T= xP 


N 
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It follows that link (M, N) = | F(x)-dx is precisely the work done by this 
M 
field on a unit magnetic pole which makes one circuit of M. 


In the next example, link (M, N) = 0, which seems surprising since the 
curves cannot really be separated. 





Vil 


Applications to 
Differential Equations 


7.1. Potential Theory 


We summarize the notation of Section 6.1. Space: E". Coordinates: 
Uy, %y, °°", L,. 


rar = 2 wide; 
“nN 
= «(rdr) = Ya «dx; = y(- 1)'"!x,da,--+dx,++-da,. 
on o 
== 
w = dz, ++: dz, = volume element of E”. 


do = nw, dt = 0. 


nl? 
V.=| we, A_,=] onl. 
I. P((n/2) + 1) ; z 


Let u be a smooth function on a domain in E", Then 


du = Tae 


du YN 
edu = ) (—1)"! dat ves dates da", 


a? 
dadu = (y san) = (Au)o, 


defining the Laplacian 
au 
Au = »y ax® ‘ 
(See Section 4.4 for details when n = 3.) 
82 
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If u and v are functions on the finite domain R, the Dirichlet (bilinear) 
integral is 


du \ / dv 
D ; = du d = a d = its J y 
Iw, 2] {, ee {, ee {> (=) (=) 


Next we have by Stokes’ theorem 


| uxdv -| d(uxdv). 
éR R 


d(uxdv) = du a «dv + ud «dv 


But 


= du a «dv + uAvo, 
hence we have 


GREEN’S FORMULA 


| uxdv = Dfu, v] + | wave 
a 


By reversing u and v and subtracting the results, we obtain 


GREEN’S SYMMETRICAL FORMULA 


| (uadyv —vadu) = | (u Av —vAu)o. 
éR R 
[One usually writes «du = (0u/dv)A where / is the (n — 1)-dimensional volume 


element on OR and du/dy is the normal derivative.] 
In case v is harmonic in the region R, Av = 0, and we have 


| (us+dv — vxdu) +{ vAuw = 0. 
oR R 


By specializing further we have this result: 
Let u and v be harmonic functions in a region R. Then 


| u *dv -| v «du, 
éR éR 


We derive further consequences by setting 


1 

0 2? 
—(n—2 

dy = "=" (ran), 


yr" 
adv = —(n — 2)t, 


dadv = 0, Av = 0. 
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The function v is defined on E” — {0}. We suppose the region R contains 9 
and we apply the formula above to the punctured region 


R —{r<e} 


with ¢ a small positive constant. We suppose uw is harmonic on all of R, 


Since 
a(R — {r < c}] = OR — {r =}, 


| usar — | wsdv = | vedu — | v du, 
éR r=e oR r=e 


] 1 
—n—2)| a + m2) | as -{ ro) du — {_ co 


We evaluate the individual terms: 


I I 1 
ut = — us = — d(uc) = — (duAng+ung). 
r=e e r=e@ é rse E rse 


Now for |x| < ¢, u(x) = u(0) + Ofe), hence 


we have 


xAU. 








| um = u(O)V, + O(e)e". 
rse 


Similarly 


| duno -| ( sit) = | O(e)w = O(e)e”, 
rSe rse Ox rse 
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hence 


| ut = nV_u(0) + O(e) = A,,_ ,u(0) + O(c). 


: d : d : d(«d 
“0% = *aUu = * 
r=e oo en? r=6 gs rse - 


1 
mm =| (Au)w = 0. 


We substitute these results and let e —> 0 to get 


0 1 . 1 «LU 
— U SS ee 
mu A,—1 éR : (n ws 2)A,~4 aR ae 


which gives the value of a harmonic function at a point in terms of the 
boundary values of it and its normal derivative. 

Special case. Let R be the spherical region of radius a centered at 0, 
R= {r < a}. In this case the second term on the right-hand side vanishes 
for the same reason that the corresponding integral taken over {r =} 
vanishes. On OR = {r =a} we have 

















1 
qrat Le 





Co 
= 
where 


] . £% 
H= My = —) (Veda, +++ daz+ ++ dx 
a 


is the element of (m — 1)-dimensional volume on {r= a}. (For a= 1 this 
reduces to ¢. Since there are  x-terms in the numerator and a = |x| is in 
the denominator, it is homogeneous of the right degree, — 1.) We have the 


Gauss Mgan VALUE THEOREM 





1 oo 
es — vrta ; 
u(0) a"-1A._, {0 | P 


which tells us that the mean value of u over the sphere of radius a is the value 
of u at the center. 
Two important properties of harmonic functions follow from this result. 


Maximum (Minimum) Princreie. Let u be harmonic on the finite region R. 
hen u never assumes its maximum (minimum) value at an interior point of 
unless u is constant. 
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For suppose u assumes its maximum at an interior point x, of R which, after 
translation of coordinates, may be taken to be 0. Let £ be any (n — 1)- 
sphere of radius a centered at 0 with a so small that {r <a} isin R. Since 
u(x) < u(0) we have 


u(0) = | om /| s | wom /| w= m0 
= zr >» x 


so we must have u(x) = u(0) for allxin&. Hence w is constant on the largest 
spherical neighborhood of x) we can draw in R. Evidently this means wu is 
constant in al] of R since we can reach any other interior point by a sequence 
of such overlapping spheres. The result for minima follows the same way. 


UNIQUENESS PRINCIPLE FOR THE BOUNDARY VALUE PRoBLEM. Let u and 
v be harmonic on a finite domain R and coincide on OR. Thenu=vonR. 


For u — v vanishes on OR and is harmonic. By the Maximum Principle, 
u—vsO0onR,usv. Similarly v < wu, hence u = »v. 

The function (1/r’~2) is ideally suited to the sphere. On other domains 
it is inconvenient because of the term involving the normal derivative in the 
expression above for u(0). Hence we introduce the Green’s function. 

Let R be a finite domain. A function v(x, y) defined for x and y distinct 
points of R is called the Green’s function of R if 

(i) For each fixed y in R, v(x, y) is a harmonic function of x for x in 
R — {y}. 
(ii) For each fixed y in R, v(x, y) = 0 for x in OR. 
(iii) For each fixed y in R, 


] 
Ix — yl"? 
is a smooth harmonic function on all of R. 
Using the same method as above one proves the following: 


If u is any harmonic function on a finite domain R and v = v(x, y) is thé 
Green’s function for R, then 


u(x, Y) ar 


—] 
uly) = | *dy0(X, y). 


In case R is the spherical domain r < a centered at 0, the Green’s funa- 
tion is 


(x, y) : = 
u(x, ¥) = ~—— 3 — 7 an 
saa ly|""?}x — wi y 
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We note that 
2 
y 
ly?? 
is the inverse of y with respect to the sphere R (reciprocal radii). For x on 
OR, |x| = a and we have 


f 


Y= 


3 ; 2 a* 
|x — y'|* = (x — y’)-(x — y’) = a? —-2—, + 
ye ye 
oly? —2x-y + Ix) ay 
= — = —_—_- 1X — F 
lyl? ly|? : 
a 
Ix—y] = — |x— yl. 
mY 
This explains why v vanishes for x on OR. 
Next 
qa" 2 
sayy = ~(n- 9] ex) — TS x -y)] 
where 
l 
t(x — y) = ——_, ) (2; — y;) #dz;,, 
Y= ay hm 
(x= y!) = (a, — yd 
U(x — = ——_—— Li — Yj) *Az;,. 
y |x — y'|" i y 


We only need dv for x on OR. Recalling that y’ = (a?/ly|)y for |x| = a we 
have 


= (n — 2) lyl? a” 
*dyv ied = 1x— y|" —yh >: ce =U) = Aa r (« = <M) "dX; 
_ —m—2)a*—ly?? _ — —(n — 2) a? — ly)? 
 ixayl a? Ix—y oa °* 


so that the representation of u in terms of its boundary values specializes to 


2 12 
pe esd di { un) 


m 
GA,-1 Jixj=alX— yl" * 


This is the Poisson InrEGRAL FormvuLA which provides an explicit solution 
formula for the Dirichlet problem on the sphere—find a harmonic function 
with prescribed boundary values. 

Returning to a general finite domain, we mention the important symmetry 
property of the Green’s function: 


o(x, y) = oly, X). 
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(That this is the case for the sphere is not apparent from the unsymmetrical 
formula above. But for the denominator of the second term we have 


a? a* 
lyl2}x — y’? = lyl2Q¢— y’)-(Q¢— y’) = iy? (bx? ype OEM) + wa) 


= |x|7|y|? — 2a7(x-y) + a*, 


which turns out to be symmetrical after all.) 
One of the many important consequences of the Poisson Integral Formula 
is the 


LiouviLLE THEOREM. Let u be a harmonic function on all of E" and u = 0. 
Then u is constant. 


We shall show that for each y in E”, u(y) = u(0). We fix y with |y| = 5 
and select anya>06. Then 








i a* — “| u(x) 
uly) = ra 
4 GA,-1 Jjxj=xal* — yl" ; 


Now 
|x — y| S [x] + ly| =a +4, 
|x — y| = |x| — ly| =a — 4, 
hence 
1 1 | 
< oe 








Ss ; 
(a +5)" |x—y|’ ~~ (a ~ 6)" 
Since u(x) 2 0 we may use these inequalities to estimate the integral: 
(a* — b’) 


[ett Suy) S aw |. Uy - 


But from the Mean Value Theorem, 


(a? — b?) 
aA, _ (a 2 b)" 


| U(X) Hy = (a"~"A,_ 1) u(0), 
{x|=a 


so we have 
2 _ p? n~2 2 __ n-2 
—— u(0) $ uly) < — u(0). 


By letting a ~~» co we have 
u(0) S u(y) S u(0), 


and so for each y, u(y) = u(0), uw is constant. 
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Remark 1. We return to the symmetrical Green’s formula 


| (uxdv — v«du) -| (u Av — vAu)ow. 
aR R 


We apply this in this situation: 
R={e<srsa}, 


v harmonic in R for all ¢ > 0, 
v vanishes on 7 = @, 
u smooth in {r < a}. 


We do not require that u be harmonic. The formula reduces to 


| usar — | waar + | vadu + | v(Au) wm = 0. 
r=a r=e r=e esrsa 


First case. 





By the methods above, 


f —(n — 2) 
| a ua = —(n — 2)A,_,u(0) + O(e), 


| v du 


é r=e 
] ] 
(= = =] [aoa = O(e?). 
Substituting these in and letting e —+ 0 we obtain 
0 ] I ( ] ] % 
= ——— ti — - -—- 
wy a" "A, r=a - (x —2)A,_, rSa ied =) ae 


This gives information about a solution u of the Poisson equation Au = f 
with boundary values of u assigned. 


Second case. 


v= 


, 


418 
S18 


] l 2x, 
xdv = ss ade, — 0 S73 9. 


One differentiates this to prove Av =0. This also follows when one notes 


that 
0{1 —(n — 2) : 
ga (3-2) & 
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The end result in this case is 


Ou n 1 (- =)(a 
SF : — t{-- u . 
02 jlo aA, Paes Up Biss i t yh a" w 


Analogous formulas for higher derivatives are possible. 








Remark 2. We have avoided » = 2. In this case the basic difference is 
that the symmetrical harmonic function with singularity at 0 is Inr rather 
than r-“~”), Using this, results similar to those above follow. 


7.2. The Heat Equation 
We consider the parabolic equation 


O7u ss ?u—ss Bu 

Fx + ay? “3° 
Suppose u is a solution, valid in a region of x, y, space which includes a 
region R and its boundary. 

First we consider 
a = (u,dy — u,dx) dt — udxdy. 
Then — 
da = (u,, + u,,)dxdydt — u,dtdxdy = 0, 


| a= | da = 0. 
oR R 
Next we consider 


B = 2u(u,dy — u,dx) dt — u? dxdy. 


hence 


Then 

dB = 2(u,dxz + u,dy)(u,dy — u,dx) dt + 2u(u,, + u,,)dxdydt — Quu,dtdxdy 
= 2(uz + us) dx dy dt. 

It follows that 


| [2u(u, dy — u, da) dt — u? dxdy] = 2| (uz + u?) dx dy dt. 
aR R 


Suppose R is taken in the special form of a cylinder T > [0, 6] where T is a 
region in the z, y-plane. We then have 


OR = (6T) KX [0,6] +T K {0} -T X {0}. 


We now assert the basic uniqueness theorem: 
If u vanishes on the base TX {0} and on the lateral surface OT XK [0, b], then u 
vanishes identically in R. For the integral formula above reduces to 


2 (u2 + uz) dxdydt + | u(x, y, 6)? dudy = 0. 
R T 
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Since everything is positive, this implies 

U,=U,=O0 in 

0: on, va {6}, 


which is more than enough to imply u=0in R. Because the heat equation 
is linear, we deduce that two temperature distributions which coincide 
initially at £ = 0 and always coincide on the boundary of R must be identical 
for all ¢ at each point of R. 





We shall now do the same thing in n dimensions, where there is an inter- 
esting sign change. Our variables are z,,---, z,, and the heat equation is 


Au = du/dt 


where as usual Au = )' 0?u/d2?. The operator * will apply to space 
variables only. 
This time we set 
B = 2u(*du) dt + (—1)"" !u*o, 


where w = dz, ‘+: dx Now 


du a (xdu) = (grad u)?@ = » (=) 20 
i 
and we have 


dB = 2(grad u)?@ dt + 2u(Au)wdt + 2(—1)""!uu,dtw 
= 2(grad u)*w dt, 


| p= 2 (yrad u)*@ dé. 
éR R 


hence 
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Now let T be a region in E*, R=T X [0,6]. Then 
oR = oT X (0, b] + (-1)"T K 0, d) 
= oT X [0,6] + (-1)"T K b+ (-1)""'TX 0. 
Suppose w vanishes on 6T X [0, b]andonT X 0. Since ds=0 on TX 8 
(i.e., £ = 6 = constant) we have 


uxdudt = 0, 
aR 
consequently 
(—1)"7! | u?@ = 2 | (grad u)*w dé, 
éR R 
that is . 


| uo + 2| (grad u)*@ dt = 0 
TXb R 


and we conclude as before (gradu)? = } (du/Oz,;)? = 0 on R, du/dz; = 0 on 
R, u is constant on R, u = 0. 


7.3. The Frobenius Integration Theoremt 


Everything is local in this section; we operate in a neighborhood of 0 in 
E". Let @ be a one-form which does not vanish at 0. We ask, under what 
conditions are there functions f and g satisfying @ = fdg? In other words, 
we seek an integrating factor for the differential equationw =0. Ifw=/fdg, 
then f does not vanish in a neighborhood of 0, hence 


dw = df ndg=df a f~'a, 
dwo=O@An0 (0=f 'df=din|f)|) 


and so 
oAdwoa=oAGAw=0. 


For a one-form w = Pdzx + Qdy + Rdz in E?, this is the condition 
P(R, — Q,) + OP, — B,) + R(Q, — P,) =0. 


We note that if w = fdg, then the equations w = 0 and dg = 0 are the same 
and hence the solutions or integral surfaces of wm = 0 are the hypersurfaces 


g = constant. 
Before passing on to precise statements we give two instructive examples. 


Example 1. Let w= yzdzx+axzdy + dz so that dw = ydzdx+xdzdy. It 


follows that 
d 
dw = (=) A@ 
Z 


+ Material in this and the next section is taken from a University of California Tech- 
nical Report of June 1957, Seminar on exterior differential forms. This was prepared 
for the U.S. Army Office of Ordnance Research Contract DA-04-200-ORD-4356. 
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which is not so useful since dz/z is singular along the z-axis. A better choice 
is 9 = —ydx — xdy and we havedw=6@”Awq. To determine the function g, 
we use the fact that each integral surface g = constant will be cut by the 
plane {x = at, y = bt} in a curve which intersects the z-axis in the solution 
z of g(0, 0, z) = constant. The equation @ =0 on the plane x = at, y = bt 
becomes 

dz + 2abztdt = 0 


with solution 
z = cexp (—abt’) 


satisfying the initial condition z(0) =c. However abt? = xy so these curves 
span out a surface 
z=ce *, 


We now think of a, b, c as variables and make the transformation 


“x=a 
Oz, y,2) 
=b th ————_ =e” : 
- a O(a, 5, c) ee 
z= ce 2 
We have 
dz =e” "de — z(adb + bda), 
which yields 


w=e de, 
or in the original variables 
w=e ~ d(ze*’) 
and the integral surfaces are 
ze*” = constant. 


It will be observed that we have arranged the function g so that g=c 
intersects the z-axis precisely in z = c. 


Example 2. This time we try the procedure on w = dz—ydx-—dy. On 
the plane x=at, y= bt, the equation w=0 becomes dz = (abt + 6)dt, 
z= tabi? + bt +c and we arrive at the surface 


z=frytyte. 
But on the parabolic cylinders x = at, y = bt? we have 
dz = (abt? + 2bt)dt, z= abt? + bi? +6, 
z=doytyte, 


a different family of surfaces. The reason for this failure to obtain integral 
surfaces is seen from 


dw = —dydz, wo vA dw = —dzdydzx # 0. 
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TurorEM. Let w = ) fdx' be a one-form which does not vanish at 0. 
Suppose there is a one-form 0 satisfying dw =@ Aw. Then there are functions 
fand g ina sufficiently small neighborhood of 0 which satisfy w = f dg. 


Note 1. Since @ is given while @ is certainly not uniquely determined, it 
simplifies the proof if we avoid explicit use of @ as long as possible. 


Note 2. The condition on w is unchanged when we replace w by a multiple 
of w. In fact, if h 4 0, then 


d(jhw)=dhnaoat+hdwa=dhawth6nw= (dh+h0) ao, 
d(heo) = [(dh)h~* + O](ho). 


Proof. Since w # 0 at 0, we may assume some one of the functions f; does 
not vanish at 0. Since neither the hypothesis nor conclusion changes when 
we multiply w by a nonvanishing factor, we may assume 


o=dz—) A,dz', <A,=A,(x, 2). 
I 


a 





a = (a’,.-., a") 


We fix any point a in x-space and consider the equation w = 0 on the 
hyperplane z! = at, (¢ = 1,---, n): 


. =) A,(at, za’. 
We solve this equation with the initial condition 2(0) =c¢. More precisely, 
we seek a function F(t, a, c) satisfying 
F(t, a,c) = ¥ Ajat, F(t, a, c)Ja' 
F(0, a,c) = ¢. 


7.3. THE FROBENIUS INTEGRATION THEOREM 95 


The usual existence theorem of ordinary differential equations yields a unique 
solution. We see that a change of scale is possible: 


I 
Fit, a,c) = (u, k a, | 


since the function on the right is again a solution to the same problem. In 
particular, setting k = 1/t, 
F(t, a,c) = F(1, ta, c). 


We introduce the change of variables 


x=U 
z= F(1, u, v) 




















with 
O(x,z)} | f 0 “9 
Au, v)l, j* ll, 
since 
C d d 
— Fil, u, = — F(1, 0, = — F(0, a, 
a we) =F POLO =F FOa0)) 
d 
Sel, ~g 
dv v=0 





Thus the new variables u, v form a local coordinate system in a sufficiently 
small neighborhood of 0. We suppose that in these coordinates we have 


w=) P,du'+ Bdv, P;=Pu,v), B= Blu,»). 
Since w vanishes identically on u = af, v = constant, we have the relation 
| >: Pat, v)a’ = 0. 


To continue, we consider the mapping @ on (f, a, v)-space to (tu, v)-space 
given by 
p(t, a, v) = (fa, v) = (u, v). 
We have 
o*w = ¥ P,(ta, v)(a' dt + tda') + Bita, v) dv 


= by tP,(ta, v)da' + Bita, v) dv, 
o*w = ¥ Pit, a, v)da' + Bit, a, v) do, 
where P,(t, a, v) = ¢P,(ta, v) so that PO, a,v) =0. The important point is 
that o* w is free of dt. 
The equation dw = 9 A w implies d(¢* w) = {¢* 0) A (b*@). We may 


set 
o* 9 = A(t, a, v) dt + other terms. 
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Then 
oP, 
d(o* wo) = > aE dida' + other terms. 


We compare the dtda' terms in this on the one hand and (#* 8) A (¢*@) on 
the other to obtain 


yl 


é = HP,. 


& 


t 


But this combined with P,(0, a, v) =0 implies, by the uniqueness theorem 
for ordinary equations, that P,=0. Hence P, = 0, 


wo = Bdv, 
the desired result. 


References. For this theorem and the generalization which follows, see 
E. Cartan [9, p. 46], [7, p. 367]. 


Example. w=xdy—ydz. Certainly w A dw = 0 since w A dw is a three- 
form. However, the form w vanishes at 0 so one does not expect that the 
integral curves of @ = 0 will span out evenly a neighborhood of 0; in fact 
these curves are just the lines ax + by = 0 through 0. We note, however, 
that dw = 8 A w is impossible in any neighborhood of 0. For dw = 2dxrdy 
so that if 6 = Adx + Bdy, then 2 = Ax + By which fails at x = y = 0. 


Remark. From the theorem we easily deduce again that a one-form w 
satisfying dw = 0 is exact. For consider 0 = dz — w where w is a form in 
x-space. Then d@ = 0, hence there is a one-parameter family z = F(x, c) 
of integral surfaces, F(0,c)=c. For each choice of c, @ vanishes on 
z=: F(x,c), ie, o=dF. (We cannot proceed without passing to one 
more dimension since @ may vanish at 0.) This trick of introducing a new 
independent variable for an unknown function is a useful one. 

We now pass to the general problem. Let w',---," be one-forms in 
r+s space, linearly independent at 0. Set Q=w'a--: Ag’. The 
system is called completely integrable if it satisfies any of the conditions of the 
following lemma. 


Lemma. The following conditions are equivalent: 
(i) There exist one-forms 6, satisfying 
doi = > 6, Aa! (@=1,---,r) (n=r-+s) 
j= 
(ii) dw’ AQ=0 ((=1,--+,7) 
(iii) There exists a one-form A satisfying 


AQD=AAQ. 
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Proof. That (i) implies (ii) is obvious (but unnecessary). Also (i) implies 
(iii) with 4 = 5° 6',. Next, (iii) implies (ii) is the case since (iii) means 
Y(-)dwinga-:: Aa na’ a ++ ar 
r 


=JA@gines* AW 


and we merely multiply by w! to deduce (ii). 
It remains to prove that (ii) implies (i). Let w’*', ---,@" be one-forms so 
that w1,---,@" form a basis of all one-forms. We write 


dw = »y Si ipo? AN o*. 
j<k 


Since dw' A Q = 0, we have 


> fipo avs ao naw aa =0, 
r<j<k 


hence f',, = 0 for r <j < k, 
r n 
do =>) (> — fio) ao’. 
j71 k=j+1 
Frospenius INTEGRATION THEoREM. Let w',---,w" be one-forms in E’, 
n=r-+ 8, linearly independent at 0. Suppose there are one-forms 0°, satisfying 


dw’ = 2 ae (imd,-++,r). 
Then there are functions f',,9/ satisfying 
tile Y fide’ Gi=1,-+-,r). 
j= 
Discussion. The hypothesis is certainly a necessary one. For if we 
write 
wo=(w',---,o0), F=f g=(',°°:.9) 


the conclusion is @ =dg F. The matrix F must be nonsingular in a neigh- 
borhood of 0 and so 


dw = —dgndF=-wna F'dF=on0 


where 
© = —F-'dF. 


Next we note the hypothesis is invariant under a linear transformation of 
the w'. In fact, if y= wA where A is an r x r matrix of functions, non- 
singular near 0, then 


dy=dwA+andA=oarAOA+aadA 
=n A (A7'94 +A 'dA). 
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We shall give two proofs of the theorem, each from a somewhat different 
point of view. The starting point is always the same. We write 


wo! = y hijdai (t= 1,-++,7). 
j=1 


Since the w! are linearly independent at 0, some r x r minor of |h' jl] is non- 
singular in a neighborhood of 0. We multiply (w', ---, aw") by the inverse 
of this minor. On changing our notation slightly then we have 


w= dzi'— >} Al(a’,--+,a%, 21, +++, 2) de! (¢=1,2,-++,7). 
jHil 


First proof. For each point a = (a',--:, a’) in x-space we consider the 
system of equations w' = 0 along the linear variety x = ta: 
dz' 


a > Ai(ta, z)a/ 


with initial conditions z‘(0) = c'. By ordinary differential equations, there 
is a unique solution in a sufficiently small neighborhood of 0, i.e., there exist 
functions F'(t, a, ¢) satisfying 


= (i, a,c) = y A‘ [ta, F(t, a, ¢)Ja/ 
Ot j=t 


Fi(0,a,e)=ci =((¢@=1,°°-,r). 


We shall write F = (F',---, F"). 
Next, we fix k and set G(t, a, c}) = F(kt, a,c). Then G(0, a, c) = c and 


agi 
ot 

hence by uniqueness, G(t, a, c) = F(t, ka, c), i.e., 
F(kt, a, c) = Fé, ka, c). 


i 
(é, a, c) = ‘ (xt, a, c) = y A! (tka, G) ka’, 


In particular, setting = 1 and then replacing k by £, 
F(t, a, c) = F(1, ta, c). 


We pass to new variables u, v according to the transformation 


=U 
z = F(i, u, v). 
This is nonsingular in some neighborhood of 0 since 
0(x, Z) _|f Of jez) 
(u, V)fo Ox ~ |avig 











* OV, 
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For 
dz! - 
dv | 


In these new variables we may write 


Gv! 


| a! 


i? 
0 


y i 
F'(1, 0, v) eg ere F 








= y Bi(u, v)dv* + y P' (u, v) du. 
k=1 j1 


The fact that each w' vanishes identically along the curve u = ta, v = con- 
stant implies 
> Pi (ta, v)a/ =0 (@=1,---,7r). 

jJ=1 
We propose to show that the functions P',(u, v) vanish identically. To do 
this we consider the cone mapping @ on (f, a, v)-space to (u, v)-space 
defined by 

o(t, a, v) = (ta, v) =(u, v). 

We have 

o* w' = > P' (ta, v)tdai + terms in du* 


= ) Pit, a, v)da/ + terms in dv* 
where P'(t, a, v) = P',(ta, v)¢ so that P',(0, a, v) =0. It follows that 


d p*w' = = =! dtda/ + other terms. 


Finally we use the hypotheses dw’ = } 6', , w*. We write 
* 6', = H'(t, a, v) dt + other terms 
and compare the coefficients of dtda/ in the relation 


d p* oo’ = )) (p*8') 0 (G* w'): 


te (t,a,v) = > A(t, a, v) P(t, a, v). 


We conclude from the uniqueness of solutions of ordinary systems together 
with the initial conditions P! j(0, a, ¥) = 0 that Pi ;=9, P: j;= 9, 


=> Bi,(u, v) dv 


as required. 
Our next proof is based on the sketch in E. Cartan [10, pp. 188, 193]. 
We begin as before with the system 


wiad?i— Y Adri (b= 1,++*,7) 
jF1 
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with the conditions dw! = ¥ 6’; A w/. We take any smooth curve from 
the origin to a point a. We solve the system w! = 0 on the cylinder this 
curve spans in x, z-space, taking some definite initial point c on the z-axis. 
We shall show that the point on this curve lying over x = a is independent 
of the particular curve we start with in x-space. The point is that in a suffici- 
ently small neighborhood of 0 in x-space, any two smooth curves with the 
same end points 0, a can be smoothly deformed, one to the other. Thus let 






XxX = x(t, 0) 
x = x(t, 1) 


x = x(é, a) be a one-parameter family of curves from 0 to a, the time variable 
t on each curve running from 0 to 1 and the parameter a taking all real values; 
we are interested in the curves x(t, 0) and x(t, 1). We are assuming 


x(0, a) = 0, x(l, a) =a. 


Fixing «, the solution of w' = 0 on the corresponding cylinder with initial 
value ¢ is given by functions F'(t, a) satisfying 


* => A! (x(t, «), F(t, 0) = 


Pe (¢=1,-+-,r). 
We reduce the problem to a two-dimensional one by considering the mapping 
@ on (é, «)-space to (x, Z)-space given by 
p(t, a) =(x(E, «), F(t, «)) = (x, z). 


Then 


o*w a as oF dn — At a 4 


= H' da 
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where 

; aF' dx! 

H' = Hit, a) = — — ¥ Ai, —. 

tH Gor LA; da 

We set 
p* 0, = P',dt + Q' da 

and compare coefficients of dt da in 


dp*w' = ¥ (p*6',) a (b* @!) 


to obtain 
oH' 
— = P'. Hi, 
ot Bey 
But 
CF’ Fala 
H4(0, a) = — — ¥ 4*.(0, a) — 
(0, a) Be ee j(O, a) Bat | smo 
d _. dxi(0, «) 
= — F'(0, — > A‘(0, a) ———— 
Ta (0, a) y) A(0, a) a 
a . d0 
=—-¢' — § A'(0, a) — = 0. 
fe 24 0,9) da 
It follows that H' = 0, 
OF' . ae 
—(t,a)= ) A’.—. 
ale (a) = D4, Ca 
We apply this in particular at t = 1; here 
dx! ad. d 
oe t, = — x l, = —qi = 0, 
Ou i, a) i=, do th, a da. 
hence 
OF' 


—— (1, a) = 0, F'(1, «) = constant. 
Ou 


We next fix the notation more precisely; we write F'(t, a; a,c) instead 
of F(t, «) so as to specify the dependence of this function on the initial 
conditions. Since F'(1, «) is independent of a we may set 


Gi(a, c) = F'(1, «; a, ©), 
and also 
F=(F',-.-,F), G=(G!,---,G@). 


Then we have the following facts: 

(i) G(0,c)=c 

(ii) For fixed ¢, a ——> (a, G(a, ¢)) is a 1-1 correspondence on a neighbor- 
hood of 0 in a-space onto a manifold V, in a, ¢-space. (For we simply take 
any curve from 0 to a and use it to define F and then G.) 
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(iii) Each @‘ vanishes identically on V.. (For w! vanishes on each 
curve on V,, since it is a one-form it vanishes identically.) 
We consider the mapping 


(a, c) —> (a, G(a, c)) = (x, z) 

on a, ¢c-space to x, z-space. Because of (ii), 
O(x, Z) 
O(a, c) 
hence we may use (a, c) for a new coordinate system in some neighborhood 
of 0. Writing w! in these new coordinates and using (iii) shows us that w! 
involves only the differentials dc’, --- , dc’, which completes the proof. 

The striking feature of these proofs is that we reduce the original system of 
partial differential equations (with integrability conditions) to a system of 
ordinary differential equations. 


— as 
0 








7.4. Applications of the Frobenius Theorem 


Example 1. We begin with a question in matrix form which is motivated 
by considerations in differential geometry centering around infinitesimal 
transformations. 

Let Q = ||~@,’|| be an r x r matrix of one-forms defined in a neighborhood 
of 0, say, in E". We ask when it is possible to find an r x r matrix A of 
functions, nonsingular, satisfying 

Q=dAA“}, 
To fix matters, let us require the initial value A, =J. It is convenient to 
set 

© = dQ — 07. 
Then the basic result is this. 

There is a matrix of functions A defined in a neighborhood of 0 such that both 

Q = (dA) A! 

af and only if 
© = 0. 
When this 1s the case, then there 1s only one such matrix A. First of all suppose 
there is a solution A. Then (dA) A~' =Q, dA = QA, and we have 
0 =d(dA) =dQA —QdA = (04+ Q?) A — Q(QA) = OA. 
Hence ©A=0. Since A is nonsingular we have OQ =0. If B is another 
solution so that dA = QA, dB = OB and A, = B, = I, then 
d(B-'A)= —B-'dBB"'A+B™'dA 
= —B~'(QB) B-1A + B71(QA) =0, 

hence B~'A is constant, B™1A =(B7!A),=J1, B= A. 
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Now we come to the existence. We pass to (n + r”)-dimensional space 
with cqordinates x!,---,2",z/ (1<1,j <r) and introduce the r? forms 
which are the coefficients of the matrix 


A =dZ — QZ, Z = \lz,)|j. 
We are assuming © = 0, hence we have 
dA = ~dQZ+ 0dZ = —-07?Z4+ Q(A + QZ), 
dA = QA. 


It follows that our system A, which is already in standard form, is completely 
integrable, hence there exists a matrix A of functions of x with prescribed 
initial values at x = 0, so that Z = A is an integral manifold of A = 0, that is, 


dA = QA. 


We remark that if QO is skew-symmetric, then A is orthogonal (provided 
the initial condition is Ay = I). For we set B = ‘A™', the inverse transpose 
of A, and have 


By = I, adB= —B(d'A)B= —B'A'OB= + QB. 
It follows by uniqueness that B= A. 


Example 2. We next consider another equation: 
dA = QA — AQ. 


Here Q is the same as before, anr x 7 matrix of one-forms in a neighborhood 
of E” and A is the unknown matrix of functions. Again we set @ = dQ — Q?. 
Let us pose the problem this way. Can we find a solution A of the above 
equation taking an arbitrary initial value A,? We seek a necessary con- 
dition by differentiating: 


0 = d(dA) =dQ A —-QdA —~dAQ-—AdQ 
= (© + Q?) A — Q(QA — AQ) — (QA — AD) Q — A(O + 22), 


which simplifies to 
OA = AO. 


Since we are assuming the initial values A, may be arbitrarily prescribed, 
the values of A at each point of a sufficiently small neighborhood of 0 will 
fill an n-dimensional domain; the commutativity of © at such points with 
so many different A evidently implies that the matrix © of two-forms must 
be of type 
O=al 
where a is a two-form. From 
dQ —- 7 =al 
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we have, differentiating, 


dal = —dQN + dQ 
= —(al + 27)2 + Olal + 22) =0. 
Thus, locally, « = do where o is a one-form. The necessary condition we 
arrive at is this: There must exist a one-form o satisfying 
dQ — Q? = (do) I. 

This can also be expressed another way: The matrix 

H=Q-ol 
satisfies 

dH — H* =0. 


Under this condition, the sufficiency is easily demonstrated. As before 


we form 
T=dZ -—-0Z+4+ 20 


in X, Z-space, and note that 

dV = -dQZ+4+0dZ+dZ0+4+ Z2d0 
= —(0? + 0)274+ QT + OZ — ZQ) +(T + OZ — ZQ)N + Z(Q? + 60) 
= Z0 —- 024+ QF + 1Q, 


hence 
aY = OT + TQ, 


which shows that the system [is a completely integrable one. The existence 
proof now proceeds as in the Jast example. Uniqueness can be handled this 
way. Since the system I = 0 is in normal form and is completely integrable 
we know there is a unique integral surface passing through a given initial 
point Z|9 = Ao. 


Example 3. We shall consider a type of system of partial differential 
equations known as a system of A. Mayer (see C. Carathéodory (6, pp. 
26-31)). 





We work in a neighborhood of 0 in Et’ with coordinates z',---, x°, 
z', ---+,2" as before and are given functions B' (x, 232 1),2 57,9 = 108° 58 
The Mayer system is 

Gz! ; 
aa = B j(®, Zz). 
We define 
6B, dBi, OB; B eB 
A Se ee a BP. 
Ke 6x" Oa +2 az" lar 
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We evidently have A‘, +A',;=0. The Mayer system is called completely 
integrable in a neighborhood of 0 provided to each choice of initial conditions 
c there exists a solution z = F(x, c) of the system with F(0,c)=c. The 
necessary and sufficient condition for complete integrability is precisely 
A‘, = 0. The reason is the following. We set 


w'=dz'— Y Bi(x,z)dei (it =1,-+-,7) 
j=1 


a system of one-forms in x, z-space in our standard form. The vanishing 
A! jk = © implies, after a short calculation, 


do! = ¥(E5 OB 


so that the system qw',:--,qw’ is completely integrable. The integral 
surfaces z= F(x, ¢) solve the Mayer system. Conversely, suppose the 
Mayer system is completely integrable. Then it is clear that the system 
w! = +--+ =q@"'=0 has integral surfaces, one for each choice of initial 
conditions c. Hence the necessary condition 


dw = y gf, A w! 
must be satisfied. On the other hand, we directly verify that 
=4) Ai, daida + ¥n', 0 w 


ac!) no 





where 
n= ek del. 
Since dz!, +--+, dx’, w'!,--+, a" are linearly independent we conclude that 
Oj Aw =) ni,nw, YY A',deide*=0, Aly= 
Note. If w',---+,w" is completely integrable, there is an r x 7 matrix of 


one-forms © = ||6',|| satisfying 
dw' = YO, na’, 
or 
doa = -waAO 
in matrix notation. However, from the solution 
w=dgF, g=(9',--'.97), F=|fiyll 
we conclude that dw = —@ A F~'dF so that we may always choose © in 


the very special form © = F~'dF. (This provides some motivation for the 


Problem in Example 1.) 
We shall see further applications of the Frobenius theorem in our study of 
local Riemannian geometry. Also cf. Problems 4-7, p. 194. 
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7.5. Systems of Ordinary Equations 


We consider a system 


dx! 
— = X(t, x}, <2) 
dx" 
— = X"(t, v1, ---, 2"). 
dt ue 


Closely associated with this system is the differential n-form 
Q. = (dx! — X' dt) --+ (da" — X" dt) 
in (¢, x)-space. By a short computation, 
dQ = (x =| dtda' --+ dz". 
We make a change of variables 
y' = y'(t, z', +++, a) (@=1,-++,n) 
and suppose that the systems 
= = X'(t,x) and wv = Y't,y) 
are equivalent under this change. We set 
Q = (dy’ — Y' dt) --+ (dy" — Y"dt) 
and propose to determine how dQ and dQ are related. Now 


dy’ _ dy’ dy" da! 
dt at Oxi dt’ 
dy! oy 
y' = — 
at Bal 
Also 
dy' = oy" — dt + 5 Hae 
“at 
hence 


; dy’ 
dy’ — Y'dt = ) — (dx’ — X/dt). 
y » 53 (de ) 
We denote the Jacobian by J: 


ay',+*+,y") -|3| 


J= 
O(xt,--+, a") | da 
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and have by exterior multiplication, 
Q) = JQ. 


We differentiate this: 
aQ = (dJ)Q +7 dQ. 


Now 
oY! 
dQ = (E55) atay! «a . y= (J rs a) ata! -+ da" 
and 
J 
(dJ)Q +SdQ = (7 d+ Dade ‘Na+a(y 5) atae' -o + dz 
od Os _.. ox! 
= — —— J — I... 
(5 + ai & + » =) dtd dz", 


hence 
ay’ lj _ ax 
Lag ala t lar} 


A function f = f(t, x) is called a first integral of the system if f is constant 
along each trajectory, or solution curve. Since the usual existence and 
uniqueness theorems guarantee a solution through each point of space where 
the system is defined, we have the condition 


df 


= 0), 
dt 


A 
i 7 rae aaa 


for a first integral. 
Suppose that each of the functions y', --- , y" in the transformation above 
is a first integral. Then 


Yi=0, Q=dy'--- dy’. 


(Such a transformation is always possible in a small region of space because 
of the existence of a general solution, one depending on arbitrarily prescribed 
initial conditions.) 
A function M = M(t, z',---, 2") is called a last multiplier of the original 
system if 
ad(MQ) = 0, 
1.e., 
OM a(MX"') 


es == {), 
Praia se ger ae 
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Using the transformation above based on first integrals, 


MO = (F)R= HO = Hay! --- ay. 


Hence 


H 
d(MQ) = (dH)dy' --- dy" = ded! . + dy" 


so that M is a last multiplier if and only if H is independent of f, i.e., 
MQ = Hy) dy' --- dy”. 
If M, and M, are two last multipliers, then 


M ,Q= H,(y) dy’ --- dy’, M,Q= H,{y)dy' --- dy”, 
hence 


M ,/M, = A,(y)/H2(y). 


It follows that M,/M, depends only on the y', hence is constant along 
trajectories (as is each y') and consequently is a first integral. This proves 
the important result: 

The quotient of two last multipliers is a first integral. 


7.6. The Third Lie Theorem 


What is known as the Third Fundamental Theorem of Sophus Lie was 
devised in order to reconstruct a continuous group given only its constants of 
structure. These concepts will be explained in Chapter IX. For our present 
purposes, we shall look upon this theorem as a result, and a rather deep one 
at that, in partial differential equations. 

We work in E”. All indices run from 1 ton. First of all we are given n 
constants c',, subject to these constraints: 


3 


C' in + c¢,;=90 
ial i J a ee 
D(C Ons + jr Ook +O jsC kr) = 0. 


Fi 

The problem is to find n one-forms a’, ---, o” which are linearly indepen- 
dent on some neighborhood of 0 in E” and which satisfy the relations 

dai = 4 > c 4.07 A of, 
The Lie Theorem asserts that this can be done. 

The quadratic relations we have assumed for the constants (c) are easily 
verified to be the same as d(do') = 0, assuming our problem is solved, hence 
they are necessary conditions. That they also are sufficient conditions will 
now be seen. The proof we give is based on those in Cartan [10, p. 239) and 
[7, pp. 280-283]. Because the proof is lengthy, we shall break it into several 
steps. 
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Step 1. We define an n X n matrix F = [f',] of homogeneous linear 


forms by 
i. by cy, D. 
Then we consider the linear initial value problem 


oH = 1+ HR H(0,x) = 0 


for an nXn matrix H=H(t,x). For each x it has a unique solution, 
defined for all ¢, and the solution is analytic in (¢,x). At x=0 we have 
F=0, so 


H 
= (0)=1, H(0,0)=0, 


hence 
H(t,0)=tl and 4A(1,0)=I. 


Step 2. We set 
w= dxH (t,x) 
so that @ is a row vector of one-forms, free of dt. Clearly 
dw=dt dx(I+HF)+z 
=dt(dx+wF) +A, 


where A is a row vector of two-forms, also free of dt. We also define a square 
matrix A=([a',] of one-forms by 


a= >ci,0" 
We note the obvious relations: 
xA=OF and dxA=—odF. 
We also note two less obvious relations : 
d(wA)=2dw AA and 2x4°=0AF. 


The first follows from the skew-symmetry of the constants c’,, in their lower 
indices. The second follows from the quadratic relations on the c’s. We 
multiply the quadratic relation (p. 108) by z’w*w* and sum: 


Ds mee. cp Bese 
Di graiai+ > fici,o wit > xa jo’, = 9, 


which in matrix language is —xA*+@AF —xA*=0. 
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Step 3. From the formula for dw that defines A, we have 
dk =dt A d(wF) =dt(dwF —wdF) 
= dt(AF—odF) =dt(AF + dxA). 


Step 4. We define 
0=A-— 4A, 


a row vector of two-forms, free of dt. We shall prove the decisive formula 


d@=dt A@F-AAA. 


We have 
d0=ddA—dw A A=dat(AF + dx A A) —[dt(dx+oF) +A]A 
=dt(AP—wFA)-AAA 
=dt(AF—-xA’*)-AAA 
=dt(AF—ijwAF)-AAA 
=dt \OF-AAA. 


Step 5. We shall now prove that 
6’ = 0. 
Since both the w/ and the ! are free of dé, so is 9’ and we have 
P=) Go ndalde, gi, =gialt, x). 


We may even assert that 
g" (0, x) = 0. 
For h' (0, x) = 0 which not only means that the coefficients of the w! vanish | 
at t= 0, but also that the coefficients of 
Oh’, Ohi, 
_ —i - da* dai 
v=4) ( aa ae) aa) . 
vanish at ¢ = 0 according to the very definition of partial derivative, 


Ohi, _ hi (0, x) x) 
oak rol OH 


From these facts, what we say about the initial values of the 9’ jx follows. 
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The result of Step 4 implies 


09" i r Ss 
at = elt Pix 





This homogeneous linear system taken together with the initial conditions, 
the vanishing of the g’s at t = 0, has the unique solution 


; oF ixlts x) — 0, 
and so @' = 0. 


Step 6. Now we can wind up this story. Since 9‘ = 0 we have 
rl — 4 » ca A a*, 
doi =4)c,0/ nwo +dtna. 
We consider this relation on the subspace ¢ = 1. 
Setting 
f= o'l|,2, = > Ail, x)dz’, 

it becomes 

do'=}4) ci,o/ a ot. 


Since h',(1, 0) = 5; (Step 1), the one-forms a", --- , o” are linearly independent 
at 0, which implies of course that they are linearly independent in some 
neighborhood of 0. The proof is complete. 


Vill 


Applications to 
Differential Geometry 


8.1. Surfaces (Continued) 


Everything in this section will be based on the local theory of Section 4.5. 
Now we have integration at our disposal and we shall discuss a few global 
results. Let Z bea closed surface in E*. For e, we take the outward drawn 
normal to &%. The mapping 

x—e, 


is a map on £ to the unit sphere $7. As x varies over L, e, varies over S$? a 
whole number of times, called the degree of the normal map (cf. Section 6.2). 
The element of area of the normal map is 


WO; = Ko,0, 
since 


Here K is the Gaussian curvature. Hence 
| K oO 19 i 4nrn 
r 


where n is the degree. The factor 47 is simply the area of the unit sphere. 
In particular, if £ is a closed convex surface, then e, covers $* exactly 
once as x covers LZ, hence 


| Ko,0, =4n 
r 


in this case. 
After this, we shall limit our discussion to closed convex surfaces. Two 
important invariants are the total area 


A -| 0104 
b 


and the integrated mean curvature 


u =| Ho,o2. 
b 


112 


8.1. SURFACES 113 


Given a closed convex surface Z and a fixed positive number a, we form 
the surface L’ parallel to X at distance a by marking off on the outward-drawn 
normal at each point x of X the distance a and taking the locus of all points 
so obtained. Thus the typical point on the parallel surface is 

y=x+ae, 
where e, always denotes the normal at x. We have 
dy =dx + ade, 

= (oe, + o2€,) + a(wye, + @2e,) 

= (0, + 4m,)e; + (02 + aw2)e,. 
It follows that the normal to the parallel surface L’ at y is again e, and that 
e, and e, can be taken as a basis of the tangent space at y. Thus we have 
with 

T =O, +90,;, 12> 02+ 402. 
It follows that the element of area of LX’ is 
™T2 = (o; + aw,)(62 + a>) 

= 06,0, + a(o,@, — 0,0,) +@7a,0, 

= (1 + 2aH + a*K)o,0, 
so that the total area of 2’ is 


A’ ~ | a -{a + 2aH + a? K)o,0, ; 


A'=A+2aM + 4na*. 


(This formula can also be proved by first doing it for a polyhedron and then 
taking limits in an approximation of % by a sequence of polyhedra. If one 
examines what the formula means when Z is.a convex polyhedron, one will 
see that H measures dihedral angle and K vertex angle.) 
By integrating with respect to a, one easily comes to a relation between the 
three-dimensional volumes V’ and V enclosed by 2’ and &, respectively: 
V'=V+aA+a°M + 4na°. 
One can also verify the relations 
M'= M + 4na, 
_ H+ak 
~ 1 +4 2aH +a7K’ 
ae K 
«1+ 20H +a?K’ 


é 
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We next introduce the support function of our closed convex surface 2%. 
This is defined by 
Pp = X-e 3° 


It is convenient to fix Z in space so that the origin O is inside Z. Then we 
have p > 0 at each point of L. 
The following method will be used to obtain several identities. Let 1 be 
any one-form on Z. Then 
| dj = 0. 
£ 


For dX = 0 and Stokes’ theorem gives us 


[a-| A= 0. 
x ax 


First we consider the form 


a= e,:(x x dx). 


Here 
da = de,:(x x dx) + e,:(dx x dx). 
Now 
de,:(x x dx) = —x:(de, x ax) 

= —x-(dx x de,) = —x:(2Ho,a,e;) 

= —2Ho,0,(x-e;) = —2pHo,0, 
and 

e,°(dx x dx) = e3-(20,0,€3) 
= 200% , 

so that 


da = 2[0,0, — pHo,o,]. 


Since the integral of da is zero, 


A -| 01,04 | pilose, 
» > = 


Bp=x-(e, x de;). 


Next we set 


By a calculation similar to that used for a, 


ap = 2[pKo,o, = Ho ,o2] 


M -| Hoo, | pk. 
> z 


so that 
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Since we have found the integrals of H and K weighted by 7 it is also reason- 
able to seek the integral over X of the form po,c,. We get this by starting 
with the vectorial area 


(9,0,)e, = 4dx x dx = (dydz, dzdz, dxdy) 
from which 
pa ,0, = (x-e,)(0,0,) = x: (dydz, dzdz, dxdy) 


= xdydz + ydzdx + zdrdy. 
Let R be the region of E* bounded by the closed convex surface and let V 


denote its volume. Then 


| P0104 -| (xdydz + ydzdx + zdxdy) 
(x=éaR) 


x 


-| d(xdydz +--+) = I) dxdydz = 3V, 
R R 


:| po,o, = V. 
E 


We close this section with the following interesting theorem: 


Let X be a closed convex surface of constant Gaussian curvature K. Then 
LX 1s a sphere. 


To prove this, we recall the relations 
dx = 0,e, + 02e, 
i = WC, + We, 
OO, = po, + qo2 
' = Oy 10; 


ae 
K=pr-q’ 


developed in Section 4.5 for any moving frame. Since Z is convex, the 


matrix 
P 
q ?r 


is positive definite, p > 0,r > 0, K > 0. (See p. 120-121 in the next section 
for details.) “Because of the arithmetic-geometric mean inequality we have 


K=pr—q Sprs(i(p tr)? =H’. 
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We also note that there can be equality, K = H?, only when ¢g=0 and 
p =r =H, which implies w, = Ho,, w, = Ho,, de, = Hdx. 


But 
{| 0,05 {| pHo,o,2 {| pJ/Ko,o, = =z || pKo,0, 
= 3 {| Ho,o,= =| /Ko,0, = {| 0185; 


where each integral is taken over X and we have exploited the hypothesis 
K = constant. Because the quantities at the ends of this chain of in- 
equalities are equal, all integrands must be equal, H = ./K. By our remarks 
in the last paragraph, this implies that de, = Hdx with HA constant, 
Hx =e, +Hec with c a constant vector, |x ~— c| = (1/H)/e,| = 1/H, Lisa 
sphere. 


8.2. Hypersurfaces 


We shall extend our study of surfaces to higher dimensions and at the same 
time motivate some of the things in the next section on Riemannian geometry. 

A hypersurface is an n-dimensional manifold M embedded in E"*!. We 
denote the moving point on M by x. Our study is local so we pick a definite 
unit normal n at each point x of M. The map x — n is a smooth map on 
M into S". (This can be done globally on a hypersurface M precisely when 
M is orientable.) The tangent space at x is an n-dimensional Euclidean 
space; we pick an orthonormal basis for it, e,,---,e,. Thus at x, the 
vectors e, ,°::,¢,, 1 make up an orthonormal basis of E"*!. Since dx is 
in the tangent space we have 


dx =0,e, +s: + 0,é,, 
where ¢,,°°:, 06, are one-forms on M. From the relations 


ee, = 0%; e,n= 0, n-n= 1, 
we deduce that 
de,-e, + e, de, — 0, 


de,;;n+e,;-dn=0, n-:dn=0 
and so 
de, = y @;;e; aa m,n 


dn = >. je; 
where (;;, @; are one-forms on M and 


Oi; + O jj = 0. 
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It is convenient to write all of these structure relations in matrix form. We 


set 
e; 
e= : ; 6 =(0,,°°°,9,), 
e,, 
Q = lia; |, w= (@,,°°*,@,). 
Then 
adx = ae 
t 
()=(o Wo )(a) 
n (1) 0 n 
, Q+'N=0. 


By taking exterior derivatives we obtain integrability conditions. (We 
shall omit the symbol “a” in what follows. All products of differentials are 
exterior.) 

0 = d(dx) = (da) e — a(de) 


= (do) e — 6(QMe — ‘wn) 
= (do — cQ)e + o'on, 


: s(2)] =62, «a'w=0; 
(2 S90)-@ “oKG 
(2 E4Q)-@ 70 


_ (dQ-2’? +'ow —‘(dw) + Q'w\ /e 
~ \dQ- a 0 n}’ 


dQ-027+'aw=0, dw=an. 
We define a skew-symmetric matrix of two-forms: 


© = ||0,,| = dQ — 02. 
We sum up our results: 


do = 6Q, 
Q + 'Q=0, 
o'o =0, 
da = aQ, 


0 + ‘aw = 0, 
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or in terms of individual elements of the matrices, 
do; =) 6;0;;; 
Oj; + on= 0, 
Yi o,0;=9, 


dw; = > O;@;;; 


6: + O;O; = Q. 
The o; form a basis for one-forms on M, hence we have relations 
; = pa bij CG; . 
Because )' o;w; = 0, the b,, must be symmetric, 
b;; = bj; < 


The mean curvature H and Gaussian curvature K are defined by 
j 
=~ bu, K = (6,,| 


Since o, °:: 0, is the n-dimensional volume element on M and a, -:: a, is 
the corresponding quantity for $", K represents the ratio of volumes, volume 
of spherical image over volume of M, due to 


Oy + Dy = (Yj 05) °° (Lng os) = [Bile °° 
= Ko,'--6,. 
Suppose one has a function v = v(y, z, -: +) of several variables where v is 
always a tangent vector to M. For example, we might assign to each point 
of a curve on M a tangent vector at that point, arriving at a vector-valued 


function of one variable. 
How does an observer constrained to M observe the motion of v? We 


write 
v= > c,e; 
where the c; are functions and have 
dv = y dc;e; -+- > c,de,; 
= )ideje; +) ei(), oie; — on) 
=) (de; + di eimis)e; — (L ea,)n. 
Our observer who is constrained to move in the hypersurface M cannot 
“see”? the motion of v which takes place in the direction normal to M; he 


sees only the tangential motion of v. Consequently he believes v is motion- 


less provided 
(de; + ) c,a,;) e; = 9; 
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that is, 
de; + ¥ ¢,a;; = 0 (G=1,:--,%). 


A vector function for which these equations are valid is said to move by 
parallel displacement. 

The following can be checked. If v = v(y,z,---) and w= wly, z,--:) 
are two such vector-valued functions which are compatible [for each point 
(y,Z,°°*) in the parameter space v and w are tangent at the same point of 
M] and each moves by parallel displacement, then v-w is constant. In 
particular, |v)? = v-v is constant. 

Let P = P(s) be a curve on M parametrized by its arc length s so that 


dP 


is the unit tangent vector. The curve is called a geodesic provided t moves 
by parallel displacement. 

There is a geometric interpretation of the matrix |jb,,|| which is quite 
fundamental. To each particular displacement dx of the position vector x 
corresponds a displacement dn of the unit normal n. Both dx and dn are 
in the tangent space at x so that we can look on 


dx — > dn 


as a linear transformation A of this tangent space. More precisely, let v be 
any tangent vector at x. Pick any curve x = x(¢) through x so that 


dx 


di = V. 


0 





We follow the normal n = n(¢) as x traverses the curve. Then 


dn 
—| =Av 
dt | o 
is our definition of A. We see that this is quite independent of the choice of 
the curve x(t) so long as it has the prescribed tangent v at ¢ = 0. 
For suppose 








v=c¢,e, +--+ +6,¢4,. 
Then 
dx =) a;@; 
SO that 
Gj 
dtl, °° 
But 
an QO; 
—_ = —| @e.. 
dt |\o dt | 
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Now 


= 2 buss Oj» 


alt oe by dt}, 


A(d c€;) = D() b,,¢;)e; 
which establishes these points: (1) A is a well-defined function on the tangent 
space at x to itself, (2) A is linear, (3) the matrix representation of A with 
respect to the basis e is |[0,,|\. 
Since the matrix ||0,;|| is symmetric, the linear transformation A on the 
(Euclidean) »-dimensional tangent space at x is self-adjoint: for each pair 
of tangent vectors v, w, 


=) b4;¢; 





so our result is 


(Av)-w = v:(Aw). 


It is clear that our definition of A depends only on the hypersurface M and 
the way it is embedded in E”*', not on our choice of the moving frame e. 
Consequently the formulas 


pi 


H = — trace (A) 


A| 


show that H and K are geometric quantities. 

Since A is self-adjoint, its characteristic roots are all real and they are 
called principal curvatures. The corresponding characteristic vectors define 
(in general) » direction fields on M called principal directions whose 
integral curves are lines of curvature. 

Convexity of M can be interpreted in terms of a definiteness condition on 
the transformation A. To make this precise, suppose M is convex and we 
choose for n the inward unit normal. 


| 


K= 
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0 


x 


We fix a point x) on M and form a corresponding normal section, the curve’ 
of intersection of M with any (two-dimensional) plane on the line ny. This 
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is a curve X = x(t) which is convex on its plane, which means that the function 


(X — Xo): 





is convex, hence satisfies 





d*f 
— = 0. 
dt? t=0 - 
We have 
df (dx(t) _ a; \. 
a7 ( dt m= (5 ) No: 
d*f d (a; a,\ (de; 
eo r5(a)e Mot) Bley No. 
Now 


pee i) nn 
a = E (Fe eee (3 i 


Since the tangent vectors e,(0) at x(0) = x, are orthogonal to the normal 
n(0) = no, the condition reduces to 
iz) 
o\dt]|o 
< 0. 


= (3) 
(a), 


But OO; = ) bi;9;; SO we have 
(=: . 3) 
it } 3 it js 


Gi 
b,:1— 
» J ( 4 
is arbitrary, we conclude from this that the matrix 
HO ssl 


is negative semidefinite so that the same is true for the transformation A 
which this symmetric matrix represents. 


<0. 














Since the direction 





122 VIII. APPLICATIONS TO DIFFERENTIAL GEOMETRY 


(The correctness of sign can be verified in the simplest possible case, that 
of a convex curve in the plane. The usual Frenet formulas are 


dx=dse, 
de,;=xKdsn 
dn = —kdse, (s = arc length) 


so that 0, =ds,w, = —Kds,b,, = —K $0.) 





We return to the general situation. The elements 0,,; of the curvature 
matrix © are curvature forms. We may write 


6,,= 4 >. R ijt On Ors Rin + Biju = 9, 


defining the Riemann curvature tensor R;,, of the hypersurface. Because of 
the relations 


0; + W;; = 0, 
and 
O;,O; = » bi 51 OF: =} Y (64.5, oe 5:15 jx) OG, 5; 
we have 
b. ob. 
R, +4 tk il = 0. 
gk b ig b 








Algebraic consequences of these formulas are the following: 
Ria t Rijn = 9, 
Raj + Ryn = 9, 
Rigg + Bar + Raj, = 9%, 
Be kt = Ryriys 


all of which follow easily. 
We shall see that the Riemann tensor is independent of how M is embedded 
in E"*! so that these relations are particularly interesting, connecting the 
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intrinsic Riemann tensor with the quantities b;,, which clearly depend on 
the embedding. 

Indeed, it turns out that the R’s are determined by the o’s alone with no 
reference to the normal n. This means that if two hypersurfaces M, and M, 
are in a one-one correspondence which preserves distance (i.e., preserves the 
Euclidean geometries in each pair of corresponding tangent hyperplanes), 
then M, and M, have the same Riemann curvature tensor. 


What we shall show is that the equations 
de = aQ, Q+'XO=0 


ij 


determine Q uniquely, so that © is completely determined by the o’s alone. 
This is more in context in general Riemannian geometry so we shall postpone 
the proof until p. 129 of the next section. Having this result, it follows that 


@ = dQ — Q? 


is also completely determined by the a’s, hence so is the Riemann tensor. 
We shall now take up the special case in which our hypersurface M is 
given in the form 


u = u(z',s++, 2") 


in z!, +++, 2", u-space. 
Uu 
n 
x 
x" 
x! 
It is convenient to set 
Ou 074 


Pim ag TU Balog! 


We have for the position vector 


x = (z!,a7,-°-,2",u) 
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and so 
dx = (dx', dz*,---, dz", du) 
= (dx',-+-, dz", > p,dz') 
=) t,dz' 
where 
t; = (011,°°* + Sins Pj): 
The vectors t;,--:, t, are tangent vectors and evidently form a linear basis, 


but generally not an orthonormal basis, of the tangent hyperplane. 
The vector 


w=(-p,,°°°, —p,,1) 
satisfies 


wt, = 0 


so it is a normal vector (with positive component in the u-direction). The 
unit normal n is given by 
w=wn 
where 
w=wew=1+) p,. 
We note that 
wdw= > pdp,=> p,r,,da’. 


We shall now determine the matrix representation of the basic linear trans- 
formation A with respect to the basis t,,---, t, of the tangent hyperplane. 
This matrix is not symmetric in general because (t,;) is not an orthonormal 
basis. However, its trace and determinant are the trace and determinant 
of A since any matrix representation yields a valid determination of these 
quantities. 

Suppose the matrix we seek is |a;;||._ Then this means 


At; =) aj;t,. 


Let v be any tangent vector. Because of these relations plus the fact that A 
is self-adjoint we have 


(Av)-t; = (At,)-v =) a,,(t;-v). 
Now going back to the very definition of A, symbolically, 
A: dx—dn, 
we see that this relation means that 


(dn)-t; = }) a,,(dx-t,). 


8.2. HYPERSURFACES 125 
It is from this set of equations that we shall determine ||a,,||. We have 
dw =dwn-+ wdn, 


dw-t; = w(dn)-t; 
since n- t; = 0. But 


dw :t; — (—dp, ar) —adp,, , O)* (Sis ea Sins Pi) 
= —dp, = -> r, ax 
so we have obtained 
l ; 
dn-t,= —— dx, 
nt; = ry 
On the other hand, 
dx:t; = (dz',---, dz", du): (55; oe Din» Pj) 
= > (Oj, + Pj Py) A", 


so we have 


l 
Y 4555 + p;p,) dx* = — = Linde, 


l 
Y A; (Oj, + Dj Pr) = — . T ik: 


Pi 
R=|r,l, p=[i], A=leyl, 
Pr 


Setting 


we may write this as 


1 1 
A +p'p)=—--—R, A=—- RI +p'p). 
Ww Ww 





Since 
P'p=)p; =v -1 
we see that 
1 1 we —1 
(I + p'p)(z— (=) 'P) =I+p'p- (—3)P'P meagre aah Oe i 
hence 


1 
+p'p)' =1- (=]p'e. 


«tale (S)o) 
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The mean curvature H is found by taking the trace: 


] 1 ] 
H=-tr(A)= - | R-— Rp'p| 
n nw w 


3 —|5 ran > (—3)z.r Ps] 


] 1 
= —|A = y = Par yPy| 


nw 


where Az is the Laplacian of w. 
The Gaussian curvature K is found by taking the determinant: 


(—1)’ 








K=(|Al= RilT si 
|A| = — [RIE + Ppl 
One finds by a short calculation} that 
[I + p'p| = 
so that 
(—1)” 
K= rea | RI. 
For the special case n = 2 of surfaces, we use the standard Monge notation 
_ Ou _ Ou aot ae 07u , ou 
Pox" a By’ x?’ Ox dy’  6y? 
and have 


w? = 1+ p* + q?, 


=i] 
H = 53 [wel + t) — (p?r + 2pqs + q7t)] 
—1 2 2 2 2 
= ol +p° +q°)(r +t) — (por + 2pgs + g't)| 


—l 
9,33 [(1 + q7)r — 2pgs + (1 + p”)é], 


and 
_ (rt — 8?) 


? 
w* 


K 


both familiar formulas. 


ft Set e; =(1,0,...,0), eg =(0,1,...,0),..., @, =(0, 0,...,0, 1). Then 


[J + p*p| = |e: + pi'p,..-,@n + pn ‘pl 
== |@1,...,@n| +> Pe-ler,..-, Ct-1, Pp, Ci41,- ++» On| 
=1+)> p? = wi. 
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8.3. Riemannian Geometry, Local Theory 


The problem here is to deal with the inner geometry of a manifold which is 
not part of a Euclidean space. Ifthe manifold were part of Euclidean space, 
it would inherit a local Euclidean geometry (distance function) from that of 
the including space, as was the case for the hypersurfaces discussed in the last 
section. However, it is not part of Euclidean space, so we must postulate 
the existence of a local distance geometry. What we do in effect is to pre- 
suppose that each tangent space possesses an inner product which is 
smooth. 

Thus we let M be an n-dimensional manifold. We suppose that an inner 
product is given in the tangent space at each point P of M. Thus if 
v and w are two tangent ‘vectors at the same point P, v-w is a real number. 
The inner product is supposed to be smooth in this sense: If v and w are 
vector fields on M, then v-w is a smooth function on M. (Precisely, at 
each point P of M, the values of the given fields v, w at P are vp, Wp, 
tangent vectors at P, and we are requiring that vp-w> be a smooth function 
of P.) 

The procedure in Section 2.5 (pp. 13-14) for finding an orthnormal basis 
may be made constructive, smooth operations at each step. We know that 
there exists on each local coordinate neighborhood on M a set of n vector 
fields, forming a basis for the tangent space at each point of the neighbor- 
hood. We convert these fields to orthonormal ones to arrive at smooth 
vector fields 

e.,°'',e, 
defined on the local coordinate neighborhood in question and satisfying 
F 


e,°e; = 0;;. 

Pretending for a moment that we are observers constrained to the manifold 
M, we set out on a somewhat symbolic voyage, hoping in doing so to motivate 
the right steps. We let P denote the moving point on M and wish to think 
of its arbitrary displacement dP as a tangent vector with differential form 
coefficients; we hopefully write 

dP =) c,e; 
with o,,°°:,6, differential one-forms on our neighborhood. (Since dP 
is in no sense an exterior derivative of anything, we distinguish this d by 
bold-face type from the usual d. The same applies to the de; below.) We 
expect to arrive at a basis o, ,:-:, 6, for the one-forms which in some sense 
is orthonormal, dual to the basis e, ,--:, e, of vectors. 

We must be guided by our experience in Euclidean space. There we 
would take.a coordinate system u', --: , u” and without hesitation write 


.[ @ 
dP = »y w'(=). 
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the vectors 


0 7) 


Out?’ Out 


being the natural frame associated with the coordinate system. But pre- 
cisely this can be done on M. The expression 


dP=} an'(=) 


is independent of the local coordinate system. Indeed, if %',---, a is 
another system, then 


7 ou’ 
di' =) aul dul, 
rh du* / 
(sa) ~ © gala): 
{0 ou du*® =f a 
vi— = —_— J] ee 
mee (=) rT a (=) 


0 0 
=). 3 hau!( 5) =) (5); 


that is, dP is the same either way. 
Having this, we express the natural frame in terms of the orthonormal one, 


Oo 
(5) = ))a4;e; 


and so 


and solve for the ¢ ik 
; rs] 


d dubayje; = )) o;e;, 
CG; — ya, du’. 


We have reached our first equation of structure for local Riemannian 

geometry : 
dP =) o,e;. 

Our next venture is to attempt an analogue to the equations for the displace- 
ments de; of the vectors of the moving frame. Here we make an essential 
departure from what we did with surfaces and hypersurfaces; we must 
restrict our attention to “‘the tangential component’”’ of de; for the simple 
reason that we are constrained to M and can see no happenings in any- 
“normal” direction. Thus we seek expressions 


de; = Y aye; 
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with one-forms w;;. We try to find such w,, so as to be consistent with these 
conditions: 


(1) de,-e, + e;:de, = 0, 
(2) d(dP) = 0. 
The reason for (1) is that d(e,-e,) = d(d,;,) = 0 and we hope to “differentiate” 
the dot product by the usual product rule. In choosing (2) we simply go 
according to the Euclidean analogue. In the usual way, (1) reduces to 
(1’) Ox + Wy; = 0. 
We explore condition (2): 
d(> CO; e;) = 0, 
> (do; — > o;0;;)e; = 9, 


so (2) is equivalent to 


(2") do; = )) 0,0;;. 
We have finally come to a well-formulated problem: given the basis 
o,,°°*, 6, of one-forms, find one-forms w,, satisfying 
(1’) Oi; + Oi, = 0 
(2') do; = )) 0 ji. 


We shall show that this problem has exactly one solution. (This completes a 
point we left open on p. 123 of the last section.) 
Since the o; form a basis we may write 


Oi; = > Vij O 


where the (unknown) functions I’; ,, are the connection coefficients, or Christoffel 
symbols. Equivalent to (1’) is 


(1) Vin + Uj = 0. 
The do, are known since the o; are known; we may write 
do, = 4) Cin j Op Cijp + yy = 9. 
We have 
do,=4) ¢:.0;0, = >, 6; @;; 
=) Cj > Dyin O, = 4 > (Cyn — Uyijlo5 oy, 
and so (2’) is equivalent to 


(2”) Die — Unig = Cijn: 


J 
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Our precise statement now is this: 
Given functions c;;, such that c;,, + Cy; = 0, the system of linear equations 
Vij, t+ Tj, = 0 
Diya — Tay = Ciik 


has a unique solution given by 
Di, = (Cy; — C jei — Ci jK)- 
For if T,,, is any solution, then we use the equations alternately to derive 
Din = “Uy = —V ey — oi 
= Diy — Cin = Uji t+ Ceiy — Cijx 


= ~—Tji t+ Gi — Cin = —T 5, — Cpe + Cees — Ci, 
hence 
20 5 5% = Cig ~ Cjei ijn 


which establishes the uniqueness of solution. It is easily verified that the 
asserted values of I’; ,, really are a solution. 
We have completed our structure equations, 


dP = ce 
de = Qe 
Q+'Q=0, 


where we have introduced the matrix notation 
ao=(0,,°°°,49,), ei: }, Q = |ia,,|| 


and already have one integrability condition, 
do = 6. 
There is no reason for believing that d(de) = 0 in any sense. We have 
d’e = d(de) = d(Me) = (dQ) e — Q(de) 


= (dQ — Je. 
We set 
© = |16,;| = dQ — 22, 


the curvature matrix which appears from the symbolic equation 


d’e = Oe 
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as representing a “‘second derivative’”—exactly how one thinks of curvature 
in elementary differential geometry. 
We derive further integrability conditions by differentiating. From 


da = aQ 
we have 
0 = d(da) = (da) Q — a(dQ) 
= (@Q) Q — a(dQ), 
hence 
a = 0. 
From 
© = dQ — O? 
we have 
d© = d(dQ) — d(Q?) 
= 0 — (dQ) Q + Q(dQ) 
—(0 +: 27)N + QO + 2), 
hence 
dO = NO — OQ, 
which comprises the Bianchi identity. 
We sum up: 

Structure equations Integrability conditions 
dP=ce do = aX 
de = Qe oO =0 

0+'Q2=0 dO = NO — OQ. 

0 =dQ — ? 


The n form o, --: 6, is the volume element of M. It is determined up to 
sign. If M is oriented, one may fix it by choosing only moving frames 
coherent to the orientation. 

The 6,; are two-forms which may be written 


OF = 4 > Rint OF 

which defines the Riemann curvature tensor. We have 

Rina t Biju = 0 

Rina + Rj = 9. 
The relation ¢© = 0, or 
» Rijn Fi O,.0; = 9, 
18 equivalent to 

Rint + Riny + Bin = 9. 
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In the special case of a hypersurface we had the symmetry condition 

Rijn = Ruri j 
as an obvious consequence of the expression for R;,,, a8 a two-round minor 
(see p. 122). Such a determinant representation is not possible in this 


general situation, but it turns out that this symmetry of the Riemann 
tensor is true anyway, an algebraic consequence of the other relations: 


(Riser — Ryrij) = Rijer + Reig + Raji) 
= (Rint Raj) — Rai 
= — Riri o? B inti 
= Riiy + (Byrn + Ryn 
= (Rig + Rijes) + Ryan 

—Ryip + By 


i 


— (Bij — Ryriy); 
and so 
2( Bs jer — Ryris) = 0 
which implies the symmetry in question. 
Those who have been through the mill in Riemannian geometry a la 
classical tensors will be anxious to see the connection. There one deals with 
a natura] frame 


0 o 
iat 


due to a local coordinate system (u',-+-, wu"). One sets 


913 = Vir Vj 
defining the positive definite symmetric matrix (metric tensor) 
= |l9; ll. 


Usually one looks instead at the corresponding definite quadratic form 


ds’? =) g,{du' du’] 
where the brackets remind us that this is ordinary, not exterior, multiplica- 
tion of differentials. This is motivated by the formula 


=|Eolz)G)] 


for the arc length of a curve u! =u'(t). 
Now one has 
dP = Y du' V;. 
We try 
dv, =) nv; 
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and then introduce Christoffel symbols by 


ni=> fA dut. 


To have 
dg;; = dv;:v,; + v,' dv, 
and 
0 = d(dP) = —) du'dy, 
requires 


Fa Liajou+ Ej} 
a ~ li} 


H= Dhow, 


Lowering indices: 


the equations become 


ee 
Lj, He] + (i, 5k] = 4 
[jf] = [j,i] 


(24 09 ix a 


with the usual solution 


Ls IH = Sak + Gal 7 oul 


and so it goes. 

This digression out of the way, we briefly consider parallel displacement. 
Let v be a tangent vector on M which is a function of one or more variables. 
We write 

v=) ce; 


where the c, are functions and have 
dv = ) de,e,; + )' ¢,0,;e; 
=) (de; + } c,@;) e;. 
A vector v moves by parallel displacement if dv = 0, i.e., 
de, + Yi c,@,; = 0. 


With this interpretation of dv one has for two compatible vector functions 
V, w that 
d(v-w) = (dv): w + v-(dw). 


Consequently if both v and w move by parallel displacement then v-w is 
constant; its differential vanishes. 
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A curve P = P(s) on M where s is the arc length is called a geodesic if 
the unit tangent dP/ds moves by parallel displacement. 


Example. We consider the upper half plane with the (Poincaré) metric 


_ [dx]? + [dy]? 
y? ) 


ds? 


Here 
dx dy 


6, >; 62 >—~—. 
y y 


Since the position vector is P = (x, y) we have 


dP = (dz, dy) =0,e, + a,e, 


so that 
e, = (y, 0), e, = (0, y). 
We have 
dxd 
do, = — = 06,0, do, = 0, 
hence Zz 
0 «6 
do, ,@0,)=(6,,0 ( ), 
( 1 2) ( 1 2) —O; 0 
~ 0 «a, 
7 Le 0 ) 


_ a 0 do,\ _ 0 1 
0©=dQ-Q = (oe, 0 =(_ 5 9)%1%- 


The only significant component of the curvature tensor is 
Riya. = 1. 


In our discussion of surfaces we wrote dw + Ko,c, =0 for the Gaussian 
curvature K. Here the right interpretation is w = @,, =0,, K = —R,> 12 
=: —1. For this reason we say that the upper half plane with the Poincaré 
metric has constant negative curvature. 

We shall show that each semicircle 


P= (a+ rcost, rsiné) O<t<n 
orthogonal to the z-axis is a geodesic. In coordinates it is given by 


x=a+rcost, y=rsint, 
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and we have for the tangent 


dP/dt = r(—sint, cost) 


= vl —sint)e, + (cost)e,] 

os : in ¢ t 

= sayy sin ye, + (cost) e,]. 
Hence for the arc length s, 
ds__ [(—sint)? + (cost)? ge 
dt sin? ¢ ~ sint 
t = dP/ds = (—sint)e, + (cost)e,. 


Along the curve we have 


@12 = 0, = —dl, 
de, = —dte,, de, = dte,, 
dt = d[(—sint)e, + (cost)e,] 


—(cost)dte, — (sint)dte, — (sint)(—dte,) + (cost)(dte,) 
= 0, 


the unit tangent moves by parallel displacement, the curve is a geodesic. 
We shall close this local study with an application of The Frobenius 
Integration Theorem of Sections 7.3 and 7.4. 
Let M be a Riemanman manifold with curvature tensor zero. Then M is 
flat : there exists a local coordinate system u',+-- , u" for which the natural frame 


0 0 


So an 


ws an orthonormal frame. 
This is proved as follows. We are assuming © = 0, ie., dQ =Q?. By 
the first application in Section 7.4, there is a matrix A of functions satisfying 


(dA)A~'=Q, 
and A is orthogonal. We define t =(t,,°°-,1,) by t= o@A. Then 
dt = d(aA) = (do) A —aodA = (0Q) A — 6 (QA) = 0. 
Kach of the one-forms 1; is closed, dt; = 0, hence exact locally, 


Tt; = du’. 
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This defines our local coordinate system (u',---,u"). On the one hand we 
have 


{0 
P=) du'|— 
d > du (=) 
and on the other, 
dP =ce=1A~'e, 


djéu' 
= A~'e. 
d/du" 


Since the frame e is orthonormal and the matrix A is orthogonal, the natural 
frame (0/du', --+ , 6/0u") is also orthonormal. 


hence 


8.4. Riemannian Geometry, Harmonic Integrals 


In this section we shall sketch the remarkable results of W. V. D. Hodge 
on the potential theory of closed Riemannian manifolds. This work pertains 
to differential forms alone so we can forget all about vector fields. In this 
spirit we had better make a fresh start and reformulate the pertinent facts 
about Riemannian manifolds which we shall need. We shall presuppose 
that the manifolds we discuss are orientable, so this will be built into the 
structure. 

Thus we have a manifold M. It is covered by a system of overlapping 


neighborhoods U,, U,,---. Oneach U, there is a basis 
omy ; eee ; 6, 
for the one-forms. Ifo,,-:-, 6, is this basis on U and G,,°°-,4, is the 
one on G, then wherever U and U intersect we must have 
G i= >, a; j oO j 
where A = |ja,,|| is a proper (determinant one) orthogonal matrix. 
The volume element o,,-°--,6, 18 an intrinsic quantity, according to 


GO, °° On = [ij |Oy °° Op = Oy °° Oy. 


Next, the star operator of Section 2.7 applies. To each p-form w corresponds 
on (n — p)-form +. Locally 


#(O,° . Oy) = Ont" . °d,. 
We recall that 
wee = (—1)PO— Pg, 


We define a new operator 6 by 


b@ = (—1)"? "ted eq. 
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The significance of the sign will appear shortly. We note that dw is a 
(p — 1)-form when w is a p-form. If w =f is a zero-form, or function, then 
6f=0. The final operator we define is the harmonic operator (generalized 
Laplacian) according to 

A=d°6+6°d. 


We henceforth restrict attention to a closed (compact) manifold M. We 
denote by t the volume element, an n-form on M which nowhere vanishes 
and which satisfies t = ¢, --: a, locally. 

We propose to turn the space of p-forms on M into an (infinite dimensional) 
inner product space. If @ and y are two p-forms then w A *y is an n-form 
and we define 


(@, n) -| WA *N. 
mM 
This is evidently linear in each variable and we have 


(o, n) = (7, @), 
a consequence of 


OAH =n A *O. 


Locally, if 
w=) ayo", 
then 
WA #0 = (> af)t 
hence 


(w, @) 2 0 
and (w, w) = 0 if and only if w = 0. 
We shall now establish the fundamental formula: 
If w is @ p-form and n a (p + 1)-form then 
(deo, n) = (a, dn). 
For we integrate over the closed manifold M the relation 


dw A *n +(—1)?o A den = d(w A *n): 


| doy nn +(-19"| orden =| dwn on) = | wA yn = 0, 
M M M aM 


(dw, n) = (—1)?~* | w A den. 
mM 
Since d «ny is an (n — p)-form, we have 


xa(d en) = #(ad ay) = (—1)?""P)d xn 
So that 
(= 1)?" tenn = (—1)P-1( 1) Pad 4) 
= (—1)"? 41 «(ad 7). 
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But 7 is a (p + 1)-form, hence 


én = (—1)"?* L)+nt lad en _ (—1)"?* Lid xn), 
and so 
(~—1)?~!d«n = *6n, 


(dw, yn) = (ert | onrndxy 
M 


-|-0 A *(6n) = (a, dn). 
M 


A form @ is called harmonic provided 
Aw = 0. 


It is clear that if the p-form q@ satisfies the two equations dw=0, dw = 0, 
then w is harmonic. The converse is also true. Indeed, if @ is any p-form, 
then 

(Aw, w) = (ddw, w) + (dda, w) 


= (6, dw) + (dm, dw). 
Now if w is harmonic, then Aw = 0, 
(6m, dw) + (dw, dw) = 0. 


But each term is nonnegative, hence each vanishes, (dw, dw) = 0, (6m, dw) 
= 0 and this implies in turn that dw = 0, dw = 0. 

The operators d, 6, act on the space of p-forms. The relation (dw, n) 
= (w, dn) may be interpreted as saying that d and 6 are adjoint to each other. 
We next see that A, which maps p-forms into p-forms, is self-adjoint, 


(Aw, n) = (a, An), 


indeed, either side is (dw, dy) + (6m, dn). Since (Aw, w) = 0 with equality 
only when Aw = 0, we are entitled to call A a positive definite (or elliptic) 
self-adjoint differential operator. 

We may now state Hodge’s main result, a deep theorem in harmonic 
analysis : 

If w is any p-form then there is a (p — 1)-form a, a (p + 1)-form B and a 
harmonic p-form y such that 


w= da + 5B + y. 


The forms da, OB, y are unique. 
The proof that a, B, and y exist is difficult. We shall only settle the 
uniqueness part. Suppose we have 


da + df + y= 0. 
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We then have d(da) = 0 and also dy = 0 since y is harmonic. Hence 


ddp = 0, 
(458, B) = 0, 
(58, 58) = 0, 

6p = 0, 
da+y=0. 


Similarly da = 0, y = 0. 
By an almost identical argument one shows that in case @ is a closed 
p-form, dw = 0, then the term 68 in the Hodge decomposition of @ is absent. 
w= da + y. 


It follows from this that if zis any p-cycle, then 


[o-fp 


that is, y has the same periods as does w. (See De Rham’s theorems, 
Section 5.9.) The result of this is that 2f w is any closed form, then there 
exists a unique harmonic form y with the same periods as those of w. 

We can also answer the following question. Given a p-form 4, when is 
there a p-form 7 such that the equation 


An=A 
is satisfied? The answer is: if and only if 
(y, 4) =0 


for every harmonic form jy. 
For suppose A = An and y is harmonic. Then 


(y, A) = (7, An) = (Ay, 2) = (0, 0) = 0. 


On the other hand, suppose / is a form satisfying (y, 4) = 0 for each harmonic 
form y. From the decomposition 


A=da+op+y 
we have, using the particular y which is part of A, 


O = (y, A) = (y, da) + (y, 58) + (, ¥) 
= (dy, a) + (dy, B) + (y ¥) 


=n (y, 7), 
hence y = 0, 
A= da + op. 
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We shall set 7 =p + v and try to solve Au = dx, Av = 6B, separatety. We 
take the first 


Ap = da. 
Decomposing «, 
a=da,+dp,+ 71, 


da —_ ad op, . 
Next, 
B, =da, + 0B, + 72, 


adop, = dé da, = (dd + 6d)(da,) — A(da,), 
da = Au with p= da. 
We find v similarly. 
Example 1. E”. We shall compute the operator A in E”. Contrary to our 
previous notation in dealing with the standard Laplacian, we shall denote the 
Laplacian by Lap, 
La y y a7 u 
Dix Us = —_——— | 
e . dx dzx' 
The result is this: if 
o=) a,dz", 
then 
Aw = —) (Lapa,) dx". 
It will suffice to establish this for the monomial 
w= Adz'---dzx?. 
We shall abbreviate the calculation by these conventions: 
(1) subscripts on A denote partial derivatives ; 
(2)a=1,2,::-,pjg=ptl,---,n; 
(3) each repeated index is summed over its range. 
We also remark that in taking the star of a monomial, the choice of sign is 


always governed by the rule 7 A *y4 = Bdz'---dz" where B>O, for 
example, 


«(da* da? *} -++ dai ---dz")= (—])sPtetats da +--+ da% +++ da? dx. 
Another point: since @ is a p-form and dw is a (p + 1)-form, 
dw = (—1)"?*"* ad xe, 


ddw = (—1)"? * 14d «da. 


8.4. RIEMANNIAN GEOMETRY, HARMONIC INTEGRALS 141 
We have 


ap = AdgPt!..- dz", 
duw = A, dx*da?t! +--+ dz", 


4N 
ad eep = (—1)"?t"t*414 da! --- dat +++ da?, 
SN 
dw = (—1)*A, dx! ---dx*---dz?, 


YN 
déw = —A,,dax! --- dx? + (~—1)**?**A, de! oo da*--- da? dx!. 
Next, 
dw = A , dei dx" vos da, 


ia 
*day = (—1))*1A ,daP*)--- da! --- dz", 
oN 
d«dq@ = (—1)?A ,,dx?** ee dat + (—1)/7 1A, dx*dx?*! see da .-- - ax", 
JN 
wd xdop = (—1)"A,,da! +++ da? +(—1)Ptetett 4 da! ---da*+++ da? dal, 


es 
ddw = —A,,dx' +++ dx? + (—1)?**A,, da’ --- dat +++ dx? dai. 
Combining these expressions, 


Aw = déw + ddw 


It 


—-[A.. + A ,;)de' epee ax? 
—(Lap A)dzx! ++ dz?. 


The minus sign seems strange in view of the relation 
(Aw, w) 2 0 


on closed manifolds. An example may clear this point. Let us take the 
zero-form sin x on the flat torus 0 < x, y, z < 2x, where numbers are identified 
if they differ by a multiple of 2x. Then 
d(sin x) = (cos x) dx 
*d(sin x) = (cos 2) dy dz 
dxd(sinx) = —(sinz) dz dy dz, 
xd xd(sin x) = — (sin 2), 
A(sinz) = dd(sin x) = +sina, 


2n 
sin? «da = (22)*2 > 0. 


(A(sin x), sin x) = eny* | 


0 
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We have used the fact that d(sinz) is a one-form, n = 3, p = I, 
d[d(sin x)] = (—1)"?*"*!(«de)d(sin x) = — «d «d(sin 2). 


Example 2. S?. If f is a function on E*, we have the spherical coordinate 
form of the Laplacian (Section 4.4, p. 40), 


rnr= rg [E rane) +5 (mo) +8 (2,2) 


Suppose that Lap f = 0 and that f is of the form 
F(r, b, 8) = r'g(, 8). 


Then g is called a spherical harmonic and must satisfy 


n(n + Ising + = [(sindag]+ (242) =o. 


ap ind 
The function g may be considered as a function, or zero-form, on S$”, the unit 
sphere. There we have 


ade. Gy=singad: 
_ Jo 

dg = 940, + (4) 02; 

on ee 


Sie (* 5) dd + (sin )g, 40, 


dxdg = = (2 =F 6 


G E ( ae )+; += = (‘sin dy) O12, 


~ sing 06 \sing 
eaedy = 5 | (ara) + ag | 


Ag = (dé + éd)g = ddg = —d dq, 
so the condition for g to be a spherical harmonic is 


Ag — n(n + l)g = 0. 








) ded + = (sin bug} a 











But 


In other words, the spherical harmonics are eigenfunctions for the generalized 
Laplacian on S*. Many of the usual facts about spherical harmonics 
follow from our calculations. We take one example, the orthogonality 
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relation: If g and h are spherical harmonics of distinct degrees m and n, 
respectively, then (g, h) = 0. 
For Ag = m(m + 1)g, Ah = n(n + 1)h, hence 


] 
G.hy=— arene 9; MS may 
a + nn +t) 
(g, h) = 0. 


8.5. Affine Connection 


We shall approach the problem of affine connection this way. We seek 
the weakest structure with which we can endow a manifold so that parallel 
displacement, of vectors along curves is possible. Considerably less than a 
Riemannian structure is required. 

Let M be a manifold. An affine frame (or simply frame) on a neighborhood 
U of M consists of » vector fields e,,---,e, on U which are linearly 
independent at each point of U. Thus at each point P of U the vectors 
(e,)p,**:, (@,)p furnish a basis of the tangent space T at P. There is a 
dual basis g', --- , o” of one-forms on U so we may write 


dP= y a'e; 


as we did for the Riemannian case in Section 8.3. 

We now want to associate with each vector field v on M a vector field dv 
with one-form coefficients. We must be able to do this for the vector fields 
e, of the basis, so we require 


de,=) w/e,, 


where the w,’ are one-forms on the neighborhood U. There are certain 
consistency conditions which guarantee that the computation of dv will be 
independent of any frame. 

Locally we may describe an affine connection as sono We are given 
U, the affine frame e, ,:--, e, pane the dual basis o',-+-, 0" of one-forms. 
An affine connection ponisiats a n* one-forms w; iar ‘6 no constraints 
whatever. 

We shall develop some of the local geometry of an affine connection before 
attacking the crucial problem of finding proper consistency conditions which 
make the definition of an affine connection over a whole manifold possible. 

By introducing matrix notation: 


ec; M11 °° * Wy, 
a= ; o =(a',--+, 0"), Q= : ; 
e,, Ont °° *° Dan 
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we may write our basic structure equations: 
dP = cee, de = Ne. 


We shall quickly point out the relation to the customary tensor formulation 
of an affine connection. First we expand each w,’ in the basis o*: 


o/ = » rio, 
defining the connection coefficients [/,. In the usual tensor formulation, 
the frame e,,---,e, stems from local coordinates: 
7) 
aut , 


and correspondingly, g'= du‘. The [/,=T/,(u',---,u") are then n4 


arbitrary functions assigned on U. 
We derive relations by differentiating the structure equations several 
times. First 


d*P = (da)e — ode = (do — aQ)e = te. 
Here 
t=(ti,---, 1") =de —0Q. 


The two-forms +! are the torsion forms. We may write 
td=ty Ti,ai not, T y+ T'., = 9, 
defining the torsion coefficients T',,. Next, 


d7e = (dQ) e — Ode = (dQ — N”)e = Oe. 
Here 
© = dQ — Q? = |/6/ I 


is the matrix of curvature forms 0,/._ The curvature tensor R,/,, is obtained 


from ; 
0, =4 > Rea o Aa, Ry + Re = 0. 


Integrability conditions are obtained by applying the exterior derivative to 
the equations t = do — 6Q and © = dQ — Q?: 


dt = (—de)Q + «dQ 
= —(+60)0+46(0 + 97), 


dt =aO — 7Q, 
and 
dO = (-dQ)N + QdQ 


= —(0 + 27)N + Q(O + N?), 
dO = D0 — ON. 
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If v is a vector field, then 
v=)fie,=Fe, F=(f',---,f"), 
where the f' are scalars; we have 
dv = (dF)e + F(de) = (dF + FQ)e. 


Suppose the vector field v is defined over a submanifold. It is said to move 
by parallel displacement if dv = 0, ie., 


dF + FO =0. 


If P = P(t) is a smooth curve on U defined over an interval tp) St St, of 
the ¢ axis, if vp) is a tangent vector at P, = P(t ), then there exists a unique 
assignment of a tangent vector v(t) at P(t) for each value of ¢ such that 
v(t,) = Vp and v(t) moves by parallel displacement. 

For we write v=) f‘e;. Along the curve the conditions for parallel 
displacement become 


+ ert (S2) = 0, 


a first order linear system which taken with the initial data determines the 


f' uniquely. 
Now we tackle the global situation. We have a manifold M in front of 
us and we must consider each conceivable moving affine frame e,,°:--, e, 


together with its neighborhood of definition U. With each one of these e 
we have an affine connection, ie., an » xX » matrix Q of one-forms on U. 
We want these to “fit together’’ whenever two such neighborhoods overlap. 

Thus let (U, e, 6, Q) be one such system and (0, é, , Q) another, where 
we assume the neighborhoods U and U overlap. The basic thing we require 
is that for any vector field v defined on the intersection of U and U, the 
computation for dv in either connection must yield the same result. 

On the overlap, where we shal] always operate in what follows, 


é= Ae 
where A = |la,/|| is a nonsingular matrix of functions. Since 


dP = g = ce, 


Ql 


we have 
GA 


@ 
It 


ee, 
cA}, 


Ql 
I 


Now 
dé = d(Ae) = (dA)e + A(de) 
= (dA + AQ)e 
= (dA + AQ)A™~'6. 
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But also dé = 9@, so we have the transformation law for Q: 
O = AQA-! 4 (dA)A7?. 


This is quite different from the transformation law for o which is forced 
by the very definition of manifold and frame. It tells us how the various 
matrices {2 we are associating with the various moving frames e must be 
related if we are to define an affine connection on M as a whole. 

From these formulas one can derive the transformation laws for t and Q: 


@=1A7}, 
© = AQA™!. 
Ifv = Fé = Feisa vector field, then F = FA~}, 
(dF + FO)S 
= (dF A~!— FA~'dAA~')Ae+(FA~1)(AQA-!4dAA~') Ae 
= (dF + FQO)e, 


so that dv is the same, either way it is computed. (We have used the rule 
d(A~') = —A~'(dA)A~", as we shall several times.) 

We shall now have a second look, only in our present context, at the 
considerations of Section 4.3. This will provide us with another way of 
looking at affine connection. What we shall do fits into a general pattern: 
quantities subject to a transformation law become absolute invariants when 
considered on a suitably extended space. 

We begin with a manifold M of dimensions. We form a new manifold 
F of dimension (n + n*). This frame manifold consists of all frames at all 
points of M. Precisely, at each point P of M consider all possible bases 
e,,°°:, e, of the tangent space Tp at P, and do this for all P. 

We obtain coordinates on F this way. First let U be a local coordinate 
neighborhood on M with coordinates u',---,u". Pick a moving frame 
e,,:°:,e,on F. Thus at each point P of U, (e,)p,---, (e,)p is one basis 
of the tangent space Tp at P. The most general basis of Tp stems from this 
given one by applying an arbitrary nonsingular transformation. It is 
f,,°-:. f, where 

f,= >, b,/-(e;)p 
with ||b,j| an arbitrary n x 1 nonsingular matrix. It is clear from this 
that the (n + 7) independent variables 


1 n j 
De 


serve as a coordinate system for the neighborhood in F consisting of all 
frames at all points P of U. 
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With the moving frame e,,-:-, e, on U goes the dual basis o!,---, o” 
of one-forms on U. The forms 6’, --- , é" defined by 
é=ocB"', 
where 


o=(o1,-:-,o"), &=(6,---,6), B= ]b/I), 


are one-forms on the part of F lying over U. More precisely, the values of 
the forms ¢!,---+ at the point f of F given by 


f,= >, b/(e;)p 


Gl¢ =(o|p)B'. 


are 


Now suppose that O is a second coordinate neighborhood, & a moving 
frame on O, etc.. Suppose that P lies both in U and in U0. We have 


é= Ae 
where A is a matrix of functions, 


¢=o0A"}, 


all worked out above. If the point f of F has coordinates B with respect to e 
and B with respect to é, then 


Be = Bé = BAe, 
B=BaA. 
Thus 
é=oB '= (@A)(BA)~? =¢@ B-', 
This implies that the one-forms G', +--+, 6" are defined on all of F and are 


completely independent of the particular local coordinate neighborhoods and 
moving frames used in their definitions. 

This is of first importance. We began with an n-manifold M. We 
constructed over it a new manifold F. On this new manifold we auto- 
matically have, free of charge, the n (linearly independent) one-forms 
g}, a 

Now suppose an affine connection is given on M. Thus to the neighbor- 
hood U with a definite affine frame e is given a matrix Q = ||w,’|| of one- 
forms. The matrix 0 corresponding to the frame & on UO is related to Q 
(on the overlap of U and Q) by a certain transformation law. This trans- 
formation law means nothing more nor less than that the n? one-forms / 
defined by 

O = ||6,/|| = BQB-* + (dB)B™', 


apparently defined only on the part of F lying over U, are defined on all of 
F, are completely independent of e and its U. 
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This statement may be verified by a calculation. According to the nota- 
tion above, one must prove the formula 


BOQB™' + dB) B71 = BQB- + (dB) B-!. 
We leave the details to the reader, but point out that the result can be 
motivated by a symbolic calculation: 
f= Be, 
df = (dB)e + Bide) = (dB + BO)e 
= (dB + BQ) B= 'f, 
similarly 
df= (dB + BO) B-'f, 
hence 
(dB + BQ)B-! = (dB + BO) B-?!. 
From the one-forms 6',---, @,/ on F, one constructs two-forms 7, 6 
according to 
t= (t', +++ 0") = de — 6X, 


© = 67] = d& — G?. 


These are the general torsion and curvature forms respectively. 


8.6. Problems 
1. Let £ be a closed convex surface with constant mean curvature H. 


Prove that 2 is a sphere. 
2. A surface X is given in the Monge form z = f(z, y), defined for all z, y. 
We suppose this surface is convex from below. Show that this means that 


r 8 
(: 
is positive semidefinite. 
3. (Continuation.) Show that the mapping 


(x, y) —> (w@ +p, 4+ 9) 


increases distance and hence is one-one. 

4. A point of a hypersurface is an umbilic if the transformation A at that 
point has all of its characteristic roots (principal curvatures) equal. Let M 
be a hypersurface all of whose points are umbilics. Prove that M is a 


portion of a hyperplane or sphere. 
5. Suppose on a Riemannian manifold M there is a scalar K such that 


6,,= —Ko;A G;- 


Prove that K is a constant. 
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6. Let M be a manifold with an affine connection given. Show that 
the two-form « defined locally by 


a=) 0,' = trace O 


is actually defined on all of M, independent of local frames. Show also that 
da =0. Itis even true that there exists a one-form J on M such that « = dA, 
but this is difficult to establish. 

7. (Continuation.) Prove 


dQ’ = QO" — 0" 
and that 
d[trace (O")] = 0. 


8. Let M be a manifold with affine connection. Given an affine frame 
e,,''',@, with corresponding connection coefficients F and torsion co- 
efficients 7’, define a new connection by specifying the new connection 
coefficients I’*: 

mf, =e, + 47 u- 
Show that this indeed defines a connection on all of M and that this con- 
nection is symmetric (no torsion). Investigate the meaning of “symmetric’’ 
for a local coordinate frame. 

9. Consider the flat torus T” which consists of all points (z,,°°-, 2,) 
where each 2, is taken modulo one. That is, (z,,---,2,) = (Y,,°°*s Yq) 
if z,— y, = integer (1=1,---,). The flat metric is given by the ortho- 
normal basis 

G, =dx,,°*', o, = dz, 


of one-forms. (In older notation, [ds]* = [dx,]* + --- +(dz,]*.) Find all 
harmonic differentials of all degrees. 


IX 


Applications to 
Group Theory 


9.1. Lie Groups 


A ne group consists of a smooth manifold G which has a group structure 
(x, y) —> xy. 
We suppose that this group operation, which may be considered as a mapping 
G x G — G, 


is smooth and also that the map z —> x~! on G — G is smooth. 
With each element x in G, there is associated a transformation L, of 
G, called left translation: 


L(y) = xy. 
A differential p-form q@ is called left invariant provided 


L.*o=@ 
for all x in G. 

Let e denote the unit element of G. The left translation L,~! = L,-: 
sendsxtoe. Ifw is left invariant, @ = LY-,w is completely determined at x 
by its value wy ate. If wo is any given p-form at e, then a left invariant 
form w is defined by 

oO, = Dyes Wo- 


These remarks serve to determine the existence of left invariant forms. 
Let us begin with one-forms. We let » be the dimension of G. Since the 
space of one-forms at e is an n-dimensional linear space, there are exactly n 
linearly independent left invariant one-forms on G. Let 


oi,+++,o" 


be such a system. Any other left invariant one-form is a linear combination 
of these with constant coefficients. 
More generally, if m is any left invariant p-form on G, 


w=) cya", 
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where the c, are constants and go? =g"---o*?, Any p-form @ can be 
expanded in this way and the coefficients c,, will, in general, be scalars on G. 
Supposing w left invariant forces each of these scalars to be left invariant. 
This means that each c, takes the same value at each point of G, hence is 
constant. 

Next, if w is left invariant, so is dw, since 


L,*(dw) = d(L,* w) = dw. 
It follows that there are constants of structure c' , such that 
da'=}4) c,0!o*, cy +ci,, =0. 
Substituting this into the relations d(do') = 0 eventually yields 
Y (Ch cts HC jC sy + Cy gCiy,) = 0. 
Particularly important is the n-form 
gi.+sg" 


which defines a left invariant volume element on G. It is clear from this 
that G is orientable. 
The right translation R, associated to a group element z is the mapping 


y—> Ry =y% 
on Gto G. From the associative law 


a(yz) = (xy)z 
we deduce that 
Lok, = R,.L,, 
hence 
Re 6 b= gk. 


Suppose w is a left invariant p-form. Then for each zx and z, 
L,*(R,* @) = R*(L,* w) = R.* @, 


hence R,* w is also a left invariant p-form. 


9.2. Examples of Lie Groups 


Example 1. 2=1. We shall determine the local structure of all one- 
dimensional groups. Let ¢ be a parameter on G, chosen so that f = 0 is the 
identity e. Let o be a nontrivial left invariant one-form; locally, 


o = f(t) dt, never zero. 


We integrate o to get a new parameter for G, 


| fit) dt. 
0 
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Thus we may assume we have started with a parameterization of a neighbor- 
hood of e by a single variable ¢ such that 


o = dt 
is a left invariant form. 

We next express the group product analytically. The product of the 
point with coordinate s with that of coordinate ¢ will have coordinate u 
given by 

u = p(s, t) 
with 
p(s,0)=s, pl0,t)=t 
according to xe = z,ey = y. In coordinates, 
L,: t—u= p(s, t). 


The left invariance of o, L,*o = o, means 


Op 
dt = — dt, 
t AI dt 
hence 
0 
T=h  pleat)=t+ oe) 
Setting ¢ = 0: 
p(s, 9) = (s), 
and so 
p(s, t) =s +t. 


It follows that the group operation corresponds to nothing more than ordinary 
addition of coordinates. As a corollary, G is abelian (commutative). 


Example 2. G is a group for which the constants of structure c’,, all vanish. 
Thus 


do! = = do" = 0 
In a small neighborhood of e, 
oi = du’, 
taken so that (0,0,:-:,0)«+e. The product of points with coordinates 


u, V, respectively, is a point with coordinates w given by 


w' = p'(u, V) 
with 
pi(u,0)=u'i, — pi(0, v) = 0". 


The invariance of co‘ under the left translation 


v—w, 3 wi = p'(u, v) 
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is expressed analytically by 


dvi=) e (u, v) dv! 


jal 


oP, Vv) = 0; 


which implies 


pi(u, v) = vo! + '(u) 

Setting v = 0 yields 
ul = '(u), 
hence 
pi(u, v) = ul +o! 

If P(u) denotes the point with coordinates u, this says 

P(u): P(v) = P(u + v) 
so that locally the group looks like a neighborhood of 0 in E”. 

Corollary. G is abelian. 


9.3. Matrix Groups 


Now we shall consider a group G which is a smooth subgroup of the group 
GL(m) of mx m nonsingular matrices. (The notation stems from the 
common name, general linear group.) 

Suppose u!,--+, u" is a coordinate systemonG in some neighborhood of J, 
theidentity matrix, and that X = X(u',---, wu") is a typical point in this 
neighborhood. The matrix dX of one-forms certainly contains n linearly 
independent one-forms because the n-dimensional group G is smoothly 
imbedded in GL(m). Consequently the matrix 


Q= X7'dXx 
of one-forms contains n linearly independent ones. But each element of Q 
1s left invariant. For if A is any fixed element of G, the left translation 
by A is given by 
xX— AX, 
while 
(AX)~'d(AX) = (X~'A7')(AdX) = X7' aX. 
This allows us to make explicit calculations in many important groups. 
Next we note an important geometric interpretation of Q. We interpret 
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each element X of G as a linear transformation on the space E” of row 
vectors Vv = (v,,°°*,v,). Thus 


v— w =vik. 


We ask, how does dw grow out of w under the group action? Here v is 
fixed and X varies over G. We have 


dw = vdX = (wX~')dX, 
dw = wQ). 


This means Q can be interpreted as an “infinitesimal group element.” 
(Cf. Section 4.2.) 

One final remark. The constants of structure can often be explicitly 
obtained from these considerations: 


Q=X-'dX 
dX = XQ, 
0 = d(dX) =dXQ+ XdO 


= (XQ)Q + XdQ. 
Hence, 
dQ +Q?7 = 0. 


9.4. Examples of Matrix Groups 


o-((5 sees 


the proper affine group on the line. One easily sees the isomorphism between 
G and the transformation group 


Example 3. 


t{— > zt + y. 


G3) 104) 
0 1 r\0 2 
] ( ~ ’) (C °) ] e dy 
Q=- =— ; 
x\O0 -2/\90 O x\ 0 4) 


Hence o' = dz/x, o* = dy/x are left invariant. The left invariant volume 
element is 


Here 


1 2 


dx d 
CG AG eo 


¥ 1 
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Since do! = 0, do? = —dxdy/x* = —o' a o’, the only significant constant 
of structure is 

C2 = -¢3, = —1. 
If we seek right invariant forms, we find them in 


@xyx-t = (9 ON Y=G = 
a\O O/\O a] 2x\0 0 


so a basis is 
dz = —ydx+ady 
a 
The right invariant volume element is 
; ,  daedy 


a AC = ’ 
x 





very different from the left invariant one. Also 


_ dady ‘ 


da? =a a a’. 





We shall compute the effect on o? of the right translation R,, where 


a b 
A = (; > 
We have 


R,(X) = XA = (; ne ) = fe ae | 


d(b bdx ldy ob 1 
(bx + y) fe es og eg 
ax an ax a 


Example 4. The step transformation group of all matrices 


x = (5 ‘) x> 0. 


0 x)’ 
1 /x —y\(dx dy ] (edz xdy— ) 
=X '4X = — = — : 
a : =(¢ ‘)( 0 L) a 0 xdx 
We may take 
d dy ~ yd 
oi = wash = d(In2), o = slat EL La = a(2). 
x x x 
The choice of new coordinates 
u=Inz, v= id 
x 


follows the procedure in Example 2. 


156 IX. APPLICATIONS TO GROUP THEORY 
Thus 
e" ve" 1 »v 
X= == @” p 
(0%) = (0 3) 


afl 6 
A=e'() :) 


1 b+4 
— path 
AX =e (, 1 ). 


If 


is another, then 


The invariant volume element is dxdy/x* = dudv. 


Example 5. G=GL(n), the general linear group of all nonsingular n x 2 
matrices. The general element is X = ||z,/|| where the z/ are independent 
variables (subject to the inequality A = det (z,/) 4 0). Set 

l 


Y = cof X, X~1=—-/Y. 
co n 


We have 
Q = X-*dX = llos*l, 
where 


1 


The n? left invariant forms o;* are necessarily linearly independent. We 


compute the volume element 
t= |[o;' 
i,k 


in two steps: 


J ‘ j; 
y* as o,* Atta o, = qe (291" a7") 2 ae 2 Oo yn! de ;*) 


1 
= a det (y;") dx," sos dak, 


From XY = AI we have det (X) det (Y) = A”, hence 
det (Y) = A"™}, 


1 
k k k 
y = —dr --- dex 5 
A? 


] A 
t= — (dx,* --+ dx, *). 
=1 


A x 
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It is clear from this that the right invariant volume element will also be 1, 
which is an unusual feature of this group since it is highly non-abelian. 


Example 6. G = SL(n), the special linear or unimodular group of all n x n 
matrices of determinant one. The special feature we shall note is 


trace (Q) = 0. 


This follows from a general formula for a matrix function X of any number 
of variables. Set A = det(X). The formula is 


d 
= = trace (X~1dX). 


This is proved as follows. Denote by 


xy? 


the jth column of X and consider 
A = det (c}, ---, ¢”) 


as a function of the columns. Then 


dA 


y A(e!, cee c/-! de/, citl cee c”) 
j=1 


y ) y;' dx, = trace[(cof X)dX] 
j=1 1 


= trace[AX~!dX]} = Atrace(X~1dX). 


For G =SL(n) we have A= 1, dA = 0, hence trace(Q) =0. For n = 2 we 
have 


u v 


ees _{ v —y\{da dy 
ah ee (_* a Po 
_ ( vdx — ydu ie Ad _ (°; *.) 


—udx+adu —-udy+2xdv a oat} 


x=(* I A=2v—yu=1, 


Differentiating A = 1 yields 
adv+vdx—ydu—udy=0, o'+o0*=0. 
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For the left invariant volume element we may take 
t=o'aorac’ = drdulvdy — ydv) 
= vdxdudy — ydxdudv. 


Example 7. G=0*(n), the proper orthogonal group of all n x matrices 
X for which 
2. Gena, Gea det (X) = +1. 


Here the superscript ¢ denotes transpose. The essential feature about 2 is 
that it is a skew-symmetric matrix, 


Q+'Q=0. 


Because the group G has dimension n(n — 1)/2, it follows that the elements 
above the main diagonal in Q form a basis for left invariant one-forms and 
their product is the left invariant volume element. We establish this 
property of 2 as follows: 


X’X=I1, (dX)'*X + X(dX)=0, 
X-1dX + (dX)'X-! =0, 
X-!dX +(X-'dX)=0, 02+'Q=0. 


For n = 2, 
cos@ sin@ 
ais Bee: ear 
_ y-17y _. {e080 —sin@\ /—sin®@ cos0\,, { 9 1 
ae ax = (tno | ow ting)? = (3 0) 


For n>2 the calculation becomes complicated and hinges on explicit 
parametrizations of G. The cases of even and odd 7 are rather different. 


9.5. Bi-invariant Forms 


We take a Lie group G with identity element e. Because G is a manifold, 
there is a coordinate neighborhood U of e with a local coordinate system such 


that the coordinates of e are (0,0,---,0). Suppose x and y are very near 
toe andin U. Then z=xzy is in U. [If the coordinates of these three 
points are (z',---, 2"), (y',---, y"), (z’, ++, 2”), respectively, we may write 


zi ziigl,--- at y,--+, yy’). 
Since xe = x and ey = y, we have 
z'(a!, +++, a, 0,°°+, 0) = ai, 
20,°°: Oy, x) =yi, 
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and because of these facts, 
z' = x! + y' + (higher order terms in the 2 and y*). 
In particular, if y = x~1, then z = e, and 
0 =z! + y' + (higher order terms). 


We apply these simple remarks as follows. Let yw denote the mapping 
W(x) =r a 
yw: GG. 
If we write 
= W(x), 
yi =y"(at, ++, 2"), 
then by the relation just discussed, 
dy' : 
ao || macs 
Ox! jo : 
This means that 
y' = —z' + (higher order terms in the 2”). 
We may also express this another way: 
w: (xt,-++, 2") —> —(z!, +++, 2") + (higher order terms). 


Since y(e) =e 1 =e, the induced mapping w* takes each differential form 
at e to another form ate. Evidently we have 


F, y*(dx') = —dz' _— (ate), 
; w*(da' --- dx?) = (—1)?(dz' --- dx?) (at e). 
Thus if w, 1s any p-form at e, then 
w*(w.) = (—1)?a,. 
For each y in G, the right translation R, was defined by 
R, (x) = xy. 
A form o is called right invariant if 
Ry*¥a@=@ 


for each yin G. Our first application of the map y is the following: 
A form w is left invariant if and only uf w* w ts right invariant. 
For 


W(R,(x)) = (ay) = (ay)! = yo ta! 
= L,-.(x~') = L,-1(w(z)), 
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hence 
W ° kh, — Ly-1 o na 
R,* o W* =p 0 
If w is left invariant, then for each y in G 
B,AY*o) = Y(L,-*@) =" 0 


hence y* w is right invariant. Similarly, if @ is right invariant then y* @ is 
left invariant. Since (v~')"'=2, wow=1 and so W*(~*w) =o. It 
follows that if y* qw is right invariant then q@ is left invariant. 

Next we shall see that using right invariant forms instead of left invariant 
ones does not give additional constants of structure. For let a',---, 0" be a 
basis of the left invariant one-forms. Then the corresponding constants of 
structure are read from the equations 


doi =4 TS cgol aot. 


The forms t' = y*oa',--+-, t” = y*o" are now a basis for the right invariant 
one-forms. But y* applied to our equation yields, since the c’s are constants, 


W*(do') = ty ip w*(a! A o*), 
d(y*a') = i>, cf (y* a) A (y* a*), 


dti=4)e,ce a tt. 


Now we pass on to the study of bi-invariant forms, i.e., forms which are 
both left and right invariant. We derive one important result. 
Let w be a bt-invariant p-form. Then 


dw = 0. 
For W*o is left invariant since w is right invariant. We know from our 
calculation on the previous page that at the point e, 
w*(w,) = (—1)’a,. 


But w and w*(w) are both left invariant, hence what is true at e is true 
everywhere, 


W*(@) = (—1)?o. 
On the other hand, dw is a (p + 1)-form, also bi-invariant, so the same 


conclusion applies: 
w*(dw) = (—1)?*' do. 


But, 
w* (dw) = d(w* w) = d[(—1)’a] = (-1)? da. 


From these equations follows dw = 0. 
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We apply this to the case in which G is a commutative group. Then the 
left and right translations are the same thing so that each left invariant form 
is bi-invariant. In particular if ¢',---,o" is a basis of left invariant one- 
forms, each do‘ = 0; the constants of structure all vanish. In Example 2 
of Section 9.2 we showed that any group with vanishing structure constants 
has the local structure of Euclidean space (and incidentally is commutative). 

Here is one more result on bi-invariant forms which goes in a different 
direction. 

Let G be an n-dimensional closed (compact) Ine group and let w be a left 
invariant n-form on G. Then w 18 bi-invariant. 

For each z in G, R,* @ is also left invariant. Assuming w # 0, we have 
R,.* w = f(x)w, where f(z) is a real number, since the space of left invariant 
n-forms has dimension one. Because R,*. R,* = R,,*, we have f(xy) 
= f(z)f(y). (Real numbers commute!) Now f(x) never vanishes since 
1 = f(e) = f(x)f(z~'). Thus f maps G into the reals R with 0 removed. 
Since G is compact the image of G under f is a bounded interval in R, 
bounded away from zero. It is also a subgroup of the multiplicative group 
of positive reals since f preserves multiplication. [Positive because the 
image f(G) of the manifold is an interval and it contains 1 = f(e).] If f(G) 
contains any real number a # 1, then a” —> oo or a” —> 0, both impossible, 
since a" must remain in the interval {(G) which is closed under multiplication. 
Hence f(G) consists of 1 alone, f(z) =1 for each x in G, R,*w =, o@ is 
bi-invariant as asserted.(We have used the fact that G is a manifold, hence 
connected, to conclude that f(G) is an interval.) 


9.6. Problems 

1. Let C* denote the multiplicative group of nonzero complex numbers. 
Find the invariant volume (area) element. 

2. Consider the 4-dimensional group of all matrices 


(0 4) 


where z and w are complex numbers, z 4 0. Determine constants of struc- 
ture. Show that the left invariant volume element is 


—(dz a di a dw a dw)/|z|*. 


Here dz = dx — idy if dz = dx + idy. 

3. Discuss other groups of complex matrices analogous to the examples 
of Section 9,3. For example, discuss the relation between unitary matrices 
and skew-hermitian ones. 

4. Extend the coordinate considerations of Section 9.5 by showing that 


zi=aity' +) a',a/y* + (terms of order three and higher). 
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Show also that 
(a! jk a’, jp 
are the constants of structure for a suitable basis of left invariant one-forms. 
Compare xy and yx. 

5. We know that each left invariant p-form can be expressed in terms of 
left invariant one-forms. Does a corresponding result hold for bi-invariant 
forms? 

6. Letc',, be constants of structure of a group and set 


i ak 
Ii = > CO in - 
i,k 
Show that g,, is a symmetric tensor. Now set 
l 
Ci = » CuI: 


Show that c,,, is a skew-symmetric tensor. 
7. Let z be a p-cycle on G. Show that for each closed p-form w and 


each g in G, 
[o- [tro 
z z 


(We must assume that G is connected, i.e., consists of one piece only.) 


Applications to 
Physics 


10.1. Phase and State Space 


We propose to study a holonomic mechanical system with a finite number 
of degrees of freedom, avoiding collision phenomena. In this section we 
formulate the geometry of such a system. 

The position space is simply an n-dimensional manifold M. 

We next define the phase space attached to M. This is the space of all 
covariant vectors at all points of M. To make this precise, we consider a 
coordinate patch U on M with local coordinates 


q yoy gl 
At a point P of U, a covariant vector is simply a one-form at P, hence is 
given by its components 
Pis***s Dns p; real 


(where the one-form itselfis }' p,dq'). 
If 

q' yo g" 
is another local coordinate system valid at P, then the components of the 
same covariant vector with respect to the q' are 
dq! 
oq’ 
The totality of all such covariant vectors at all points of M constitutes the 
(2n)-dimensional phase space P. To each coordinate neighborhood U on M 


p=) D; 


with local coordinates g', -:-,g" corresponds the coordinate neighborhood 
U 4 E" with local coordinates 
QO Pir ts Dae 
It follows that the one-form 
; a= > p,dq 
is a one-form on P, entirely independent of local coordinates. We have 
da =} dp,dq' 
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so that the phase density (see Section 2.3) 


dp, --- dp,dq' «+= dq” 


is a 2n-form on P, never zero, defined by 


+ (da)” = (n!)(dp, a dp, dq' - -dq"), 
and serves us as a volume element on P. 
We shall derive some useful] relations from the transformation of co- 
ordinates 
q' = g(q', ore »g') 
: éq/ (1=1,--++,n) 
Pi = y Pj ag 
valid on the overlap of local coordinate neighborhoods U and O. 
We set 
P, 
q= (q’, tee g”) and p= 
Py 
and define q and p similarly. Then @=4(q) implies 


dq=dqA, A= bet = A(q). 
Since a= dq p= dq p, we have 
p= Ap, that is, p= Ap, 
where A= A~ is also a Jacobean matrix. From 
dp=dAp + Adp and dq=dqA 

we deduce first that (0p,/ap,;) = A = (dq//0q"). _ 

To continue, we note two relations. From A=A~ we have 

dA=-A-'dA A“*=—-AdA A. 

From dq@=dqA we have 0=dq dA, which easily implies 





da', _ da’, 
dq! 7 agt 
Now 
dp=dAp+Adp= — AdAA p+ Adp 
= —AdAp+ Adp. 
Also 


dp =d( Ap) =dAp + Adp. ‘ 


Thus 0p,/dq* is the coefficient of dg* in the 1-th row of — A dA p and dp, /agi 
is the coefficient of dg' in the k-th row of dA p. Now 
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Z 
iaee (S 0 sar) 
dq* 


and 





Therefore the coefficient of dg* in the i-th row of — Ad Ap and the coeffi- 
cient of dq' in the k-th row of dAp are respectively 


OD. Oa. a al. 
ia~ Sas, and oP x = » eh. 
dq dq 0q; dq‘ 


We conclude that 0f,/dq* = — dp*/dg', so we have proved 


OR; _ OP 
dq" og! 
Op; = dg! 
6p, ag’ 


Finally, the state space is the product 


S=P\ E! 
a (2n +1) dimensional space. We think of E’ as the time axis. Local 
coordinates for S are 


gy-c'y g's Dig 8" sDes t. 


10.2. Hamiltonian Systems 


We wish to consider a dynamical system in Hamiltonian form. We begin 
by tracing the evolution of this from Lagrange’s equations of motion, which 
in Euclidean coordinates reduce to Newton’s law of motion. We deal only 
with conservative holonomic systems. 

The treatment first of all is local. We deal with a coordinate patch in 
q',--+,q@" space. For each instant of time, there is a point (position) 


(q'(t), +--+, 9°) 


which represents the trajectory of the system. As is customary, we set 
u=du/dt. The kinetic energy is a function 


(iC rr a ce i) 
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which is supposed to be a positive definite quadratic form in the variables g!, 
The potential energy is a function 


V=V(q',--+, 9", #) 
and the Lagrangian function, or kinetic potential, is 
L=T-YV. 

The differential equations of motion are then 

d (oL OL 

a Peas tion 5 p= 1,+++ in). 

du (=) aq" : . 
For the first term we have 

d (dL a7 aL . OL. . 

Silat = re + Pf 59 + “3 45 qf 

dt \dq agi dt = dg’ dg aq! ag" 
so that the Lagrange equations are a system of n second order ordinary 
equations for the unknowns q!, --:,¢". We now convert these to a system 


of 2n first order equations in 2n unknowns. 
We introduce the generalized momentum components 


Prs**' > DPr 
by 

_ aL of 

tT agt ag! 


Because the quadratic form 7’ is definite, the transformation of variables 


Ch A ee i) one ( ee Ane cc 


is a smooth one both ways. 

To reach the Hamilton form, we shall follow tradition and use a rather 
confusing notation. The matter was better expressed in Section 3.5. 

The function 7' is always considered as a function of the 2” variables 


Ce Ak a 
The function V which involves the q' (and #) alone may be considered as a 
function on the space of variables g', ---,q", g’,-::, @", € or on the space of 
variables q', +++, 9", P1,°°*s Dao t- 


We introduce the Hamiltonian 
H=H(q',-++. 0 Pis** Par th => vig — L, 
always considered as a function on the space of variables 


Qt Qi Pis t's Darl 
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Since 7' is homogeneous quadratic in the g' we have 


OT 
2T = {== y*, 
» dg’ > pi 


hence H = 27 —- L=27T —(T —D), 


H=T+YJ, 


and so H represents total energy. 


From 


we have 


But 


Subtracting, 


From this, 


Thus 


2T = > pid 


2dT => p,dg' + ¥ q'dp;. 
aT = (San! + Fae 


=D Soa +E pai! 
iT = -L = a! + adp, 


ss 


av | 
dH =dT +dV = z(-5 +S a)aa! + Dade, 


OL 
= “Lait + Lee 


dH OL sd (OL . 
égi ag’ St (i) = ae 
oH 
Op; ¢ 


This gives us the equations of motion in Hamilton, or canonical, form: 


ae oH 
Op; 
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We next check what functions H qualify as Hamiltonians. We write 


T=} » a,,(q)q' ¢ 


where ||@,,(q)|| is a symmetric positive definite matrix function of the position 
variables g. It is convenient to set 


lo(q) il = las(@)l-* 
also symmetric, positive definite. Then 


oT 
Pi = ag = » A 


which we invert to 
q' = 1 b* p,. 
This gives us 
T=}) a,b bpp, 

=} > bpp, . 

From this, 
H(q, p, t) = $ )) b4(q)p,p, + V(q, #). 

This shows us the form of any Hamiltonian function. 


We now wish to formulate Hamiltonian mechanics globally. To discover 
the correct approach, we compare two Hamiltonian systems 


H=T-V H=T-Y 
T =43Y b(q)p,p, T = $55") 3,8, 
V = V(q,t) V = Vq, t) 
gi= ees j= oH 
Op; Op; 
0H dH 
| | oq Bi = - og 


which are supposed to be defined on intersecting regions U, O and be 
equivalent on the common part U () O of U and U. The coordinate 
transformation 

(g, p, t) > (G, B, t) 


is given as in Section 10.1: 
q = q‘(q', are’ i 
aq! 


Bim L Pian 


10.2. HAMILTONIAN SYSTEMS 169 


We have 
bg Se ee 
» P; = q = » a ¢ 
0g' 0H aq' “s 
= La ip, op, =2 ag* 
Hence 
7 5 bg _ Og! ki 
= 4) 3B, P;= td Bizae P 
=) p,b"'p, = T, 
and so 
T=T 
on the common part of U and O, hence also 
H-H=TV7-YSJ. 


Using the symmetry relation, 


Op, , IPs 
qt * qi 


derived in the last section, and the remaining equations of motion: 
oH ap, Fi oP 
sa Le +e 
dp, 0H ie oH 
=~ Lag ap, Day ag 
q op, q' dq" 
0H 
agi’ 


6H _ OH 
ag! ror ag! 


= 0, 


From this, 


o 
af vee (t= 1,+°++, %). 


It follows that V — V is a function of ¢ alone, 
V=V+f(), 
H =H + f(t), 
this taking place on the intersection of U and 0. 
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Having this, we can formulate what we mean by a global Hamiltonian 
system. 

We begin with a position space M and form its derived spaces, phase 
space P, and state space S. We are give a function 7’ on P such that over 
any local coordinate neighborhood U on M with coordinates q’, ---, ¢" we 
have 


T =} b4(q)p,p,, 


a positive definite quadratic form. Here (g',-::,q",p1,°°*,D,) are the 
derived coordinates on the neighborhood U  E" of P lying over U. 

Now let U,, U,,--- denote the various local coordinate neighborhoods 
on M. For each one of these U, we have a function 


V,= V,(q, t) 


U,X E’. 
Whenever U, overlaps U,, then on the common part U, 1 U, we have 
Ve V5 = faalt), 
a function of t alone. (Clearly f,. + f;, = 0, and 
Sap t+ Spyt+Syn=9 on U,90U, 0 U,.) 
On the part of state space S lying over U,, ie., on U, X E" X E’, set 


defined on 


H,=T+V,. 
The equations of motion are given by 
4, _ oH 
Op, 
oH 
Li aq! 


on U, D4 E” >< E!, where (q') are local coordinates on U,. These are 
independent of the local coordinate systems (consistent) and define a motion 
(or flow) on all of S, always moving forward in the é, or time, direction. 


Remark. It is customary in mechanics to exhibit the potential function 
so we have set down the functions VY, in our formulation. Actually the 
equations of motion merely require the gradient of the potential which is an 
intrinsic quantity on S. Precisely, there is a differential form mw on S so 
that locally 





o=) 


OV ge 25 OV 
*dgi + ¥ —* dp,. 
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We can free ourselves altogether of reference to the functions V., by requiring 
that there be given a one-form m on S satisfying 


(1) @ is free of dt, 
(2) da = 0. 


By the converse of the Poincaré Lemma (Sections 3.6 and 3.7, especially the 
last remark on p. 31) this implies the existence of the functions V, . 

A trajectory of the motion is any particular solution of this system of 
differential equations. From the theory of ordinary differential equations 
we know that there is a unique trajectory through each point, so that state 
space S$ is smoothly filled with these curves. Along each one, ¢ steadily 
increases, but not necessarily to arbitrarily large values. (For example, a 
particle may run off to infinity in finite time.) 

Finally we note the energy law: 


Along any trajectory, 
aH = oH 
dt at 
For 
dH. oH oH 0H 
—=H=y— G+ ¥—9,+— 
a Lait tla e+ % 


_ _. Of 
= )(-b9) +) eb + = 


oH 
“Ot * 

Remark. Whether or not the functions f,,(¢) can be removed altogether 
by redefining each V, so that there will be a single potential function V on all 
of M X E? depends on two things, the manner in which the applied forces 
vary with time, and the topology of M. 

The topological difficulties are easily seen from the standard example of 
the steady magnetic field in the manifold M, which consists of E*> minus the 
z-axis, due to a steady electric current in the z-axis. The trouble is that 
there are closed loops in M which are not boundaries of surfaces. (Cf. the 
Situation in De Rham’s theorems, Section 5.9.) 


10.3. Integral-invariants 
Over a local coordinate neighborhood U =U, in position space M we 
consider the one-form 


w=, =) p,dq' — Hat. 
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This is defined on the portion U < E” x E' of state space which lies over 
U. These forms w, do not necessarily fit together to make a one-form on § 
because on an intersection U, 7 U, we have 


H, —H, =V,—Vy=fap(t)- 


If the V, can be chosen so that all f,,(¢) = 0, then we do have such a one-form 
on all of S. This is exactly the case for a globally conservative system, the 
case in which the external forces are derived from a single potential function 
V. While this cannot be expected in general, we do see that 


dw = dw, = dp,dq' — dHdt 


is a 2-form on all of S$, independent of local coordinates. This simply means 
that 
dw, = dW, 


on U, 1 U,, which is true because 
dH, dt = d(H, —f,,)dt = (dH, — f,,dt) dt = dH, dt. 


We shall call this 2-form dw, even though there is no one-form @ on all of $ 
which it is the “d”’ of. 

Suppose we have on some portion of $ an r-parameter family of solutions 
of the equations of motion. Let us denote by z,,---,z, the parameters. 
What this means is that we have a mapping ¢ on a region W of (t, x) space 
of the sort indicated, a cylinder in the ¢ direction with top and bottom 
curved r-chains, 


o: WS 
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where in local coordinates 
qgi=fit, x,, S52.) 
op: P = gilt, %y,°°*, 2,) 
b=, 


For each (z,), this represents a trajectory, hence 


of' OH 

— = — (f, ,t 
al ap, 6? B 4) 
O91 


oH 
= — —(f, g,#). 


The mapping ¢ is supposed smooth and one-to-one, so that an (7 + 1)-dimen- 
sional region in $ is evenly filled up by these trajectories. 
Now we compute ¢* (dw). We have 


dw = > dp,dq' — ES aaa — Lg, apd 
-r( oe tara) 


Op; 
Now 


#* (ag! - = at) = (F at + OF a) - Fat as, 


and similarly 


so that 


090") 2 
du) = 4D Fe sal ha 


= y A*(x, t) da ,dx ko 


which establishes our first point, ¢* (dw) 1s independent of dt. 
Since d(dw) = 0, we have d[¢* (dw)] =0. But 


a[p* (dw)) = dx, + (terms in dx,dx,dz,) = 0 











A 
Ot 
We conclude that 

oAs 


= Q, 
ot 
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each AJ* = AJ*(x) is independent of t, and we may write 

A differential form « of degree s on the state space § is called an (absolute) 
integral-invariant (historical terminology) if for each r-parameter family of 
trajectories given by such a mapping ¢, ¢* (a) is an s-form on the x-space 
alone (no ¢ nor df terms) and if in addition dx = 0. 


Each of the forms 
dw, (dw)*, +++, (dw)" 


is an integral-invariant. For 
d(dw)* = 0 
and 
p* (deo)® = [p* (deo) }° = [) AM (x) dx, dx, } 


is independent of ¢ and df. 





Consider in state space S a small piece c, of surface which is filled by a 
one-parameter family of trajectories. We may describe c, analytically by 


f =q'(t, y) 
Pp; = pilt, y) 
aly) St S oly) 
| Yo SYSN- 


Our reasoning above shows that the two-form we get by substituting these 
expressions for p and q in dw is a two-form in y and dy only, hence vanishes. 
This means in particular 
| dw = 0. 
€2 
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Next, suppose one has a piece of volume, or three-chain c,; in § which is 
the span of a two-parameter family of trajectories: 


i = qitt, Y; 2) 
i= pilt, y; z) 
= z) StS bly, z) 


(y, z) in a domain D. 


P 


x, 


Then 
0c, =X, —-2%y+c¢, 


where £,, Lp are the terminal and initial surfaces respectively and where the 
lateral surface ¢, is spanned by a one-parameter family of trajectories 
corresponding to the parameter point (y,z) on dD. Now 


| dw -| d(dw) = 0, 
OC; €3 
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| dw = 0, 

€2 

| dw -| dw. 
ry Xo 


We have established the very striking property of the form dw: 
If Xo 13 any 2-chain in § transversal to the trajectories, and if one displaces 
each point of Xp any amount along its trajectory to form a new surface XL, , then 


[> 


It should be clear that we have not used any special properties of dw in 
proving this result so that any integral-invariant satisfies a corresponding 
property. 

Let a be an integral-invariant of degree r on S. Let ¢, be any r-chain on § 
transversal to the trajectories. Let c', be a second such r-chain so that the points 
of c, and c’, may be put into one-one correspondence with corresponding points 
on the same trajectory. Then 

[a= [+ 
Cc, C’r 


It is possible to reverse our steps to prove that the property expressed in 
this result actually characterizes integral-invariants. This is done in 
Cartan [8]. 

We pass on to relative integral-invariants. 

An r-form « on § is a relative integral-invariant provided da is an integral- 
invariant. The basic result about relative integral-invariants is this. 


and we showed above that 


hence 


Pp 
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Let « be a relative integral-invariant of degree r. Let b, and b’, be two 
r-dimenstonal boundaries which are in one-one correspondence in such a way 
that corresponding points are on the same trajectory. Then 


fe-|e 
r b’, 


To prove this, we select (r + 1)-chains ¢,,,, ¢’,4, 80 that 
b, = OC, 41 ? b’, 7 OC'n+1 


and do this in such a way that c,,, and c’,,, correspond one-one with 
corresponding points on the same trajectory. Then 


[e-| e=| ax = | da = | a~ [ix 
b’, BC’ p+ 1 C’r+ 1 Cy+34 OCe + i b, 


The differential form 
oO= > Pp aq’ —-H dt, 


which is defined locally, is a relative integral-invariant of degree one where 
it is defined. Consequently so are the forms 


wdw, w(dw)*, +--, w(da)" 


since 


d[w(dw)"] = (dw)"*". 


We shall specialize by considering chains which exist at a single instance 
of time. Let « be any differential r-form and ¢ an r-chain in S lying in a 
hyperplane ¢ = constant. Then clearly 


[2 = {lane 
ce ce 
P 
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Let us apply this to dw in particular. We have 


do) | ge=0 = » dp, dq’, 
which leads to this result: 


Let X, be a 2-chain in S$ at t = ty and X, the 2-chain obtained by moving each 
point of Xp along its trajectory totimet=t,. Then 


In this result we may think of LX, and Z, as 2-chains in phase space. 
If we apply the procedure to the (27)-form 


(dw)" = + ni(dp, ---dp,dq* «++ dq") + Adt 
which gives us the phase density 


p= dp, ++: dp,dq' --- dq’ 
we have the 


LiouviLLE THEOREM. If a 2n-dimensional region Dy in phase space at 
time ty moves to a region D, at time t, , then 


Jot Jat 


We shall close this section with a result which will be needed in Section 
10.5. It shows that the integral-invariant dw completely determines the 
equations of motion, in itself an important mechanical principle. 

Let 

* = A'(t, q, P) 
PB, = Bit, q, P) 
be a system of equations on a region of state space S which has 
dw =) dp,dq' — dHdt 
as an integral-invariant. Then 
H H 
_@ B 0 


— 


A'= op, a 
For let 
i = f(t, %1, °° +, ay) 
Pi = Git, Xy, °°", Taq) 
be a general solution, so that dw must be expressible in the terms of dz, 
alone. From this, and the differential equations 


of' i 09; 
a 
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we deduce 
dq’ = A‘dt + A’, dp,= B,dt+ pn,, 


where A’, 4; are one-forms in the dz, alone, and so 


=) dp,dq' — dH dt 
oH CH 
= (L Alm — Balla Dad! — (La + DS a) at 
Since dw is free of dt, 
CH oH 
Aine BAH YY —7'= = 0, 
y ui > i Dag D Gp, Ky 
CH oH 
AP = SS B —)ji'=0 
( 5M z ( +55); : 
oH 
E(4 ' - =) dp. - B,dt) - -¥ (B+ )tda! - Atay) = 0, 


oH oH , 
E (4 — =) dps — (Bs + $a) da + (yd = 


We conclude that 


oH 0H 


as asserted. 


10.4. Brackets 

In the transformation theory of classical mechanics one uses the bracket 
expressions of Poisson and Lagrange. In this section we shall show how 
these expressions relate to differential forms. 

Before doing this, it is a good idea to digress on the subject of Lie brackets. 
We take any differentiable manifold N and recall the definition in Section 
5.3 of tangent vector and the definition of vector field at the end of that 
section. We may consider a vector field v on N as an operator which takes 
each function on N to another function on N: 


v: F°(N) — F°(N). 
If z!, ---, 2" is a local coordinate system in which v has the representation 
7) 
-_ i ads 
v=) a'(x) aqh? 
then 


vin = Tain & 
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shows, locally, what v does to a function f. One sees from this formula, or 
from the very definition of v as the assignment of a tangent vector (directional 
differentiation) at P to each point P of N, that 


v(f-9) = v(f)-9 + f-v(9) 


for any two functions f and gon N. If v and w are two vector fields on N, 
we define the Lne bracket of v and w by 


[v, Ww] =Vow— wo. 
This is another vector field on N. If in local coordinates 
veyac, wad, 
then a computation shows that 
b al 
(V.w]=> (ras — po) —_ 
(The main point to notice is that the second partials cancel each other.) 
The following algebraic identities are easily established: 
(i) [v, v] = 0. 
(ii) [w, v] +[v, w] = 0. 
(iii) [v, + v2, Ww] =[v,, Ww] +[v2, w]. 
(iv) [[u, v], w] +[f[v, w],u]+[{w, u],v]=0 (Jacobi identity). 


Now we return to mechanical systems. As before, let P be the phase 
space associated to a position space M. We denote by « the differential form 


a= > p,dq' 
on P so that 
da = >} dp,dq' 


as in Section 10.1. 


Poisson Brackets. To each pair f, g of real functions on phase space P 
we associate a new function (f, g) defined by 


n(df dg)  (da)"~* = (f, g)(da)”. 
In local coordinates 
a 
(da)"~* = [(n — 1)!] }) (dp, dg") -- + (dp, dq’) -- + (dp, dq") 


and. 
(da)" = (n!)dp, dq’ --- dp, dq", 
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from which we deduce the local expression for (f, g): 





a(f, 9) 
(f,9)= 
J 2 O(p i: q') 
From the definition one derives these relations: 


(i) (7, f) =0. 
(i) (f,9) +9, f) = 09. 
(ili) (f,91 +92) = (f, 91) + (Ff, 92). 
Using (af) d(9192) = 91 (dfdg2) + 92(dfdg,), one has 
(iv) (f, 9192) = 91° (fF 92) + 92° (fs 91). 


The identities (iii) and (iv) taken together may be expressed by saying that 
for fixed f, 


is a vector field Vv, on P: 


vs(9) = (f, 9). 
The basic connection between the Lie and Poisson brackets is given by 
(V) Vy, 9) = [Ye » VYo]- 


One has 
Vor, 9)() = (fF, 9); h) 


[v,, v,](4) = v,[v,()] — v,[v,(A)] 
= vy((9, 2) — v4((f, 4) 
= (f, (9, 2)) — (9 (Ff, )) 
— ((9,4),f) -(a Sf), 9) 
so that (v) is equivalent to Jacobi’s relation 
(vi) ((f,9), 4) + (9,4), f) + (Af), 9) = 0. 


One proves this relation [(v) or (vi)] either by a lengthy direct calculation, 
or by the following more sophisticated argument based on the fact that a 
vector field is completely determined locally by its effect on each member of a 
local coordinate system. 


and 


ll 





First of all, 
i, ph) ol oF of 
0) = Lae eh Xap, aa, 
0 5 


of 


i 
(f,¢') Bp,’ 
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which may be interpreted as 


0 
Vi ap! ; 
Similarly 
of 7] 
Pi=-— d =—. 
(f P,) dq' an Voi aq' 


Because of these relations (vi) easily follows when both g and h are taken from 
the set of coordinate functions {q’,---,p,}. This means that 


Vos,x) = [v, Vx] 


when z is one of these coordinate functions since the vector fields on both 
sides agree when applied to any coordinate function q/ or p, Hence for 
any h, 

((f, 2), h) + ((z, h), f) + ((h, f), x) =). 
This is now established for any functions f and h and any coordinate function 
x. But this may be interpreted as saying 


Von, ryt = (Va» Vy le 
so that the vector fields 
Vu. and [Vn ¥y] 
must agree since they agree on all of the coordinate functions. Hence 
Vin, p)9) = (Ys Ve lg 


for all f,g,. This completes the proof of (v) and (vi). 

One applies brackets to a pair of functions on the state space S by simply 
treating ¢ as a parameter. 

A function f on $ is called a first integral of the equations of motion if it 
is constant on each trajectory. If this is the case, then 


_ af Le 2 
mae +e 

of aH —, of oH af 
= 2 iat ap, -2 ap, oq! * at 
=(,f) +2 


so that the partial differential equation for first integrals is 


of, H), 


10.5. CONTACT TRANSFORMATIONS 183 
LAGRANGE Brackets. Let 
@: E*—P, 
where E* is the Euclidean plane with rectangular coordinates u,v. Then 
o* (da) is a 2-form on E? (where as before da = )' dp,dq') and we write 
b* (da) = [u, v]dudv 
defining the Lagrange brackets which have the local expression 
y 0(P;, 7’) 


u,v) = : 
lu, 2] O(u, v) 
10.5. Contact Transformations 

Here we can but touch briefly on an extensive topic. For simplicity we 
shall only treat the local problem and shall look on contact transformations 
as coordinate changes. 

First we take the case of a time independent change. As usual, let q', p, 
denote local coordinates in phase space P. We consider new coordinates 
q'; p i? ; A 

q = 97'(q Macias (28 Se: De) 
P; = pq’, "49", D157") Dy) 
which are unrelated to the old except for one requirement. We set 


a=) p,dq, a=) p,q 


da = da. 


and require 


This defines what we shall call a homogeneous contact transformation. 
Since d(a -- «) = 0 and we are only working locally, the condition may be 
expressed as 
g=a+dd 


where ¢ is a real function on P. In the relation da = dx, 
. dp, dq = >, dp, dq’, 


one replaces dp, and dg' by their expressions in terms of dp, and dq* to obtain 








OB; 7’) 6(B;, 7) ;, 7’) 
: = Q, — ee 0, = 1}. 
(5 a(q’, 1") y OD; > Px) d Ap;g)  * 
(These relations may be expressed in terms of Lagrange brackets.) 
Similarly, 
_ ogi Sy . oF =| 
—— +=, Pa ee 
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If we have equations of motion 


oH 
ap, 
0H 
q' 


oi 


p= 


we shall show that they transform into equations of the same type. For 
suppose after changing the coordinates according to our contact transforma- 
tion we arrive at 


g! — A! 
B; pe ee B; 
so that 
i_ y OY OH -y oy’ OH 
dq) dp; ~ ap, aq! 
0p, OH 0p, oH 
= Diag éq/ dp, +2. Op, agi” 
If we multiply the first by 


Op Op 
dp,= > ag Ot + 7a Pap 


and the second by dq similarly and sum, we find after some simplification 


_ : oH oH 
> Aldp, + }) Bi ag a2 At te on, dp, 
oH 
H —— dt. 
ot 
Hence 
oH oH 
Ai=—, B = —— 
Op; og’ 


so that the equations of motion in the new coordinates are precisely 


0H 
ap, 
oH 


i= — sa 


gi= 
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Now we pass to the general situation. We begin with a mechanical 
system 
_ OH 
Op, 
oH 
ag 


. 
tt 


Pi= 
on the state space § and consider a coordinate change on §: 


git, qty. qd", Dis ' +> Da) 
P; Bilt, g's 0°59", P15 ° "Prd 
t. 


ba] 
ll 


We set 
O= > p,dq' = H dt 

as usual. The coordinate change is called a contact transformation if there is a 
function H and a function ¢ so that if 

O= YB dg’ = H dt, 
then 

@=a-+ dd, 
or what is the same thing, 
dw = dw. 

The first basic result is that the equations of motion in the new coordinates 
are precisely 


, oH 

1 a, 
aH 

Pix — FR 


For whatever the new equations of motion are, they admit 

da = ¥ dp,dg' — dHdt 
as an integral-invariant. But the final result in Section 10.3 tells us that 
the only equations with this integral-invariant are the stated ones. (Com- 


pare this slick proof with a direct computation!) 
A particular type of contact transformation is obtained as follows. Let 


Pt,z2,°-- ayer y’) 
be a function of 2n + 1 variables satisfying the independence condition 
a°o 


Ga! Ay! ne 
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Settd=dg° °F. g), 


_ op . op 
Pi agi? Pi= og’ 


For fixed g and p, the g are determined by the first set of equations (since 
the determinant above does not vanish) and then the j are determined by 
the second set of equations. 


We have 
i op =i 
iia dt + 5 aa + a agi 4 
— Cx i = i 
~ Be Yi p,dg' — > H,47', 
wo — @ = —(Y p,q’ — Hat) + (¥ p,dq' — Hadt) 
= +49 Pars Hat nat 
and so 
Go=a-—dd 
provided we set 
H= H+o. 


The most important case is that in which ¢ is a solution of the Hamilton- 
Jacobi equation 
op 


op 
t, : -= H t, > eee > ee one — 0. 
ae | q,9) + ( q qi 


In this situation H = 0 and the new equations of motion are simply 
B, = 9 


with solutions 
g ‘ = constant, p,; = constant. 


The original system is said to be transformed to equilibrium. 

One other point we shall notice is this. If a contact transformation is 
stationary, i.e., independent of time, then it is equivalent to a homogeneous 
contact transformation. For suppose 


g' = 7'(q, P) 


Pi = Bq; P) 
and 


(¥ #,dq' — Ht) = (¥ p,dq' — Hat) + dg. 
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Equating coefficients: 


_ og op 

LPs = Pit 3G 
. oY d¢ 
dP Cp, - Op, 

Ho-n-@ 


Ot! 
Since q, p, J, p are independent of t, we deduce from the first two equations 
that 
2 2 
Cc ey 
dt dg! at op, 
Hence 0¢/é¢t is a function of ¢ alone. This evidently implies 


P(t, 9, p) = f(t) + Wg, P) 


and so 
» B:4q' = ¥ p,dq' + dy 
which means we have a homogeneous contact transformation as asserted, 
Let us briefly examine what are called infinitesimal contact transformations. 

We ask when 

g=q' t+ ef! 

Pi = Di + 89; 
is a contact transformation up to first order termsiné. We restrict attention 
to the homogeneous (stationary) case. We have 


> p,q — ¥ p,dq' = edd, 
> (p, + 9,)(dg' + edf') — } p,dq' = edd, 
e() 9:4q' +) pdf’) = ed@, 


so the condition is 

¥ 9:dq' + ¥ pdf! = dd. 
If we set 

yp=o- Y pif! 

this becomes 

»» g, dq’ = > fidp; = dy, 
or 
ow oy 


a qi’ Pea 
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We finally have 


op 
=! i 
=| =f 
Op; 
4 
Di Pr a 
oq 


Here yw is called a generating function for the infinitesimal contact trans- 
formation. 


10.6. Fluid Mechanics 
We consider a fluid moving in a region of E*. The position vector as 
usual is 


x = (x, y, z) = (z!, x, v3). 
At each time ¢, the velocity at x is v, 
v = v(t, x) = (u, v, w) = (v!, v?, v?). 
The density of the fluid is a scalar 
p = plt, X). 
In this section we shall denote the vectorial area element of a surface by a, 
o = (dydz, dzdz, dxdy). 


(See Section 4.5, p. 43.) 
If c, is a three-dimensional region which is fixed in space, the change in 


mass at each point x of ¢, per unit time is 


- dx dy dz 


and so the total time derivative of mass in Cc, is 


op 
— da dy dz. 
|, Zardvae 


We assume conservation of matter, so this must result from flux of fluid over 


the boundary, hence 
0 
| oP da dy dz = -| pv-a. 
C3 ot C3 


By Gauss’ theorem, 


| pra =| div (pv) dx dy dz. 
oc, C34 
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By taking c, arbitrary we deduce from the equality of these integrals the 
continuity equation 


dp 
rr — + div (pv) = 


a necessary condition the flow must satisfy. We shall deduce some con- 
sequences of this. We set 


Q = p(dx' — v' dt)(dx? — v? dt)(dx* — v° dt). 
To compute dQ, we set 


B = (dz! — v! dt)(dx* — v? dt)(dx* — v° dt) 
so that 


3 
dp =~ ade da! dtdx* dx? + dz! (ss dx? at) dx? — da! dx? a dx? dt 
dx} Ox? 6x3 
= (div v) dt dx’ da? dx’, 
dQ = d(pB) = dp a B+ pap 


-(2 dt + a, ep ; dx’ A B + p(div v)(dt da! dx? dx) 


-| + > v' — + p(div »)| (dt da’ dx? dx). 
Thus the  leiae equation is equivalent to the relation 
dQ =0. 


Suppose we express the flow in terms of initial conditions (or other 
parameters) by 
x= x(t, a!,--+, aw) 


so that the a! are the parameters and 


Thus 


o 
(da? — v'dt) = (= dt + >a) — v' dt 
a 


so that 
O(x!, x”, x3) 
A(a', «7, «3) 


= A(t, ada! da? da. 


Q=p da) da* dat 
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Since dQ = 0, we deduce that ¢0A/dt = 0, 
Q = A(a) da! da? da?. 
This means that © is an integral-invariant for the flow as explained in 
Section 10.3. We consequently have the following result. 
Let ¢,, ¢', be three-chains in the four-dimensional (t, x) space which are in 
one-one correspondence in such a way that corresponding points lie on the same 


trajectory of the flow. Then 
| Q- | 0. 
C3 c's 


In particular, if all the points of c, exist simultaneously, i.e., at a fixed 


time f,, then 
| a= | Ql =19 -| pdx dy dz. 
3 €3 C3 


Thus if we take a region c;°°) at time ¢, and follow it to c," at time ¢, , we 


have 
| pardy as = | pdx dy dz 
369) e3()) 


which says that mass is preserved in the flow, another form of the conserva- 
tion of mass. 

We now proceed to the dynamic situation. We suppose our fluid is 
nonviscous so that the pressure is a force per unit area at each point acting 
normal to any surface element through the point, always with the same 
magnitude. Let 

p = p(t, x) = pressure 
F = F(t, x) = body force per unit mass. 

Let c, be a fixed region in space. The total acceleration of all matter in 

C; is 


dv 
p — dxdy dz. 
ea 


At the instant of time in question, this must equal the total force on the 


matter in ¢; which is 
| pFdzdydz -| pa. 
C3 ec3 


By a variantt of Gauss’ theorem, 
| pa -| (grad p) dx dy dz, 
0€3 c3 


t The usual simple proof is to dot both sides into a constant vector a and apply Gauss’ 
theorem plus the fact that 
div (pa) = (grad p)-a. 
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hence 


dv 
| (> Ti pF + gradp) dzdydz = 0. 
c; 


We conclude that 


dv fe 1 q 
di raid P, 


the Huler equation of motion. 
Here the interpretation is 
dy ov n ov, 
dt ot Saat 
Let us suppose that the body force F is conservative, 
F = — grad V 
where 
V = Vit, x) 
is the force potential. 
We shall add the hypothesis that p and p are functionally related, i-e., 


dp Adp=0, 


as is the case, for example, with an isothermal motion. In this case we 
can define a function g = q(t, x) by 


(t, x) dp 
i-("2 


0 p 
so that 
dp 
dq =. 
p 
The equations of motion may be written 
d 
<= —gred(V + 9) 
or 
du 0 
—=—-—(V te., 
1.€., 
Ou Ou Cu Ou 0 
ee uae a ae a ee : te. 
Ot og ag ee 3a | +9) iad 
We set 


E=}(v-v)+V+q 
=f(w +07 4+w7)+V +49, 
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the energy per unit mass. Finally we define the vorticity 
Ow dv du dw dv du 
r=(nO=(S-2 S-S 2-S). 
Oy oz dz Ox Ox Gy 
We compute 
ol LR eae (oS + w 2) jig ie 0 (V 44) 
dz at du \ dx On) | Ox * 
= (oS + w 2) (oS + w} 
 N ée Ox oy Gz 
= vf — wy 
and similarly 
OE % Ou vt 
gx Ot ee 
OE av 
ey + a we — ul 
o£ n Ow P ve 
ar 
(These are equivalent to the vector formula 
Ov 
dzE+4+—= , 
grad E + ry VX.) 
Now we consider the differential form 
w= udx + vdy + wdz — Edt. 
We have 
dw = (Edydz + ndzdx + Cdxdy) 
+dt(u,dx + v,dy + w,dz) 
—(H,,dx + E,dy + E,dz) dt, 


dw = (Edydz + ndzdx + C dxdy) 


+ dt[(v¢ — wn) da + (we — ul) dy + (un — v6) dz] 


Actually, dw is an integral-invariant so that q is a relative integral-invariant 
One sees this indirectly by making a comparison to Hamiltonian systems 
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For if one thinks for the moment of z, y, z, u, v, w, t as independent variables, 
the equations of motion are 


; 0 
L=U Be eg) 

: 7) 
y= v= —>(V +9) 

oy 
zZ=w is pa 

Oz 

or 

/(%= 0E/ou u= —0E/Ox 
y = OE /dv o = —d0E/oy 


z= 0E/dw w= —dH/dz 
which makes our assertion clear. It follows that if c,, ¢’, are two-chains in 
(t, x) space which are in one-one correspondence so that corresponding points 
are on the same trajectory, then 


| dw -{ dw. 
c2 e’2 


In particular if ¢,°° is a 2-chain in E? at time ty and it moves to ¢, at 
time t, according to the motion, then 


| (Edydz + ndzdx + Cdxdy) = | (Edydz + ndzdx + Cdady). 
c 20) 


29) 
This is the Helmholtz theorem on conservation of vorticity. In the first 
integral we must understand € = €(ig, x), ete., and in the second, € = E(¢,; , x), 
etc. Animportant consequence is this further result of Helmholtz. Suppose 
at a fixed time ty, the vorticity vanishes identically. Then tt always vanishes. 
For by the formula above, 


| (Edydz+-:-)=0 
e2()) 


for each 2-chain c,“) at each time ¢,, which evidently means that the 
integrand must vanish. 


10.7. Problems 
1. On p. 54 and again on pp. 179-180 we defined a vector field v on a 
manifold M as a mapping which takes each function f on M to another 
function v(f) and satisfies the rules 
Vf + 9) = v(f) + v(9) 
v(f9) =f-v(g) + v(f)-9- 
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We know that locally 


7) 
v = Laz) =. 


We have alluded several times to the duality between vector fields and one- 
forms. (Cf. pp. 127, 143 for example.) Now we shall make this more 
explicit by associating with each vector field v and each one-form o on M 
a scalar (a, v). Locally, if v has the representation above and 
o =) bx) dz’, 
then we set 
(a, v) =} a'(x)b,(2). 
Show that 
(0, + 02, Vv) = (6,, ¥) + (62, Y), 


(0, ¥, + V2) = (0, V,) + (2, V2), 
(go, V) =9°(0, V) = (0, gv) 


for each scalar g. Show also that this operation has the intrinsic characteri- 
zation 


(dg, v) = v(g). 


2. Show similarly that two-forms may be paired with two-vectors by 


(0,,V;) (64, V2) 
(0; AG2,¥, A V2) = . 
(02,431) (02,2) 


3. Now prove the formula 
(dco, v A w) + (a, [v, w]) = v{(o, w)} — wi{(a, v)}, 


which relates all these things to the Lie bracket. 

4. Combine this with the Frobenius integration theorem (Section 7.3) to 
establish this result: Let v,,---,v, be vector fields on a neighborhood of 
0 in E” which are linearly independent at each point. Suppose that for each 
t and j, [v,;, V,] is a linear combination of v,,:--,v,. Then there is a 
coordinate system z', --- , x" of some neighborhood of 0 such that 


r ; CO ; 
v; = Laie) et (@=1,---,r). 


5. Letv,,-:°:, Vv, be vector fields on a neighborhood of 0 in E” which are 
linearly independent at each point. Suppose that each bracket [v,, v,] 
vanishes. Prove that there exists a local coordinate system z',:--, 2" 
such that 
0 


“i= Ox! 


10.7. PROBLEMS 195 
6. Let u,,-°:,u, be functions of z!, +--+, x" and set 


the Jacobian matrix. Write J = A + B where A is symmetric and B is 
skew so that A = |la,,||, B = ||6,,|l 


1/du; du; l/du,; du, 
Qi; = (a0 + =), b;; = 5(= ~_ =): 


(i) Suppose B=0. Prove there is a function f satisfying df =) u,dz'. 
(ii) Suppose A =0. Prove that J is constant so that the u, are linear 
functions of z!,---,2". (This comes up in the theory of strain.) 
7. Let a;; be n* functions of z',---,2". Show that the integrability 
conditions 
07a,, 


Ox* Gx! = Aart Ax! ~~ Ax* dx! = Ax! Ox! 


are necessary and sufficient for the existence of functions u,,---, u, 
satisfying 
l/du, du, 
ass =slagt a): 
4 2\ae/ Ga! 


Investigate n = 3 and =2. (Due to Saint-Venant. See Love [18, p. 49].) 
8. Let v be a vector field in E*. Show that there exists a function ¢ and 
a vector field w with div (w) = 0 so that 


v = grad (%) + curl (w). 


(Hint: use the mechanism of harmonic differentials. Assume that if g is 
any function in E*, the equation Lap(f) = g has a solution. See Abraham- 
Becker [1, pp. 37-38] and Love [18, p. 47].) 

9. Consider the transformation 


Gg =F'(q', oe 5Q") 


Bi= YP 5. 
i > 4 ag! 


The problem is to give a new proof of the relation 


Op; , OP, 
dg" + ag: 
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derived in Section 10.1. Set 


og 
and prove 
(i) YA, A dq =0 
(ii) A; =) on 49" 
from which you conclude that g;, = 9,;. Show that this is equivalent to the 
desired relations. 
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Glossary of Notation 


This is a list of special terminology, definitions, and symbols, especially 
from the early chapters. 


A. Spaces 


E” Euclidean n-space. 
R The set of real numbers, also considered as E', the Euclidean line. 


U,V,-:- Open sets (in E”, or on a manifold). 
L,M,--- Vector spaces. 

PL The space of p-vectors on L. 
M,N,::: Manifolds. 


F?(U) The collection of all p-forms on U. 
Cartesian product. If S and T are arbitrary sets (collections of objects), 
their cartesian product is the set 


sx T 
consisting of all ordered pairs (s, t) where s belongs to $ and 


t to T. 
7S This is the cartesian product $  S of S with itself. Similarly 


K?*S=S xs KS, ete. 
SMT _ This is the intersection of the sets $ and T. For example, if 
S = {1, 2, 5, 7} and T = {2, 3, 7, 9}, then Sr T = {2, 7}. 
[| The unit interval 0 <i < 1. 


B. Functions 

Mapping. A mapping is a smooth function ¢ from one space M to another 

N. We write 
od: M—N. 

Composite mapping. If ¢: M—>N and wy: N-—-P, then we may 
form the composite mapping w.od: M—P. It is defined 
by 

(yo P)(x) = plp(z)) for x in M. 
See pp. 3, 22, ff. 
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Linear functional. 


GLOSSARY OF NOTATION 


A linear transformation on a linear space L to the one- 
dimensional space R of real numbers. 


Jacobian. Ifu!’ = u‘(x',---, 2")(i=1,--+, 2), the Jacobian of this mapping 


o* 
Px 


‘A 


[Aj 
dim L 


ot 


lsh, <h,< ree CAS, 


H' 


sgn 7 


* 


(x, B) 
a 


YN 
dt} n+++andzia 


is the determinant 
| Ou! /da? |. 


The mapping on differential forms induced by the mapping @¢ 
between spaces, p. 23, ff. 

The mapping on chains induced by the mapping ¢ between 
spaces, p. 71, ff. 


C. Special symbols 


The transpose of the matrix A, obtained from A by inter- 
changing rows and columns. 

The determinant of the linear transformation (matrix) A, p. 21, ff. 

The dimension of the linear space L. 

Here H = {h,,---,h,}, a set of indices in increasing order, 


ota ot aga: ag, 


This is the complementary set of indices. For example, if 
n = 8 and H = {2, 3, 5, 6}, then H’ = {1, 4, 7, 8}. 

If x is a permutation on {1, 2,---,}, then sgnz=1 if x is 
effected by an even number of interchanges (of two numbers) 
and sgnx = —1 if z is effected by an odd number of inter- 
changes. 

The star operator, p. 15. 

The inner product, p. 12, ff. 

The boundary operator, p. 58. 


and 


i+1 


+++ A da” means dx} A+++ Adz’ a daitl ja cee a dat. 


The circumflex indicates omission. 
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Cycle, 63 
Cylinder construction, 27 ff 


D 


Debever, R., 198 

Degree of a mapping, 77 ff 
De Rham, G., 197 

De Rham’s theorems, 66 ff 
Dirichlet integral, 83 


E 


Electromagnetic field, 16, 44 ff 
Euler equation of motion, 191 
Exact form, 67 

Exterior derivative, 20 ff 
Exterior multiplication, 8 ff 


F 


Faraday’s law, 45 

First integral, 107, 182 

Flanders, H., 198 

Flat Riemannian manifold, 135 
Fluid mechanics, 188 ff 

Flux, 43 

Frame manifold, 146 

Frenet formulas, 122 

Frobenius, G., 4 

Frobenius integration theorem, 92 ff 


G 


Gallissot, F., 198 

Gauss, C. F., 43 

Gauss integral, 79 

Gauss mean value theorem, 85 
Gaussian curvature, 42, 118, 126 
General linear group, 156 
Geodesic, 119, 134 

Goldberg, S. I., 198 

Golomb, M., 198 


204 SUBJECT INDEX 


Goursat, E., 4, 26, 198 

Grammian (Gram determinant), 12 
Green’s formulas, 83 

Green’s function, 86 


H 


Hadamard, J., 198 

Hamilton-Jacobi equation, 186 

Hamiltonian, 166 

Hamiltonian systems, 165 ff 

Hamiltonian equations of motion, 167 

Heat equation, 90 ff 

Helmholtz law of 
vorticity, 193 

Hodge, W. V. D., 15, 136, 138, 198 

Homogeneous contact transformation, 
183 

Hopf, H., 197 

Hopf invariant, 79 

Hypersurfaces, 116 ff 


conservation of 


I 


Infinitesimal contact transformation, 187 
Infinitesimal transformation, 37 
Integral-invariant, 174 

Integral surfaces, 92 

Injection (imbedding), 53 


J 


Jacobi identity (relation), 180, 181 
Jordan-Brouwer theorem, 77 


K 


Kahler, E., 4 

Kinetic potential, 166 
Klein bottle, 71 
Kronecker integral, 77 


L 


Lagrange brackets, 183 
Lagrangian, 166 

Lamb, H., 198 

Laplace expansion (determinant), 10 
Laplacian, 38 ff, 44, 82 

Last multiplier, 107 

Left invariant, form, 150 
Left translation, 150 

Lie, 8., third theorem of, 108 
Lie bracket, 180 

Lie group, 150 ff 


Line of curvature, 120 

Linking number, 79 ff 

Liouville theorem (harmonic functions), 
88 

Liouville theorem (phase space), 178 

Local coordinate neighborhood, 49 ff 

Lorentz inner product (or metric), 12. 
46 

Love, A. E. H., 195, 198 


M 
Manifold, 49 ff 
Maximum principle, 85 
Maxwell’s equations, 16, 44 ff 
Mayer, A., system of, 104 
Mean curvature, 42, 118, 126 
Minimal surface, 44 
Misner, C. W., 198 
Momentum components, 166 
Monge notation, 126 
Moving frames, 4, 32 ff 


N 
Nickerson, H. K., 198 


0 


Ordinary differential equations, 106 ff 
Orientable manifold, 51 

Orientation (linear space), 15 
Orthogonal coordinates, 39 
Orthogonal group, 158 

Orthonormal basis, 13 ff 


P 


Parabolic equation, 90 

Parallel displacement, 119, 133, 145 
Parallel surface, 113 

Periods, 66 ff 

Phase density, 10, 164 

Phase space, 163 

Poincaré, H., 4 

Poincaré lemma, 2, 4, 27 ff 
Poincaré metric (half-plane), 134 
Poisson brackets, 180 

Poisson integral formula, 87 
Position space, 163 

Potential theory, 82 ff 

Poynting vector, 47 

Principal curvatures, 42, 120 
Principal directions (hypersurfaces), 120 
Projective plane, 70 


Projective space, 73 

Proper affine group, 154 
Proper orthogonal group, 158 
p-vector, 5 ff 


R 


Relative integral-invariant, 176 
Riemann curvature tensor, 122, 131 
Riemannian geometry, 4, 127 ff 
Right invariant form, 151, 159 
Right translation, 151 


S 


Saint-Venant, B. de, 195 
Seifert, H., 77, 198 
Signature (inner product), 13 
Simplex, simplices, 57 ff 
Sommerfeld, A. J. W., 198 
Special linear group, 157 
Spencer, D. C., 198 

Spherical coordinates, 34 
Spherical harmonies, 142 
Standard n-simplex, 60 

Star operator, 15 ff 

State space, 165 

Steenrod, N., 198 

Step transformation group, 155 


SUBJECT INDEX 205 


Stokes’ theorem, 2, 64 ff 

Submanifold, 52 

Support function, 114 

Surfaces, differential geometry of, 40 ff, 
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